We used R to run all the analysis \cite{R20} with the added functionality of the tidyverse packages \cite{wickham19}. We wrote the dynamic occupancy model in Stan and used the package rstan to work from R  \cite{rstan20}. We ran four Hamiltonian Monte Carlo chains of 1500 iterations, discarding the first 500 as adaptation.