I have been using brms to make a mixed-effects model, that tries to calculate the hospital random effect intercepts, to determine the effect of each hospital on the mortality of patients.
The variable Risk Adjustment just represents the predicted percentage mortality of the patients obtained from a previous model.
brms_glmer <- brm(Mortality ~ offset(RiskAdjustment) + (1|Hospital), data = HLM_data,
family = bernoulli (link = "logit"),
cores = 16,
chains = 4,
iter = 2000,
seed = 1234
)
My question: is there anyway to specify the variance of the means of the hospital intercepts? Or some way to increase the variance of the hospital mean intercepts generated by the model (to get a wider distribution).
Thanks for your help in advance!