powered by
test.mh(dl, psi = mean(dl$gene.map$psi.hat), omega = mean(dl$gene.map$omega.hat), num.cores = getOption("DL.num.cores", max(parallel::detectCores() - 1, 1)), iterations = 1000, thin = 15)