- For each trajectory we sampled an initial population size and life history parameters from the posterior distribution of the dynamic occupancy model.
- Within each trajectory, and using a step size of one year, we update the number of harriers in the different age categories according to the parameters sampled in step 1.
- We compute population quantiles (90%, 75%, 50%, 25%, 10%) at each time step across population trajectories.