Hello everyone
I hope you are doing well
I have run a bayesian multilevel model. The equal variances model, and specific effects model worked well. However, the unequal variances model seem to take a very long time. I waited for two days (few times), but it did not make any progress. Can I rely only on equal variances and specific effects models for my study without running the unequal variances model?
here is my brms model:
library(brms)
fit_participation1md ← brm(
bf(
mvbind(sign_pet, join_boycott, attend_demo, join_strike) ~
happiness * var_scaled +
health * var_scaled +
free_choice * var_scaled +
life_satis * var_scaled +
financial_satis * var_scaled +
Employed + Retired + Student + Unemployed + Male + careful_trust + (1 |p| country),
disc ~ 1 + (1 |p| country)
),
data = final_dataset,
family = cumulative(link = “logit”),
control = list(adapt_delta = 0.99, max_treedepth = 15),
chains = 4, cores = 4, iter = 2000, init_r = 0.1
)
Respectfully,
MD