我有一个混合效应模型,我想看看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
3条答案
按热度按时间zlwx9yxi1#
对于R²,可以使用
‘MuMIn
包中的r.squaredGLMM(fit1.lme)
,它将返回边际和条件R²。对于p值,您可以使用
summary
和lmerTest
包来查找它们。有关混合模型p值的详细信息:http://mindingthebrain.blogspot.ch/2014/02/three-ways-to-get-parameter-specific-p.html
xkrw2x1b2#
我添加了一个非常小的臭氧层层次模型演示,其中模型承认它随月份变化。您可以在下面找到比较。我只能在
MuMIn
包中找到R squared
术语。MuMin封装
其中我们比较了线性回归和混合效应模型,即分层回归模型。
聚合物测试
nnt7mjpx3#
您可以尝试sjPlot或sjstats包,第一个包帮助从lme4分析中创建APA样式的表,第二个包用于提取拟合统计数据。
您只需要简单地编写代码:
它将输出APA表,包括估计值斜率、截距、CI、p值、方差、残差、观察数、ICC、边际和条件R平方等。
看起来像这样: