简体   繁体   中英

Warnings calculating confidence intervals with mixed models

I'm fairly new to R, but have been running models using lmer. Now I'm trying to calculate the confidence intervals for the intercept and the condition.

For example:

Model1 <- lmer(Response ~ Treatment + 
               (1+Treatment|Hospital), data=Data) 
confint(Model1)

I receive the following output and warnings such as:

              2.5 %     97.5 %
.sig01         5.6742125  6.7891223
.sig02        -0.9234187 -0.5258181
.sig03         2.7687916  8.7684621
.sigma         1.8734195  2.3471127
(Intercept)   26.7283448 28.9332127
TreatmentDrugx 12.1278192 20.4275981

Warning messages: 1: In optwrap(optimizer, par = start, fn = function(x) dd(mkpar(npar1, : convergence code 1 from bobyqa: bobyqa -- maximum number of function > evaluations exceeded

2: In optwrap(optimizer, par = start, fn = function(x) dd(mkpar(npar1, : convergence code 1 from bobyqa: bobyqa -- maximum number of function evaluations exceeded

When I try:

confint(Model1, parm="beta_")

I do not get the warnings but only get the confidence intervals for the intercept and for the level of treatment.

Are these intercepts okay, or should the previous warnings make me wary to accept those as CIs?

Sorry if this is a dumb question, again, do not have to much in depth statistical knowledge or coding experience.

it probably won't make a huge difference, but you should try adding

control=lmerControl(optCtrl=list(maxfun=10000))

to your function arguments.

The technical post webpages of this site follow the CC BY-SA 4.0 protocol. If you need to reprint, please indicate the site URL or the original address.Any question please contact:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM