add function to compute bayes factor estimate

over replicates in parallel