我正在使用 R 包“survey”、“svy_VGAM”和“mitools”对复杂调查数据集的多重插补运行回归模型。我在 svyimputationList 上使用 svy_VGAM::svy_vglm()(在 'mitools' 中创建),然后将结果与 mitools::MIcombine 结合,得到一个 MIresult 对象。这些模型在模型族(累积、sratio、acat、bernoulli)和链接函数(probit、logit、cloglog)以及预测变量方面有所不同。
现在我需要比较模型。软件包的文档都没有讨论如何做到这一点。其他软件包中的诊断不接受“miresult”对象中的池化 svy_vglm 输出。例如,当我将未组合的 svy_vglm 模型从survey::with 转换为“mira”对象以进行池(使用mice::pool())或评估(使用mice::D1() 进行Wald 测试或mice:: D3() 用于似然比),我收到错误“svy_vglm 类的对象没有整洁的方法”
有哪些工具可用于在survey::with 或mitools::MIcombine 的输出上创建绘图和运行诊断?更广泛地说,任何人都可以建议如何比较多重估算的回归。重新加权使用该软件建模的复杂调查数据?
(请注意,我就同一研究的不同主题提出了不同的问题。对于出现冗余,我深表歉意。)