[英]R-Squared of lmer model fit
I have a mixed effects model and I would like to see the R²- and p-value.我有一个混合效应模型,我想查看 R² 和 p 值。 I thought this is acessible by summary() but it's not.我认为这是可以通过 summary() 访问的,但事实并非如此。 Or at least I don't realize it.或者至少我没有意识到。
> 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
For the R², you can use r.squaredGLMM(fit1.lme)
from 'MuMIn
package. 对于R²,您可以使用'MuMIn
包中的r.squaredGLMM(fit1.lme)
。 It will returns the marginal and the conditional R². 它将返回边际和条件R²。
For the p-value, you can find them by using summary
with the lmerTest
package. 对于p值,您可以使用lmerTest
包中的summary
来查找它们。
For more information on p-values with mixed models: http://mindingthebrain.blogspot.ch/2014/02/three-ways-to-get-parameter-specific-p.html 有关混合模型的p值的更多信息,请访问: http : //mindingthebrain.blogspot.ch/2014/02/three-ways-to-get-parameter-specific-p.html
I add a very small demos with hierarchial modeling for ozone layer where the modeling acknowledges that it varies by month. 我添加了一个非常小的演示,其中包含臭氧层的分层建模,其中建模确认它随月份而变化。 You can find comparisons below. 你可以在下面找到比较。 I could find the R squared
term only in MuMIn
package. 我只能在MuMIn
包中找到R squared
项。
MuMIn package 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
where we compare the linear regression and the mixed effect model aka hierarchial regression model. 我们比较线性回归和混合效应模型aka层次回归模型。
Linear regression 线性回归
> 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 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 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
You can try packages sjPlot or sjstats.您可以尝试包 sjPlot 或 sjstats。 The first package help create APA-style tables from lme4 analyses, and the second package for extracting fit statistics.第一个包帮助从 lme4 分析创建 APA 样式的表,第二个包用于提取拟合统计数据。
You just need to simply write code:你只需要简单地编写代码:
tab_model(fit1.lme)
It will output an APA table including estimates slope, intercept, CI, p-value, variances, residuals, observation number, ICC, marginal and conditional R square etc.它将输出一个 APA 表,包括估计斜率、截距、CI、p 值、方差、残差、观察数、ICC、边际和条件 R 平方等。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.