testlist <- list(actual = c(-1.59004957435928e-243, -2.66740479864599e-80, 3.33788380378054e-219, -1.33487386280268e-285, -4.07198978077703e-61, -1.42894563434117e+52, -1.71564460696604e+159, -1.28496133219007e-105, 4.31522135049808e-18, -1.45192531721821e+61, -5.43770437835069e+57, -1.44725859719644e-269, -0.733829340261768, -1.87225692964603e+209, NA, 5.84330094535224e-273, -2.05899621539983e-264, -7.04167440213677e-275, 3.50035709865439e+251, -1.52776061389597e-281, NA, 4.80808783195699e-133, 2.60726192716523e-38, -8.57205235344328e-279, 1.33514607168992e-17, -5.26535217823004e+209, -3.06548007257039e+151, NaN, 9.83649214540975e+199, -1.91460520194034e+216, 1.95348506435662e+125, -Inf, 3.12736019985354e+36, -6.60135728085622e-12, -2.72486047687397e+125, -3.39312424515957e+186, 2.29267547691833e+51, 1.47245123537047e+53, -1.19332538926796e+74, -Inf, 1.20000510219992e-81, 0), predicted = structure(c(1.5167956999958e-285, 3.08762902166214e+249, 1.1751831615252e-14, 2.22846958868524e-38, 9.37820611740662e+178, 5.80270469624338e-267, 2.68156787160659e-89, 6.17277696088128e-178, 7.76082303047441e+70), .Dim = c(3L, 3L ))) result <- do.call(ModelMetrics:::mlogLoss_,testlist) str(result)