[R-sig-ME] confidence intervals from lmer() with random correlation parameters

[Sol Lago]

Hi there,

Lately I keep encountering the same problem and I'm wondering whether other people have been able to get around it. I'm running a mixed effects model using lmer(). My model has by-subject and by-item intercepts and slopes, and random correlation parameters between them. Since the current version of lmer() does not have MCMC sampling implemented, I cannot use pvals.fnc(). I get this message:

> Error in pvals.fnc(m, withMCMC = T) : MCMC sampling is not implemented in recent versions of lme4 for models with random correlation parameters

pvals.fnc() is also the function I always use to get confidence intervals (HPD95lower and HPD95upper were two columns in the pvals.fnc output). Does anyone know of an alternative way of getting confidence intervals for the fixed effects estimates in the model? Or does using models with random correlations means that we can no longer get CIs from R?

Thanks!

--Sol