LRT 比较随机效应模型和嵌套逻辑回归模型

机器算法验证 r 物流 混合模式 随机效应模型 似然比
2022-04-18 18:37:55

我有结构的逻辑回归模型,以及结构y ~ x1 + x2的具有随机截距和随机斜率的广义线性混合模型 (GLMM) y ~ x1 + x2 + (1 + x2 | id)我的目标是确定是否需要混合模型。这两个模型是否被认为是“嵌套的”,它们可以与似然比检验 (LRT) 进行比较吗?我将如何确定该测试的自由度?

如果可以将它们与 LRT 进行比较,那么如何在数学上陈述零假设?,其中是 GLMM 模型的随机效应向量。适当的零假设是吗?G=cov(bi)bi2×1H0:G=0

如果假设的上述设置是正确的,那么如何识别结果检验统计量的分布。我知道测试的 50:50 混合,因为零假设在边界上的参数空间。但就我而言,我想测试整个矩阵(即)是否为零。H0:g22=0χ12χ22G2×2

1个回答
  • ,则混合模型会简化为更简单的模型(这与相同,因为如果方差为零,则协方差必须为零,但是根据上的联合条件来说明它可能更容易理解。)σ12=σx22=0G=0{σ12,σx22}
  • 通常形式的似然比检验不能正常工作——它是保守的——因为似然比检验的推导取决于空参数周围对数似然的泰勒展开,如果空参数是在可行模型空间的边界上(不能围绕σ2=0,因为这意味着您在扩展中包含负方差值)。这在很多地方都有讨论(Self 和 Liang 1987;Stram 和 Lee 1994;Goldman 和 Whelan 2000;Pinheiro 和 Bates 2000)。对于简单模型,通常的零分布有一个已知的校正因子。例如,如果您正在测试由单个方差参数不同的模型(例如随机截距模型与非随机截距模型),则的空分布为,其中2Δ(logL)0.5χ02+0.5χ12χ02是零点质量;这里的底线是名义 LRT p 值应该除以 2。对于更复杂的模型,通常很难推导出来,人们通常通过参数引导来计算 p 值。GLMM FAQ 有一个关于这个的部分......

特别是,Stram 和 Lee (1994) 讨论了一些更复杂情况的几何学(我已经很长时间没有阅读它了......)形成零分布 s 的特定混合可能是解析可推导,但根据我的经验,人们通常会放弃并通过模拟找到零分布。下面的例子来自 Pinheiro 和 Bates (2000) p。87(通过谷歌图书):他们通过计算表明,特定比较(即vs. )的空分布约为然后他们或多或少地说他们继续使用天真的轻轨,因为它更容易。χ21|Worker1|Worker/Machine0.65χ02+0.35χ12

如上面链接的 GLMM FAQ 部分所示,您可以使用pbkrtest::PBmodcomp()通过参数引导来获取有效的 p 值......

在此处输入图像描述


Stram、Daniel O 和 Jae Won Lee。“纵向固定效应模型中的方差分量测试。” 生物识别 50,没有。4 (1994): 1171-77。