added gamm analysis as a better alternative to nlme so that smoothing can be asymmetric on julian day or day of the year