I fit a GAM to turtle growth data using mgcv: >m1 <- gam(growth~s(mean.size, bs="cr")+s(year,bs="cr",k=7)+s(cohort,bs="cr")+s(age,bs="cr"), data=grow, family=quasi(link="identity"))
The errors are skewed (and seem to be correlated with age) (code and plots here: https://docs.google.com/fileview?id=0B1cQ7z9xYFl2OWY2MzNlYzUtMWZmMS00YjExLTk0MjMtMGU3MmRhYzI5N2Vm&hl=en). I tried fitting a oversimplified GAMM, >mm1 <- gamm(growth~s(mean.size, bs="cr")+s(year,bs="cr",k=7)+s(cohort,bs="cr")+s(age,bs="cr"), data=grow, family=quasi(link="identity")) which dealt with some of the correlation... but... I don't know what/how to specify a structure for the errors. Any help/suggestions greatly appreciated, Thank you, Claudia [[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.
