new version of multisample, with added bayes factor calculation