Hi

I am working in some missing data examples with brms using the model imputation with mi()

The model works fine, but I am wondering in the multivariate formula method, the formulas are technically independent. So, with multiple variables that have missing data as predictors, I would like to include as many of the model variables in missing data argument to help it.

In imputation methods this usually means that the same variables are both outcomes and predictors of missing data across variables, in methods like the one use by mice.

So, if I do the same in brms, would this be an appropriate model as the equations are independent, or am I creating some misspecified feedback loop?

With something like the following example, where for example mi(rumination) predicts depression | mi(), and also mi(depression) predicts rumination | mi()

```
bf_dep <- bf(depression | mi() ~ mi(rumination) + mi(effortControl) + mi(anxiety) + sex + age)
bf_rum <- bf(rumination | mi() ~ mi(depression) + mi(effortControl) + mi(anxiety) +sex + age )
fit_imp2 <- brm(mvbf(bf_dep, bf_rum, rescor = F),
data = dat_mis,
chains = 3,
iter = 3000,
family = "gaussian")
```

Thank you