Question

我正在使用lme4包进行线性混合效果建模

混合效应模型如下：

fm01 <- lmer(sublat <- goal + (1|userid))

上面的命令返回一个名为fm01

的S4对象

此对象包括固定效果及其OLS标准错误（如下）

Fixed effects:

            Estimate Std. Error t value
(Intercept)   31.644      3.320   9.530
goaltypeF1    -4.075      3.243  -1.257
goaltypeF2    -9.187      5.609  -1.638
goaltypeF3   -13.935      9.455  -1.474
goaltypeF4   -20.219      8.196  -2.467
goaltypeF5   -12.134      8.797  -1.379"

但是，我需要提供强大的标准错误

是否有人知道使用lme4返回的S4对象执行此操作的代码或程序包？

谢谢！

艾德丽安

Answer 1

我认为这就是您要寻找的：https://cran.r-project.org/web/packages/robustlmm/vignettes/rlmer.pdf

这是robustlmm程序包，它具有rlmer功能。

“在需要时，对象的结构和方法的实现应与具有坚固性特定扩展名的lme4尽可能相似。”

fm01_rob <- rlmer(sublat <- goal + (1|userid))

R的lme4封装中混合效应模型的稳健标准误差

1 个答案: