samp_bf function

Sample multiple datasets and compute the Bayes factor in each