Code to Analyze Occupancy Data for Dixie Valley Toads, Anaxyrus williamsi in Churchill County, Nevada from 2018 to 2021

-R code to create and fit dynamic occupancy models to Dixie Valley Toad survey data in JAGS and produce useful summaries of model results. This script will perform the Gibbs Variable Selection (GVS) using the initial set of environmental covariates, and then fit the "final" model using only important environmental covariates. The dynamic occupancy model is adapted from the model published by Duarte et al. (2020).