To help me understand the model and fix bugs in my code, I'm building the model one layer at a time and running it on simulated data. Yesterday I added multiple years and individual effects to the simulated data so I could add individual effects to the model. I ran a mixed model with lme4 in R just to try and it works better than a model without individual and year effects - most especially it's better at estimating the steepness of the curve.