testlist <- list(a1 = -2.06316751081532e-192, b = 1.50758697244514e-201,      n = 113097773L, resids = c(8.53510522205566e+40, -2.11608220844914e-286,      1.49502318067595e-80, 2.76932739972246e-111, 9.85767125196701e+119,      -1.22173618763461e+54, -2.25867464534669e+73, -99279963.9237618,      3.16395718553144e+42, 7.26928186503983e+50, -1.7545381372656e+183,      4.13252411068885e-50, 3.17562122814866e-266))
result <- do.call(bartBMA:::update_sigma,testlist)
str(result)