testlist <- list(gh = structure(c(1022.16495558765, 8.17651269005385e+251, 1.18557932922793e-164, 1.32136921557247e-306), .Dim = c(1L, 4L )), lp = c(-2.12316316700847e+76, -2.12316316700847e+76, 7.02983088517237e+124, NaN), phi = c(NA, Inf), tau2 = c(-1.18057813378067e+37, -1.83330193130014e+113, 3.60379930989372e-21, -1.14976641599168e-180, -1.63432115105052e-34, -3.09901150179119e+272, 1.4764013208548e-34, 1.1612035211291e-167, -2.92227247008808e+39, 3.30645975531507e+47, 3.26444443034013e+125, 3.84607436453304e-286, 1.85905499831123e-208, Inf, 2.1086116312455e-73, NA, -8.43264431930967e+291, 1.37704952437755e+292, 1.07750378732186e+248, -1.00800499448322e+275, -9.03285374431099e-181, -4.42531426422316, -1.65996166193112e-251, 1.78622116561591e-105, 5.55125378945135e-253, -1.29020037367523e-284, NaN, -1.25445385845547e-242, 4.85235959044992e+263, -5.11714529446299e+22, -3.22549497267183e+188, -0.00736304351370888, -2.59439360412874e+95, -1349913.44122967, 8.07718299385497e-114, -1.63299872082322e-122, -7.12837701153695e-239, 5.3630781763729e-39, -1.378326026943e+169, 3.81208662135388e+100, -2.21288319987785e-163, 3.34227427104536e-139, 1.45813597778465e+35, NaN, -6.73641384092395e-48, 1.3603483735641e+286, 0), y = c(-1.399203775067e+41, NaN, -3.96152151305233e-91, -2.52067914084328e-138, -6.16526036387348e-36, 1.03634467518301e-28, 1.53896677218411e+185, 17102780259293784, 1.15417090967161e+232, -9.02213642044339e-187, -1.29423606977447e-265, Inf, -6.58303688038407e+114, 5.92840128038097e-295, -1.72121662039347e+262, -5.05326121402449e+75, -Inf, 1.19154342744768e+219, -6.17362032525506e+192, 5.03772717045048e-178, 5.24891359239829e-133, -108259.147925947, 1.04204092621642e-215, -1.61343500540468e-240, -1.02612402337211e-187, -1.57368328428706e-233, 4.28951558111186e+87, -1.91546224840815e+92, Inf, -1.61614597728681e-83, 3.79371330527601e-86, -1.048371684886e+40, -5.95089393925176e-305, -6.14598392817378e-27, 1.44009611772388e+276, -5.9688844444255e+68, -1.83154432291392e+292, -8.92988655890896e+142, 1.03392593660272e+150, 0)) result <- do.call(robmixglm:::llrandgammacpp,testlist) str(result)