Re: [R] What computer power for GAMM models

2013-10-25 Thread Gavin Simpson
1) The model is far too complex. I have trouble imagining a system where the trend potentially uses 49 degrees of freedom and yet the unexplained variance within years follows an AR(7). It looks like you have time points within years over many years; is that why you are trying to use such a complex

Re: [R] What computer power for GAMM models

2013-10-25 Thread Bert Gunter
... and, I should have added, your maxiter parameters in the call are set way too high. Lowering them to something more sensible will also get your error messages faster. -- Bert On Fri, Oct 25, 2013 at 6:19 AM, Sylvie Martin wrote: > > > Hello, > > I use GAMM function (MGCV package) in R softw

Re: [R] What computer power for GAMM models

2013-10-25 Thread Bert Gunter
Sylvie: I do not know whether your system lacks sufficient power, but its user -- you -- may. The models you are trying to fit are complex and nonlinear, and possibly depending on your data -- do you have a lot of missingness? How large are your data sets? Can they support your models -- you may

[R] What computer power for GAMM models

2013-10-25 Thread Sylvie Martin
Hello, I use GAMM function (MGCV package) in R software to study relationship between pollen and pollinosis. My models include autoregressive terms. Here is an exemple: CupAR7_2<-gamm(nb_rca_2~s(Trend,bs="ps",k=49) +Semaine             +s(TX03,bs="ps",k=8) +s(UN03,bs="ps",k=4)  ..+polle