testlist <- list(p = structure(c(1.30613572389549e-300, 3.13250230186967e+301, 4.14467815817155e-52, 7.44454355882871e-262, 4.43429504733332e-31, 3.37114523626541e+103, 3.74883401087728e-119, 7.81100138067908e-10, 1.91918387411323e+207, 6.28509003523093e-108, 4.57527140520858e+84, 5.7975105097069e+140, 1.71850040668305e+176, 3.66533613049169e+264, 5.7440465360135e-207, 4.77042141991774e-277, 1.38138443354758e-255, 6.38527321716769e-297, 7.78489137915683e-32, 4.78626638744197e-49, 1.99268053418731e-224, 2.9321077422219e+284, 4.01374126919666e-196, 4.55769487761222e-183, 5.6306742796553e+135, 1.47251938414639e+37, 5.30202068066352e-119, 8.77362320842782e+216, 9.57460156214806e+41, 1.87494576589683e+23, 4.29962994621425e-190, 2.53197754195931e-15, 1.78367782816586e+304, 1.52158613637624e+81, 3.61657884101299e+81, 4.26546181368945e-120, 7.95111840517066e-116, 9.8574313890341e+140, 1.48259426103242e-205, 1.5727114142785e-193, 7.87207510185534e-272, 1.69241397174968e+246), .Dim = 6:7), pi_inv = structure(c(8.35983597255897e+183, 12433685310.5968, 1.88458982056426e-83, 2.23994031531595e+229, 3.19721198215496e+282, 1.92928173328364e+130, 9.98058235350412e-122, 2.87578537616475e-19), .Dim = c(1L, 8L)), ref_order = structure(c(5.27068094429231e-37, 1.64521889151021e-68, 8.27875944152203e-70, 8.88429034480326e-268, 2.39811641606243e-190, 1.7572694841654e-172, 2.63185615575886e+217, 5.38375671452751e-115, 5.42583662088686e+306, 0.000487701596956602, 1.49910982461869e+106, 6.18633246976901e-180, 3.68391043901549e+107, 1.21282005443679e-109, 1.12829454704595e-188, 7.76622176989471e+278, 3.40388471093251e+253, 2.24873694164131e-87, 5.82853312361407e+222, 4.48981346078649e-41, 1.3496300135629e+293, 6.48395768149842e+180, 7.19937540315476e+287, 1.79760659115635e-148, 8.26737095193579e-29, 7.81967151890548e+132, 1.05252849016904e+245, 1.54863440815549e+175, 1.18087509744805e+229, 2.56314636567102e-177, 3.73982361828618e+53, 3.34575651062865e-302, 4.14538150793955e-129, 1.05083428201464e+232, 2.15324492596084e-43, 1.7173166825214e-74, 1.00693318697063e+267, 1.45294014331008e-08, 1.04615380650751e-161, 5.72971392802784e-145, 4.86363037161214e+154, 6.26973104679236e+31, 3.81591116500471e-204, 1.54776947079393e-257, 9.37364909229576e+154, 3.07919623867538e-24, 1.88673723951606e-58, 4.0627703243747e+110, 8.98249415358993e-134, 1.29676974102549e-244, 9.71395382456439e+148, 1.22431600491305e-220, 1.91520563704629e+32, 1.14421785779038e+60), .Dim = c(9L, 6L)), weights = c(1.32769666429078e-157, -3.11928183261913e+245, -1.86752067654946e+154, 3.24358230171322e-51, 1.87012097329086e+274, -4.73334455611919e+122, Inf, NA, -4.94790193947868e+281, 1.39866550655787e+253, 1.83055397731575e+88, 1.58249991451309e+305, -4.61089981522954e+157, -7.59555576528487e+161, 2.69991322252273e+30, Inf, 1.64130971007667e+230, -14533813584.7817, -Inf, -2.0546669045724e+260, 7.55438756680574e+53, 6.17657859271204e+95, -1.10506417849925e-124, -8.16688339212269e+98, -1.55169573295347e-154, 1.2577968277305e+50, 3.89993919211752e+169, 5.98883577194747e-30, 0)) result <- do.call(PLMIX:::Estep,testlist) str(result)