我想从一个lmermod对象中提取一些结果
require(lme4)
(fm1 <- lmer(Yield ~ 1|Batch, Dyestuff))
Run Code Online (Sandbox Code Playgroud)
这产生了
Linear mixed model fit by REML ['lmerMod']
Formula: Yield ~ 1 | Batch
Data: Dyestuff
REML criterion at convergence: 319.6543
Random effects:
Groups Name Std.Dev.
Batch (Intercept) 42.00
Residual 49.51
Number of obs: 30, groups: Batch, 6
Fixed Effects:
(Intercept)
1527
Run Code Online (Sandbox Code Playgroud)
特别是,我想提取2个随机效应的标准偏差,即42.00和49.51.我想可能有一个内置的方法来做到这一点,但我找不到它,所以我认为我应该能够通过使用找到它str(fm1),这产生了这个:
Formal class 'lmerMod' [package "lme4"] with 13 slots
..@ resp :Reference class 'lmerResp' [package "lme4"] with 9 fields
.. ..$ Ptr :<externalptr>
.. ..$ mu …Run Code Online (Sandbox Code Playgroud)