有没有办法为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的总和。
使用arm包，您可以使用se.fixef和se.ranef（以及se.coef获取固定和随机效果的标准错误，但这只是一个连接在两个第一产出中。）
我的问题是：系数的SE是简单的coef（ranef）+ coef（fixef）吗？或者随机和固定系数的总和没有真正的标准误差？