我有一个具有固定和随机效果的mer对象.如何提取随机效应的方差估计？这是我的问题的简化版本.
study <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)
study
这样可以提供长输出 - 在这种情况下不会太长.无论如何,我如何明确选择
Random effects:
Groups Name Variance Std.Dev.
Subject (Intercept) 1378.18 37.124
Residual 960.46 30.991
部分输出？我想要自己的价值观.
我长期看看
str(study)
那里什么都没有!还检查了lme4包中的任何提取器功能都无济于事.请帮忙!
我在R中使用nlme包的lme函数拟合线性混合模型.我的目标是找到固定效果部分的线性系数组合的p值.
Value Std.Error DF t-value p-value
(Intercept) -1.5897634 0.019578359 41547 -81.20004 0.0000
group_k 0.0866493 0.027662455 594 3.13238 0.0018
x1 0.1621744 0.006963385 41547 23.28960 0.0000
例如,在此结果中提供假设$\beta_1 = 0 $的p值(0.0018).但是如果我想找到一个假设的p值,例如$\beta_1 + 3*\beta_2 = $,我该怎么办？是否有任何现成的命令来进行此计算？