我想使用 lmerTest R 包测试固定和随机效应。使用 lmer() 定义模型后,我使用summary() 来获取结果。对于我的结果报告,我想显示更多小数点的 AIC 和 BIC、随机效应中的方差和标准差、估计值和标准误差(参见图片: )。
我已经尝试过但没有帮助的事情:
options(digits=3)
summary(lmer, digits=3)
coef(summary(lmer), digits=3)
round(summary(lmer), digits=3)
round(coef(summary(lmer)), digits=3)
实现此目的的一种方法是加载
jtools
包并使用 options( "jtools-digits" = 3 )
。调用模型结果时,请使用 summ
函数查看它们。
目前无法调整 AIC 等打印的数字设置。您可以使用内部功能对此进行破解:
lme4:::.prt.aictab(summary(model)$AICtab, digits = 10)
AIC BIC logLik deviance df.resid
1763.9393445 1783.0970856 -875.9696722 1751.9393445 174
全局
digits
选项应该会影响 coef(summary(.))
的打印:
options(digits = 10)
coef(summary(model))
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 251.40510485 6.632122742 18.00113477 37.907185169 1.263751938e-18
以及随机效应方差:
VarCorr(model)
Groups Name Std.Dev. Corr
Subject (Intercept) 23.7797596
Days 5.7167985 0.081321
Residual 25.5919070
要提取混合模型的组件,另请参阅 this Question 以及
glance()
包中的 tidy()
和 broom.mixed
函数。
emm_options(opt.digits = FALSE)