In R, I fitted a dynamic occupancy model using the unmarked package. The best model included one interaction term between a site covariate and a
unmarked
site covariate