testlist <- list(gh = structure(c(1.66423450559976e-137, 1.47923599166225e-119, 1.37730436774025e+47, 8.81076968802032e+300, 5.84058871364184e+144, 9.28997407860819e-120, 2.0181295680432e+214, 6.61838458259232e+247, 4.273679258969e+85), .Dim = c(9L, 1L)), lp = c(-3.28862537274078e+138, 1.1875549927329e-142, 4.9206178857461e-122, Inf, 7.23481753349106e+108, -1.87070164228496e+161, -1.83597221414369e-141, 2.91136540596289e-276, 1.25617878902246e-306, 9.99945350734315e-115, 5.8316694984436e+134, 1.70869662838666e+189, -3.95014293295922e+142, 3.071625878611e+182, Inf, -3.61075659605699e+264, -1346.06711005202, -5.34005725049706e-16, 5.8316694984436e+134, -1.26720491978839e-83, 2.60220268485927e-188, 4.79334118015902e+211, 1.14696942417758e-158, 8.992295177772e+295, -9.23287137036662e-207, -Inf, 6.51205528365423e-202, 1.01965139772613e+248, -2.05403948307095e+276, -6.01076038841146e-214, -1.58547787017787e+293, 0), tau2 = c(NaN, 2.33992638335868e-17, 1.63839473377141e+277, -9.04414724421948e-120, NaN, 1.03002629081612e-25, NA, 5.05460263997121e-66, 1.13900018685178e-175, 3.13211940382004e-253, 8.28640699349771e+234, 8.18412733935947e-58, NaN, 4.64511378903821e+300, 1.20946187376266e-179, 5.21387897312266e-42, 1.44067678989322e-294, Inf), y = structure(c(1.19171081484103e-254, 7.00002790216982e-62, 2.23133658012291e-200, 4.5154730509857e+225, 5.74018005687698e+238, 4.08462408462852e+248), .Dim = c(1L, 6L ))) result <- do.call(robmixglm:::llrandbinomcpp,testlist) str(result)