Altered the model so that it estimates year-specific covariate effects and baseline occurrence
sorted the dataset by year, basically