2

有没有办法获得包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.fixefse.ranefse.coef

我的问题是:系数的 SE 是否只是 coef(ranef) + coef(fixef)?还是随机系数和固定系数的总和没有真正的标准误差?

4

0 回答 0