testlist <- list(prob = c(3.04019542923065e-47, 4.27116816578505e+255, 4.83272307438135e-168, -4.53603492193817e-296, -1.36240510310382e+101, 2.2448192363738e-237, -4.16939288998871e-302, -1.20774141549914e-151, 8.17076869196912e+223, 1.35966401640482e+191, 2.89570887577557e+77, 5.58090179482791e-287, 1.16790440721096e+248, -2.58641095696587e-225, -4.20067039531804e-224, 1.80509692244071e+216, 9.67629032843785e-100, -3.62904999789855e-95, 1.19981919581742e+67, 2.40623379107328e-261, 1.7546665350544e+103, 3.7812424901922e-109, 3.46098985818454e+22, 1.08858523712453e+178, -1.84613288735162e+209, -4.68616429940763e-61, -1.10511022542159e+51, 2.2257943099129e+291, 4.66000113105346e+190, -1.97655454417052e-29, NA, 3.67911305542372e+293, NaN, NA, NA, 5.39858484411856e+243, 2.65875615145775e-295, -6.20959448307596e-86, 8.35222282581511e-44, 0)) result <- do.call(SpatialEpi:::ProbSampleReplace,testlist) str(result)