问题描述
我正在处理健康数据，其中缺少大约 360 个身体质量指数 (BMI) 变量值。我已经使用 mouse 包估算了这些值，现在我希望运行逻辑回归模型来模拟一个人患高血压（虚拟变量）的可能性，其中 BMI 是一个自变量。在插补之前，我已经使用 logitmfx 函数编写了我的代码，因为它似乎是评估逻辑模型的最直接方法，其系数以均值 (MEM) 评估并且具有稳健的聚类标准误差。这是我的代码对于不使用插补数据的模型的样子。
logit_margin3<-logitmfx(HYPERTEN~AGE+scale_BMICALC+FAMSIZE+ BELOW_TWICE_POVERTY
+IMMIGRANT+FEMALE+REGION_1+REGION_3+REGION_4+YEAR_2015+YEAR_2016+YEAR_2017+YEAR_2018+HLTH_INS+SMOKE+MARRIED+NON_WHITE+DIABETES,data=df,atmean=TRUE,robust=T,clustervar1="PSU",clustervar2=NULL)
使用估算数据，我希望能够运行这样的程序。 imp.data.master 是插补数据集的名称
hypertension <- with(logitmfx(HYPERTEN~AGE+scale_BMICALC+FAMSIZE+ BELOW_TWICE_POVERTY
+IMMIGRANT+FEMALE+REGION_1+REGION_3+REGION_4+YEAR_2015+YEAR_2016+YEAR_2017+YEAR_2018+HLTH_INS+SMOKE+MARRIED+NON_WHITE+DIABETES,clustervar2=NULL)))
summary(pool(hypertension))
但是，根据我在网上阅读的内容，一旦我估算了 BMI 值，似乎无法使用 logitmfx 函数输入估算数据，因为我收到错误：“logitmfxest(formula,数据，平均，稳健，clustervar1，clustervar2，： 数据参数必须包含 data.frame 对象"
我在网上读到可以将推算数据输入到我可以做的 glm 函数中。然而，margins 包似乎并没有让使用强大的 SE 和在 MEM 上进行评估变得那么简单，所以我想先在这里检查一下，以防有人知道使用 logitmfx 的方法。谢谢！
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）