我有一个固定和随机效果的mer对象.如何提取随机效应的方差估计?以下是我问题的简化版本.

study <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)
study

这会产生很长的输出,在这种情况下不会太长.无论如何,我如何明确地 Select

Random effects:
Groups   Name        Variance Std.Dev.
Subject  (Intercept) 1378.18  37.124  
Residual              960.46  30.991  

部分输出?我想要价值观本身.

我看了很久

str(study)

那里什么都没有!还判断了lme4包装中的任何提取器功能,但均无效.请帮忙!

推荐答案

lmer返回一个S4对象,所以应该可以:

remat <- summary(study)@REmat
print(remat, quote=FALSE)

其中打印:

 Groups   Name        Variance Std.Dev.
 Subject  (Intercept) 1378.18  37.124  
 Residual              960.46  30.991  

...一般来说,您可以查看"mer"对象的printsummary种方法的来源:

class(study) # mer
selectMethod("print", "mer")
selectMethod("summary", "mer")

R相关问答推荐

插入指示行之间时间间隔的新行

为什么stat_bin在R中的ggplot中显示错误的数字?

为什么以及如何修复Mapview不显示所有点并且st_buffer合并一些区域R?

在数据表中呈现数学符号

selectInput不返回ALL,并将因子转换为shiny 的数字

用相同方法得到不同函数的ROC最优截断值

RStudio中相关数据的分组箱形图

gganimate在使用shadow_mark选项时不保留所有过go 的标记

如何使下一个按钮只出现在Rshiny 的一段时间后?""

R s iml包如何处理语法上无效的因子级别?'

条形图顶部与其错误条形图不对齐

如何移除GGPlot中超出与面相交的任何格网像元

如何在R中改变fviz_pca_biplot中圆的边界线的 colored颜色 ?

错误包arrowR:READ_PARQUET/OPEN_DATASET&QOT;无法反序列化SARIFT:TProtocolException:超出大小限制&Quot;

向R中的数据帧添加一列,该列统计另一列中每个唯一值的二进制观测值的数量

有没有办法通过str_Detect()或其他字符串匹配函数来连接两个长度不等的数据帧?

如何显示准确的p值而不是<;0.001*?

在r中整理图例和堆叠图的问题

将美学添加到ggploy中的文本标签

如何编辑被动式数据表?