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

[Steve Walker]

Hi Sol,

Ben Bolker has put together a help file (?pvalues) for reference about 
p-values and confidence intervals in the new lme4.

Hope this helps,
Steve

On 2013-10-19 2:02 PM, Sol Lago wrote:
> 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
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>