1

为什么不是confint.default基于渐近正态性的模型不适用于lmer模型？

fit <- lmer(y~(1|operator)+(1|part),data=dat)

Linear mixed model fit by REML ['lmerMod']
Formula: y ~ (1 | operator) + (1 | part)
 Data: dat
REML criterion at convergence: 409.3913
Random effects:
Groups   Name        Std.Dev.
part     (Intercept) 3.2018  
operator (Intercept) 0.1031  
Residual             0.9398  
Number of obs: 120, groups:  part, 20; operator, 3
Fixed Effects:
(Intercept)  
     22.39  


 confint.default(fit)
 Error in as.integer(x) : 
   cannot coerce type 'S4' to vector of type 'integer'

错误在说什么？如何根据 lmer 模型的渐近正态性获得置信区间？

4

2 回答 2

4

不使用confint.default()，只使用confint()。对于不同的模型类型，计算置信区间的方法是不同的。你可以看到不同的方法methods(confint)。函数的“正确”版本是根据传递给函数的第一个对象的类调用的。直接调用其中一种方法通常不是一个好主意。

有用于计算lmer对象边界的选项。查看帮助页面?confint.merMod以查看该模型类型独有的选项。

于 2015-08-01T01:12:56.300 回答
2

@MrFlick 是正确的，但可能值得添加的是confint.merMod()默认情况下给出可能性配置文件 CI；confint(.,method="Wald")将给出基于渐近正态性的置信区间：

'“Wald”'：基于似然面的估计局部曲率近似置信区间（仅限固定效应参数；所有方差-协方差参数 CI 将返回为“NA”）；

（这在帮助页面上很明显，但可能值得在这里重申）。

于 2015-08-01T01:44:28.257 回答