R Code for Inclusion Bayes Factors across Matched Models

ldh_stats

Hello,

I have been using JASP to run Bayesian ANCOVAs to obtain inclusion Bayes Factors (BFincl) across matched models (following: https://www.cogsci.nl/blog/interpreting-bayesian-repeated-measures-in-jasp). Given that JASP uses R to run analyses, is there available code to run these exact analyses within R? 

Thanks in advance!

EJ

I've passed on the issue to the JASP team member with the relevant expertise.

EJ

ldh_stats

Thank you!

ldh_stats

@EJ I was wondering if there is any update on the R code for running Bayesian ANCOVAs to obtain inclusion Bayes Factors?

kuba22

https://projektowanie.art/

Yep, you can run similar Bayesian ANCOVAs in R using the `brms` package (which works with Stan) for Bayesian models. You’ll just need to set up your models and run them with the `bf()` function, and then extract Bayes Factors using `bayes_factor()`. It’s pretty straightforward once you get the hang of it.

EJ

Yeah but the complication with brms is that the prior distributions need to be set up with care.

I have reminded our expert :-)