Displaying 1 result from an estimated 1 matches for "totn".
Did you mean:
ton
2011 Dec 15
0
corCompSymm in gamm()?
...o be able to use AR-1 (I assume that I have too little
data for autoregression models of higher orders (ARMA)). If I rename weeks
as week 1-week 30 to get unique identifiers, I loose the seasonal and
year-specific effect in the final model.
M2c.gamm <- gamm(het ~ s(LN.DIN, k=3) + s(LN.totn, k=3) + s(Ncell,
k=3) + s(LN.biom) + s(temp, k=3) + s(week, k=3)
+ fstation + fyear, method = "ML",
weights = varIdent(form=~1 | fstation),
data = data1,
correlation = corCompSymm(form = ~ week|year)) #seems not to work in
a gamm()
Thank you for your time!
Anna Z...