testlist <- list(a = -1.46036036480052e+48, lambda = -1.30600670447582e+305, mu = 1.71257617800484e-283, nu = -9.08327952567724e+305, obs_to_nodes_temp = structure(c(2.38345740568314e+154, 2.96436221675297e+236, 8.34546541608167e+217, -Inf, NA, 8.22382407453839e-197, -5.23158502952447e+73, 6.1281077234895e-174), .Dim = c(4L, 2L)), treetable_temp = structure(c(5.74433911083787e+70, Inf, 6.22697960364416e+119, 4.37746848029226e+95, 1.79161656269254e+266, NA, 1.01679418241044e+123, 6.61354695804704e-270, Inf, 9.07764854231226e-32 ), .Dim = c(5L, 2L)), y_temp = c(-4.40542167393296e-65, 9.93765015009146e+292, 2.32002580293826e+83, -3.79864241792839e-121, 1.8883108834625e-213, 2.45400007991628e-127, -5.16833556972955e+138, -6.51850214476773e+234, 5.70876483084426e+39, -1925.98874452245, 3.40619744136121e+177, 3.79124675916177e-141, -3.78688578449402e-26, 6.67997172957005e+201, 1.9529638921179e-225, 3.6712239261411e+174, -6.02246826322777e+297, NA, -3.64068544838652e+264, -1.02667353736304e+221, -4.93898503461023e-292, 1.18829283671081e+83, 4.54940249892308e+83, 6.14545433376592e-207, 2.68266245762969e+268, -1.49404189804505e+247, 7.68738535365778e-208, NaN, -1.21095565819795e-290, -9.56764834309804e-120, -1.23782442780946e+306, -6.58679767209119e-128, 6.87917021540263e+136, -4.09753193092625e-228, 5.93333407886065e-215, 2.89107199791954e+289, -1.646287653183e+219, 2.77718146285877e+36, 9.30500924864227e-126, -8.40171409844756e-48, -8.72773219496197e-222, 4.27064912662926e+243, 7.56179477482992e+267, 1.4852464645177e-139, 4.58275012087099e+262, NA, -8.2398052672162e-124, -1.61625624077693e-196, Inf, NaN, 3.46444855761134e-66, 0 )) result <- do.call(bartBMA:::likelihood_function2,testlist) str(result)