- Operating System: Windows 10
- brms Version: 2.9.0
I’m trying to run a multi-level analysis with a 3 category outcome variable. However, even with the null model, the analysis runs slowly, about 30 minutes for one chain on my machine, and much longer when predictors, etc. are introduced.
I’ve attached the relevant data, and the analysis can be run with:
data <- read.csv("testdata.csv") m <- brm(data = data, family = categorical, dv ~ 1 + (1 | grp), chains = 1 )
This obviously uses the default priors, I have tried some different options (including gamma, normal, and different parameters to student_t), none of which have helped performance. The results (regardless of priors) look reasonable, and rhat/effective samples looks OK. Plot of results with 4 chains:
I am new to Bayesian analysis, so it’s highly likely I’m doing something wrong. In particular, I’m not confident I’ve used reasonable priors. I’d appreciate any help you can give me with this.
testdata.csv (158.0 KB)