testlist <- list(est0 = structure(c(4.11945619609922e-261, 1.73747098328106e-41, 1.55819047957312e-257, 2.32409607243715e+23, 1.5141956870927e+296, 1.14815939103714e-21, 4.13264647308931e-208, 5.63993782138768e-162, 8.50245340521839e-28, 1.85794219917169e-180, 2.59150550635976e-144, 6.50424409531643e+115, 2.91345292979472e+92, 4.43386127810476e+213, 2.09698527056237e+215, 4.60530644481075e-152, 5.15216825771006e-69, 2.49573595602146e-258, 4.63190318661114e+126, 0.00139191715221442, 1.9378260399672e-70, 1.07751837710965e+131, 3.96947274531481e+136, 5.4911777326431e-220, 5.09434568597649e-246, 3.6545355249627e+28, 1087030835.19346, 1.70720776028565e-52, 5.19320996611745e-233, 1.69968642276295e-229, 1.84669130404175e+226, 5.02216810523442e+284, 9.20651795170416e+164, 1.36152278402067e+154, 1.25387370290657e-147, 2.14899299098352e+277, 4.64274589665337e+167, 9.1384709239327e-170, 1.16370208410278e-176, 4.45211802319474e+226, 2.18719618265509e+263, 5.82608891568927e-136), .Dim = 6:7), freq = c(-5.31109705587248e+295, NaN, -3.67619536836605e+248, -1.2223517612915e+220, Inf, -4.40989902893833e-97, 3.73613573697047e-92, -4.70930082326761e+109, -2.25792172691128e+216, -Inf, Inf, 3.39312226414051e+160, 0), genotype = structure(c(4.65752502202863e+228, 1.6761852271416e-138, 6.37285481713109e+175, 2.45057953873745e-99, 4.50263899575477e+252, 6.81556654178975e-47, 3.69071792574955e-201, 5.06045097099388e-109, 1.96378321214019e-44, 1.83369127091923e+156, 3.50749219530831e+135, 1.68045078281712e-259, 1826418857117.27, 1.34015805461918e-196, 2.3035219418568e-72, 1.70275753844619e-165, 2.79177158789627e-160, 1.17280775336274e-101, 8.96403949657468e-195, 5.24052089234381e-102, 7.25858050791017e+42, 9.59586576521786e-99, 1.23712227860646e-288, 1.1182579045061e-275, 4.04913811081941e-108, 3.42374227827782e+238, 8.48836867965598e+123, 1.67804000165913e-180, 2.67724275046118e-260, 8.12805865925918e+203), .Dim = 5:6)) result <- do.call(rres:::twostep_grm_cpp,testlist) str(result)