glmer 摘要输出中缺少级别答案

【问题标题】：missing level in glmer summary outputglmer 摘要输出中缺少级别
【发布时间】：2015-10-26 09:25:20
【问题描述】：

关于 glmer 模型的摘要输出，我有一个令人难以置信的问题。

在下面的m0.1 中，我想知道Listgp 的第一级在哪里。它应该是 ListgpT。

m0.1 <- glmer(match ~ Listgp + (1|stimulus) + (1|Listener), data = PATdata, family = "binomial")
> summary(m0.1)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) [
glmerMod]
 Family: binomial  ( logit )
Formula: match ~ Listgp + (1 | stimulus) + (1 | Listener)
   Data: PATdata

     AIC      BIC   logLik deviance df.resid 
  5218.3   5253.4  -2604.2   5208.3     8203 

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-21.9276  -0.2804  -0.2059   0.2740   9.4275 

Random effects:
 Groups   Name        Variance Std.Dev.
 Listener (Intercept) 1.676    1.294   
 stimulus (Intercept) 4.949    2.225   
Number of obs: 8208, groups:  Listener, 228; stimulus, 12

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)  -1.3754     0.6792  -2.025   0.0429 *
ListgpTA      0.2284     0.3073   0.743   0.4572  
ListgpTQ      0.1432     0.2513   0.570   0.5687  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
         (Intr) LstgTA
ListgpTA -0.235       
ListgpTQ -0.288  0.636

从上面的输出可以看出，只显示了 ListgpTA 和 ListgpTQ，但没有显示 ListgpT。

这是否意味着ListgpT的结果与Intercept的结果相同？

【问题讨论】：

标签： r

【解决方案1】：

其实这很正常。在那种情况下，我猜测 Listgp 只有 3 种可能的情况：ListgpT、ListgpTA 和 ListgpTQ。所以 ListgpT 不需要输出，因为它被认为是参考电平，这意味着 ListgpTA 和 ListgpTQ 是与 ListgpT 比较来表示的。

您可以阅读this，它在R 中显示输出，您可以看到在 4 个 Cancer Stage 中只有 3 个在输出中（因为假设第一阶段是默认值，就像 ListgpT 一样您的示例中的默认值）

【讨论】：

非常感谢大卫。这真的很有帮助。