lmer模型擬合的R平方

Question

我有一個混合效應模型，我想查看 R² 和 p 值。 我認為這是可以通過 summary() 訪問的，但事實並非如此。 或者至少我沒有意識到。

> summary(fit1.lme <- lmer(log(log(Amplification)) ~ poly(Voltage, 3) + (1 | Serial_number), data = bdf))
Linear mixed model fit by REML ['lmerMod']
Formula: log(log(Amplification)) ~ poly(Voltage, 3) + (1 | Serial_number)
   Data: bdf

REML criterion at convergence: -253237.6

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-14.8183  -0.4863  -0.0681   0.2941   9.3292 

Random effects:
 Groups        Name        Variance Std.Dev.
 Serial_number (Intercept) 0.008435 0.09184 
 Residual                  0.001985 0.04456 
Number of obs: 76914, groups:  Serial_number, 1270

Fixed effects:
                    Estimate Std. Error t value
(Intercept)         0.826745   0.002582     320
poly(Voltage, 3)1 286.978430   0.045248    6342
poly(Voltage, 3)2 -74.061993   0.045846   -1615
poly(Voltage, 3)3  39.605454   0.045505     870

Correlation of Fixed Effects:
            (Intr) p(V,3)1 p(V,3)2
ply(Vlt,3)1 0.001                 
ply(Vlt,3)2 0.002  0.021          
ply(Vlt,3)3 0.001  0.032   0.028  

Answer 1

對於R²，您可以使用'MuMIn包中的r.squaredGLMM(fit1.lme) 。 它將返回邊際和條件R²。

對於p值，您可以使用lmerTest包中的summary來查找它們。

有關混合模型的p值的更多信息，請訪問： http ： //mindingthebrain.blogspot.ch/2014/02/three-ways-to-get-parameter-specific-p.html

Answer 2

我添加了一個非常小的演示，其中包含臭氧層的分層建模，其中建模確認它隨月份而變化。 你可以在下面找到比較。 我只能在MuMIn包中找到R squared項。

MuMIn包

> data(airquality)

> MuMIn::r.squaredGLMM(lme4::lmer(data=airquality, Ozone ~ 1 + (1|Month)))
     R2m       R2c
[1,]   0 0.2390012
> summary(lm(data=airquality, Ozone ~ 1 + (1|Month)))$r.squared
[1] 0

我們比較線性回歸和混合效應模型aka層次回歸模型。

線性回歸

> summary(lm(data=airquality, Ozone ~ 1 + (1|Month)))

Call:
lm(formula = Ozone ~ 1 + (1 | Month), data = airquality)

Residuals:
   Min     1Q Median     3Q    Max 
-41.13 -24.13 -10.63  21.12 125.87 

Coefficients: (1 not defined because of singularities)
              Estimate Std. Error t value Pr(>|t|)    
(Intercept)     42.129      3.063   13.76   <2e-16 ***
1 | MonthTRUE       NA         NA      NA       NA    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 32.99 on 115 degrees of freedom
  (37 observations deleted due to missingness)

lmer4

> summary(lme4::lmer(data=airquality, Ozone ~ 1 + (1|Month)))
Linear mixed model fit by REML ['lmerMod']
Formula: Ozone ~ 1 + (1 | Month)
   Data: airquality

REML criterion at convergence: 1116.5

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.7084 -0.6269 -0.2669  0.4121  3.7507 

Random effects:
 Groups   Name        Variance Std.Dev.
 Month    (Intercept) 270.6    16.45   
 Residual             861.6    29.35   
Number of obs: 116, groups:  Month, 5

Fixed effects:
            Estimate Std. Error t value
(Intercept)   41.093      7.922   5.187

lmerTest

library(lmerTest)

> lmerTest::lmer(data=airquality, Ozone ~ 1 + (1|Month))
Linear mixed model fit by REML ['lmerModLmerTest']
Formula: Ozone ~ 1 + (1 | Month)
   Data: airquality
REML criterion at convergence: 1116.544
Random effects:
 Groups   Name        Std.Dev.
 Month    (Intercept) 16.45   
 Residual             29.35   
Number of obs: 116, groups:  Month, 5
Fixed Effects:
(Intercept)  
      41.09  
> summary(lmerTest::lmer(data=airquality, Ozone ~ 1 + (1|Month)))
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: Ozone ~ 1 + (1 | Month)
   Data: airquality

REML criterion at convergence: 1116.5

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.7084 -0.6269 -0.2669  0.4121  3.7507 

Random effects:
 Groups   Name        Variance Std.Dev.
 Month    (Intercept) 270.6    16.45   
 Residual             861.6    29.35   
Number of obs: 116, groups:  Month, 5

Fixed effects:
            Estimate Std. Error     df t value Pr(>|t|)   
(Intercept)   41.093      7.922  4.096   5.187  0.00616 **
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Answer 3

您可以嘗試包 sjPlot 或 sjstats。 第一個包幫助從 lme4 分析創建 APA 樣式的表，第二個包用於提取擬合統計數據。

你只需要簡單地編寫代碼：

tab_model(fit1.lme)

它將輸出一個 APA 表，包括估計斜率、截距、CI、p 值、方差、殘差、觀察數、ICC、邊際和條件 R 平方等。

看起來像這樣：