testlist <- list(alpha = c(8.21083511161007e+163, 1.10688144538961e+126, -8.69789758027628e-52, -8.3972444924628e+186, 3.14646635030145e+234, 5.31603299873875e-223, -337233663.773242, NA, -3.9029839244136e+155, -1.02071570419833e-272, 7.10798261430425e-46, 1.27274258542366e-144, 4.3092349856493e+97, -4.30095413405714e-146, 1.88583789258482e-188, 1.33232887448596e-288, 2.47910575804432e-94, -2.30338706700092e+88, -1.66938394026492e-96, 2.16682471779353e+228, 6.29121214822243e+40, 7.22274072916462e+190, 9.17429601042768e-81, -1.57617635527491e+124, -54146171262928461824, 1.57233921267359e-199, -2.58010589285433e+29, -1.09666583018852e+192, -1.30334599788956e-114, 1.10876539286829e-217, -2.61919194548597e-06, -3.03810190673057e-83, -3.90905637900774e+21, 5.77524080248576e+280, -1.39647699963017e-142, 7.57965997496759e-206, 1.98979990574981e+90, -1.07443545334824e+29, -0.0253855879133046, 80919.4917994781, NA, -2.32687456681762e-84, 0.976449103364431, 6.39139606658221e+183, 2.85869968805415e-36, 63624651704919719936, 1.39394044627572e+82, 5.34752537209056e-269, -1.81822992922295e+192, 9.92942280738504e-149, -1.28351599539849e+40, 4.42589329271277e-91, -1.13337472106739e+73, -6.46758239094001e-77, -3.31444017638099e+79, -1.33339927039387e-28, 8.33605794215744e-137, 4.91005599842903e-305, -4.6484290796306e-208, 9.70093351989796e-45, 1.89120685072928e-157, -1.55602457200731e+292, -1.13964579110389e-283, 6.43231133718401e-227, 2.56889960753463e-270, 3.58471969852271e+226, -1.86167417672263e+21, 8.84550205243477e+199, -2.70251817798242e-56, -1.26564888726293e-63, 2.3249440746566e-181, -5.09453558014039e+296, 4.40892695063235e+105, -2.17912322306063e+196, -Inf, -2.11181908662958e+209, 3.55212703739538e+27, -2.92516674299643e-103, 3.03626246329939e-86, 9.83557849153725e-195, NaN, -5.99173908602851e+203, 0), data = structure(c(1.18002797449058e-76, 4.57528287727422e+98, 2.08682241377383e+180, 4.40227796480162e+304, 4.61999803557341e+128, 1.04004027501579e+144, 1.23787989969351e+185, 4.64613976164796e-11, 3.81770479378804e-209, 9.26273057394056e+249, 1.13415503480368e-154, 2.24347383891449e+48, 6.26811209759216e+229, 3.79972223736023e+72, 1.3780337647634e+172, 1.18902870403196e-217 ), .Dim = c(8L, 2L)), lambda = c(-3.73799158832719e-71, -4.72107567274611e+116, -1.81905895936073e+100, -3.33382281692081e+22, 7.70336837497828e+165, 5.25559666932077e+246, -Inf, 8.96165642986678e-121, -2336320087160872448, -1.74743010589876e+135, -1.15551450922664e-150, -8.81869635415787e+196, -5.79393572449354e+39, 5.32307690738095e+256, -9.09826577161754e+285, NA, 9.25455114391888e+191, -1.02680819525692e-47, -3.95522577587436e+125, 161535349953.319, 7.91841371482742e-36, 3.54238446083392e-76, 7.5529571146842e+103, NA, 3.00226602165172e+172, -2.66940339669364e+201, -2.98526446972642e+284, 2.20942694505107e+52, -2.02098433426802e+64, -4.18169364391881e+279, NA, 1.97935992706659e+149, -1.57932320891757e-220, -8.18694155632965e-171, 0), pi = c(35153972759995532, 1.33005626801739e+133, 6.63381342296374e+291, 1.86308008142289e-301, -6.95145940764727e-185, 4.35239755333997e+150, 1.05310441261224e+87, -2.51982634500719e+91, Inf, NaN, 4.69606756114336e+61, -1.36083780683473e-120, -1.95259536606243e-149, -2.18110565083858e+295, NaN)) result <- do.call(mixR:::loglik_gamma_g_C,testlist) str(result)