updating code to get results and figures for talk with new data
added gamm to look for variance explained by jday with year as random effect