我正在使用 R 包lme4、merTools 和 Amelia来运行具有多重估算数据的多级模型。merTools 函数“print.summary.merModList”给出了模型结果,包括模型的固定效应部分(它给出了固定效应的估计值、标准误差和 t 值)。但是,我需要报告固定效应的 p 值。我了解多级模型的 p 值问题,并将在我的解释中关注置信区间,但我需要报告不符合惯例的 p 值。
我看到有几个选项可以在 lme4 包中显示 p 值(使用 Satterthwaite 近似值得出)。但是,当模型结果采用 merModList 对象(来自 merTools)的形式时,我似乎无法找到一种方法来处理多重插补数据。如果置信区间自动显示,而不是必须计算它们,这也会很有帮助。有谁知道这些事情是否可能?