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) ......+pollen01
,correlation=corARMA(p=7,q=0,form=~Trend|Annees),
niter=40,control=lmeControl(maxIter=100000,msMaxIter=100000)
,family=quasipoisson, data=donnees,
na=na.omit)
Each model need 1/2 hour to 1 hour to get a result and I frequently get the
following error messages:
Erreur dans lme.formula(fixed = fixed, random = random, data = data,
correlation = correlation, :
nlminb problem, convergence error code = 1
message = singular convergence (7)
Erreur dans na.fail.default(list(Xr.4 = c(0.00374377951791214,
0.00438373736920182, :
missing values in object
Does my computer lack power,(if so,what is required?) or is R limited in
speed of execution or is my syntax not ok. Here are the features of my computer:
processor: 3.10 GHz i5-2400 CPU
RAM: 4 Go (3.87 Go)
OS: 64 bits
Windows 7 professionnel (Pack 1)
Another question: how can we add a smoothing parameter in the model GAMM to
help convergence ("sp =" is not accepted)?
Thanks for your help
Best regards
Sylvie Martin
SEPIA-Santé
Bureau d'études en épidémiologie,
biostatistiques, environnement
31 rue de Pontivy
56 150 BAUD
FRANCE
tél : 02 97 28 80 38
fax: 02 97 28 81 10
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
