search for: station4

...cv). 1. In have, say 5 time series. Monthly data, 20 year. The 5 time series are from 5 stations. The data are in vectors, so I have fitted something along the lines of: tmp<-gamm(Y ~ s(Year,by=station1)+s(Year,by=station2)+ s(Year,by=station3)+s(Year,by=station4)+ s(Year,by=station5)+ factor(station)*factor(month), correlation=corAR1(form=~MyTime|station), famliy=gaussian) station is just a long vector with ones, twos, threes, fours and fiv...