R中的Bucky包-多个协变量的全局测试？

问题描述

我在 R 中使用 Bucky 包为多重插补数据集生成稳健的标准错误。我想知道是否有某种方法可以对我的结果进行似然比检验或 Wald 检验（类似于 D1() function 中的 mice）。我需要比较两个模型（一个有交互项，一个没有交互项；交互项有多个类别，我需要对所有类别进行全局测试）。

我尝试将我的模型转换为 mira 对象，如下所示，但每次都得到相同的错误代码。有没有办法解决？谢谢！

ad.model.dep <- bucky::mi.eval(
  glm(
    w7.dep.bin ~ w6.shootmean + w5.dep.cont + w5.gad.cont + w5.pd.cont +
      sex + age + ethnic + lunch + w3.cyberbully + w5.connect + w5.db +
      w5.adhd + w5.alc + w5.nicmjdrugs,family = binomial(link = "logit"),data = HH.imputed
  ),cluster = HH.imputed$data$school
)

int.ethnic.model.dep <- bucky::mi.eval(
  glm(
    w7.dep.bin ~ w6.shootmean*ethnic + w5.dep.cont + w5.gad.cont + w5.pd.cont +
      age + sex + lunch + w3.cyberbully + w5.connect + w5.db +
      w5.adhd + w5.alc + w5.nicmjdrugs,cluster = HH.imputed$data$school
)

D1(as.mira(int.ethnic.model.dep),as.mira(ad.model.dep))

Error: No glance method for objects of class numeric

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

r r robust