r - 比较 svy_vglm 回归模型的结果（'survey','svy_VGAM','mitools','mice' in R）

Question

我正在使用 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 的输出上创建绘图和运行诊断？更广泛地说，任何人都可以建议如何比较多重估算的回归​​。重新加权使用该软件建模的复杂调查数据？

（请注意，我就同一研究的不同主题提出了不同的问题。对于出现冗余，我深表歉意。）

标签： rlinear-regressionsurveyr-micevgam