MICE R的回归分析

Question

我最近更新了MICE程序包和一些我以前无法破解的代码... 我不太会使用MICE，因此非常感谢您提出修复建议

imputed_scabies_model <- with(coded_imputed_scabies,glm(impetigo_active~factor(scabies_infestation)+factor(gender)+factor(agegroups), family = binomial()))

  

将模型拟合到估算数据集

imputed_scabies_model_summary <- (summary(pool(imputed_scabies_model)))

  

获取摘要统计信息

imputed_scabies_model_OR <-exp(cbind(imputed_scabies_model_summary[,"est"],imputed_scabies_model_summary[,"lo 95"],imputed_scabies_model_summary[,"hi 95"]))

  

对估算值求幂，将95％CI的上下限设置为“ OR”和95％CI

imputed_scabies_model_summary <- (cbind(imputed_scabies_model_OR,imputed_scabies_model_summary))

  

将所有内容放在一起-所以我要有一个表，其中包含OR，95％CI和模型的原始输出。

我认为问题在于

的输出

imputed_scabies_model <- with(coded_imputed_scabies,glm(impetigo_active~factor(scabies_infestation)+factor(gender)+factor(agegroups), family = binomial()))
imputed_scabies_model_summary <- (summary(pool(imputed_scabies_model)))

这通常用于提供95％CI的估计值，上下限作为单独的列。 估计值仍然会出来（现在称为估计值，而不是“ est”，因此很容易修复这部分代码），但是现在，摘要似乎只是具有SE而不是CI范围。

![output of imputed_scabies_model_summary <- (summary(pool(imputed_scabies_model))) ] 1

所以我想我需要求幂然后手动生成CI-但想知道使用MICE的人是否能够提供一些指导。

我希望我的查询清楚。如果有帮助，我可以提供数据集。