Inconsistent results between regressionBF and lmBF
Hello! I am computing Bayes Factors in R and have noticed that when I use the regressionBF function in order to compare all additive models (i.e. y ~ .), the regressionBF function does not provide me with the same Bayes Factor for a given model as using the lmBF function. The regressionBF function consistently produces higher BF than the lmBF function. I am just running simple linear regressions, no interactions or categorical variables. Moreover, the top ("best") model using regressionBF actually has a lower Bayes Factor compared to the second-best or third-best models when using the lmBF function. Why is this inconsistent and which method is most reliable?
I'd greatly appreciate any input - thank you very much!
Comments
I'll attend Richard Morey to your question.
E.J.
Great, thank you very much!
Thank you very much for getting in touch - we figured out that the problem was due to the need to remove NAs, which led to inconsistencies when regressionBF was used for a subpool of the predictors. All fixed now, but greatly appreciate the responsiveness!