在 lme4 中查看标准误差和参数估计
Viewing standard errors and parameter estimates in lme4
我刚开始使用 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
但是,当我拟合同一个模型时,
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
我不确定我做错了什么,如果有的话,或者如何查看方差分量的参数和标准误差估计。有人可以给我一些指示吗?非常感谢!!
试试看
summary(fm06)
。这给出了所需的输出。
我刚开始使用 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
但是,当我拟合同一个模型时,
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
我不确定我做错了什么,如果有的话,或者如何查看方差分量的参数和标准误差估计。有人可以给我一些指示吗?非常感谢!!
试试看
summary(fm06)
。这给出了所需的输出。