混合效应逻辑回归的固定效应解释


10

我对UCLA网页上有关混合效应逻辑回归的陈述感到困惑。他们显示了一个通过拟合这样的模型得到的固定效应系数表,下面的第一段似乎完全像正常的逻辑回归那样解释系数。但是,当他们谈论比值比时,他们说您必须根据随机效应来解释它们。是什么使对数奇数的解释与其指数值不同?

  1. 要么不需要“保持其他所有内容不变”?
  2. 从该模型解释固定效应系数的正确方法是什么?我总是给人的印象是“正常”逻辑回归没有任何变化,因为随机效应的期望值为零。因此,无论有无随机效应,您都可以解释完全相同的对数奇数和优势比-仅更改了SE。

估计可以基本上一如既往地解释。例如,对于IL6,IL6的单位增加与预期的缓解对数几率中的0.053单位降低相关。同样,已婚或已婚的人有望获得的缓解几率比未婚者高0.26。

许多人喜欢解释赔率。但是,当存在混合效果时,这些将具有更细微的含义。在常规logistic回归中,优势比与所有其他预测变量固定的预期优势比成正比。这是有道理的,因为我们经常对统计数据进行调整以适应其他影响(例如年龄),以获得结婚的“纯”影响或任何主要的主要预测指标。混合效应逻辑模型也是如此,此外,保持所有其他特征不变包括保持随机效应不变。也就是说,赔率是指保持年龄和IL6恒定的人以及具有相同医生或具有相同随机效应的医生的人的条件赔率


我可能是错的,但对此表示怀疑。对于赔率比而不是对数赔率差异,没有特殊考虑。保持其他一切不变意味着要以剩余的固定效应和随机效应为条件。“预期已婚或已婚人士的缓解率比单身人士的对数几率要高.26”,应增加“如果他们具有相同的年龄,ILS和随机截距值”。这是一个简单的旧方程式。
异方体Jim 18'Sep 8'3

Answers:


19

确实,在混合效应逻辑回归中,并且由于使用非线性链接函数将结果的平均值与线性预测变量联系起来,固定效应系数的解释取决于随机效应。

一个简单的例子如下:假设您有一个多中心临床试验,其中每家医院的患者被随机分配到两种治疗方法,即A或B。还应说明所关注的结果是二元的(例如,患者需要手术,是或否)。为了说明该试验的多中心性质,我们将混合效应逻辑回归与每家医院的随机效应进行拟合(即随机截距模型)。从该模型中,我们获得了治疗变量的回归系数,即。这个是来自同一患者的两种疗法之间的对数比值比ββ医院。现在,如果您使用广义估计方程(GEE)方法分析了相同的数据,那么您将获得带有边际解释的系数。继续上面的示例,从GEE 得出的估计系数将是医院患者的两种治疗之间的对数比值比-换句话说,对整个医院平均的对数比值比。β

有多种方法可以从混合效应逻辑回归中获得具有边际解释的系数。有关更多详细信息,请参阅我的课程笔记的5.2节。对于这种方法,以获得与来自GLMM边际解释,检查函数的系数为r的执行marginal_coefs()GLMMadaptive包; 更多信息也可以在这里找到


感谢您的明确答复!您的笔记看起来很棒,我希望讲座能在线进行!
B_Miner

您能否确认,这些解释是否也适用于线性混合模型(不仅仅是glmms)
B_Miner

1
在线性混合模型中,系数同时具有边际解释和特定主题解释。
Dimitris Rizopoulos '18

谢谢。这是否意味着只要不对系数进行变换(例如,取幂),就可以一目了然并且特定于主题?因此,对于逻辑混合模型,只要对系数的解释是对数奇数,就可以同时对它们进行双向解释?
B_Miner

2
不,即使您不求幂,对数赔率仍具有特定于主题的解释。也就是说,在混合效应逻辑回归中,您可以对建模。如果对随机效果的分布进行期望,则会得到,即固定效果部分。但是,这是边际对数赔率。logPr(Y=1|b)1Pr(Y=1|b)XβEb{logPr(Y=1|b)1Pr(Y=1|b)}=XβlogEb{Pr(Y=1|b)}1Eb{Pr(Y=1|b)}
Dimitris Rizopoulos '18
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.