# NOT RUN {
library(mice)
imp <- mice(nhanes2)
# initially fit the model
fit_imp1 <- brm_multiple(bmi ~ age + hyp + chl, data = imp, chains = 1)
summary(fit_imp1)
# update the model using fewer predictors
fit_imp2 <- update(fit_imp1, formula. = . ~ hyp + chl, newdata = imp)
summary(fit_imp2)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab