[R-sig-ME] nlme : how to obtain confidence intervals of fitted values

Ben Bolker bbolker at gmail.com
Fri Oct 31 16:07:33 CET 2014

On 14-10-31 10:50 AM, Laura Riggi wrote:
> Hi
> I used lme (from package nlme) to fit a model containing two fixed
> factors, one random factor and a variance function
> (weights=varIdent(form=~1|as.factor(year_collected)). This is my
> "best" model.
> I would like to plot the confidence intervals for the predicted
> fitted values. Something similar to se.fit=True in predict.lm but for
> lme models containing random effects and variance function.
> Does anybody knows how to do this with R ?
> Thanks Laura

   http://glmm.wikidot.com/faq#predconf has a recipe for doing this, as

  Ben Bolker

More information about the R-sig-mixed-models mailing list