testlist <- list(alpha = c(-1.82345480893601e-134, 7.02932507688266e-145, -8.78031930250735e-227, 1.80853666575964e-35, -Inf, -5.15129689039746e+128, -2.51243585603214e+75, -2.20659967285933e-218, 9.72521801287326e-151, 3.09284972667018e+282, 3.84688775539034e+199, 4.96544452797543e+275, NA, -3.86949137960274e-63, -3.9686709273731e+162, -2.36742817649074e-145, 3.49051345224257e+284, -Inf, -1.08668495760388e-63, -3.58792751784694e-298, Inf, NA, 0), data = structure(c(2.88847595190092e+209, 14089863137281778 ), .Dim = 2:1), lambda = c(5.13379589103216e+147, -6.00714912133061e-230, -3.87513376508027e+65, Inf, -1.20900589856212e+21, -1.05137092786561e+200, Inf, -2.23198166279758e-273, -1.59988685996208e-161, -2.62926872681595e+281, -3.9384975798735e-272, 3.17437447992959e-302, 8.26589957370324e-227, NaN, -1.96111948367253e+79, 3.42219086703652e+291, 1.73563800728723e-205, 3.94789160457658e+20, -1.45900258945393e-167, -6.96701993564468e+147, -8.2090881081728e-249, -Inf, 0)) result <- do.call(mixR:::EXgamma_C,testlist) str(result)