testlist <- list(k1 = 7.60578074503397e+105, k2 = -5.36624935815783e-226, k3 = 8.86680574344948e-84, kappa_opt = 0.00130059547935935, log_const_vmcos = -1.06121397865829e-21, log_I0_kappa_opt = -1.68797193534406e-189, logK = -5.90831401418114e+59, mode_1 = -1.98758422395412e+217, mode_2 = -3.81848449929176e-199, mu1 = 1.77803481973129e+225, mu2 = 3.67822021100512e+71, n = -388027017L, unifpropn = -2.97532338484649e+69, vmpropn = -7.72721898928751e-16) result <- do.call(BAMBI:::rcos_bimodal,testlist) str(result)