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
nlme : how to obtain confidence intervals of fitted values
2 messages · Laura Riggi, 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