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

[Ben Bolker]

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
does
http://stackoverflow.com/questions/14358811/extract-prediction-band-from-lme-fit
...

  Ben Bolker