我有相关数据,并且正在使用逻辑回归混合效应模型来估计感兴趣的预测变量的个体水平(条件)效应。我知道对于标准边际模型,使用 Wald 检验推断模型参数对于似然比和分数检验是一致的。它们通常大致相同。因为 Wald 很容易计算并且在 R 输出中可用,所以我 99% 的时间都在使用它。
但是,对于混合效应模型,我很感兴趣地看到固定效应的 Wald 检验与 R 的模型输出中报告的 Wald 检验与“手动”似然比检验之间存在巨大差异,该检验涉及实际上拟合简化模型。直觉上,我可以看到为什么这可能会产生巨大的差异,因为在简化模型中,随机效应的方差被重新估计并且可以显着影响可能性。
有人可以解释
- 如何在 R 中计算固定效应的 Wald 检验统计量?
- 混合效应模型中估计的模型参数的信息矩阵是什么?(并且计算 Wald 检验统计量的 mx 是否相同?)
- 在我描述的案例中,两个测试的结果在解释上有何不同?哪些通常是有动机的,并在文献中用于推理?