r - 在 lme4 中查看标准误差和参数估计

标签 r lme4 mixed-models

我刚开始使用 R 中用于线性混合模型的 lme4 包。我开始阅读 Douglas Bates 的书,以了解该过程。在书中,使用包中的“sleepstudy”数据集拟合模型的输出如下:

Linear mixed model fit by maximum likelihood [✬merMod✬]
Formula: Reaction ~ 1 + Days + (1 + Days | Subject)
Data: sleepstudy
AIC BIC logLik deviance
1763.9393 1783.0971 -875.9697 1751.9393
Random effects:
Groups Name Variance Std.Dev. Corr
Subject (Intercept) 565.52 23.781
Days 32.68 5.717 0.081
Residual 654.94 25.592
Number of obs: 180, groups: Subject, 18
Fixed effects:
Estimate Std. Error t value
(Intercept) 251.405 6.632 37.91
Days 10.467 1.502 6.97
Correlation of Fixed Effects:
(Intr)
Days -0.138

但是,当我拟契约(Contract)一个模型时,

fm06 <- lmer(Reaction~1 + Days + (1 + Days|Subject),sleepstudy,REML =FALSE)

我得到了不同的输出,没有标准误差估计。这是我的输出的样子:

    Linear mixed model fit by maximum likelihood  ['lmerMod']
Formula: Reaction ~ 1 + Days + (1 + Days | Subject)
   Data: sleepstudy
      AIC       BIC    logLik  deviance  df.resid 
1763.9393 1783.0971 -875.9697 1751.9393       174 
Random effects:
 Groups   Name        Std.Dev. Corr
 Subject  (Intercept) 23.781       
          Days         5.717   0.08
 Residual             25.592       
Number of obs: 180, groups:  Subject, 18
Fixed Effects:
(Intercept)         Days  
     251.41        10.47 

我不确定自己做错了什么(如果有的话),或者如何查看方差分量的参数和标准误差估计值。有人可以给我一些指示吗?非常感谢!!

最佳答案

尝试 摘要(fm06)。这给出了所需的输出。

关于r - 在 lme4 中查看标准误差和参数估计,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/30338971/

相关文章:

r - R中的置信区间计算

r - 如何使用 gtsummary::tbl_regression() 报告随机效应的方差估计

使用之前的估计值重新启动混合效应模型估计

r - 对于 lmer 模型,我可以提取整个模型的 'y' 的拟合值吗

r - Shiny 的应用程序 : How to dynamically change box title in server. R?

R 检查两个列表中的排序顺序

r - 如何调用使用 magrittr 管道创建的对象的元素?

r - 如何从 lme4 中仅提取公式的随机效应部分

r - 将时间序列合并到 R 中的混合效应模型中(使用 lme4)

r - 使用lme4设置rhoend参数