有没有办法获得包coef.merMod
方法的标准误差lme4
,或者计算系数的 SE 是无意义的?
例子:
library(lme4)
fit <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
coef(fit)
$Subject
(Intercept) Days
308 253.6637 19.6662579
309 211.0065 1.8475828
310 212.4449 5.0184061
330 275.0956 5.6529547
331 273.6653 7.3973914
332 260.4446 10.1951153
333 268.2455 10.2436615
334 244.1725 11.5418620
335 251.0714 -0.2848731
337 286.2955 19.0955699
349 226.1950 11.6407002
350 238.3351 17.0814910
351 255.9829 7.4520288
352 272.2687 14.0032993
369 254.6806 11.3395026
370 225.7922 15.2897506
371 252.2121 9.4791309
372 263.7196 11.7513157
这coef.merMod
只是每个分组因子的每个级别的每个解释变量的ranef
+总和。fixef
使用该软件包,您可以使用andarm
获得固定和随机效果的标准错误(和,但是,这只是两个第一个输出的串联)。se.fixef
se.ranef
se.coef
我的问题是:系数的 SE 是否只是 coef(ranef) + coef(fixef)?还是随机系数和固定系数的总和没有真正的标准误差?