However, depending on value of degree of freedom (K) chosen , the result
I run the following code to fit the GAM
fit <- gamm(Mydata ~ s(Age, k =25), correlation = corARMA(form = ~ Age, p =
The package documentation provides some support in "choose K" section by
exploring p-value and k-index (gam.check function), however, I have gotten
no p-value and K-index either.
Could anyone help me to find the best degree of freedom?
Thank you in advance for any support.
Jackson M. Rodrigues
Department of Palynology and Climate Dynamics
Albrecht-von-Haller-Institute for Plant Sciences
Untere Karspuele 2
Tel.: 0049 (0) 176 8186 4994