[R-sig-ME] predict GAMM

Rafael Wüest rafael.wueest at gmail.com
Wed Nov 26 11:17:54 CET 2014

Dear ME experts

I'm trying to use GAMMs in order to allow for a correlation structure in a GAM. (A simplified version of) My model looks 
as follows:

m1 <- gamm(y ~ s(p1) + s(p2), correlation = corSymm(form = ~ 1 | p3), data = data.me)

I have two related questions:

1) How do I get the fitted values of this model? I suspect I have to use predict(m1$lme) or fitted(m1$lme)...

2) Is it possible to predict with new data? predict(m1$gam) allows this via a newdata argument, but I suspect the 
correlation structure wouldn't be included then...

Many thanks for any hints,

Rafael Wüest
rafael.wueest at gmail.com

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