尝试使用HMI（hierch。多重插补）拟合两级插补模型......
我正在使用的模型是这个（我只想随机拦截）：
glmer(pica_yn ~ 1 + visit_c+visit_c2 + geo_child + hhloc + diar_c + hemo_c + (1|pid))
我一直收到这个错误：
buildZ中的错误（rmodel.terms [r]，data = data，nginverse = names（ginverse））：找不到对象id
似乎HMI更喜欢指定的公式也具有随机斜率。
是否有人适合 BINARY 响应的多级插补模型？
以下是您可以运行的示例，它将获得相同的错误：
data("sleepstudy", package="lme4")
sleepstudy[sample(1:nrow(sleepstudy), size = 20), "Reaction"] <- NA
sleep_formula<-Reaction ~ Days + (1|Subject)
hmi_imp <- hmi(data = sleepstudy, model_formula = sleep_formula, M = 5, maxit = 1)