spOccupancy

Fits single-species, multi-species, and integrated non-spatial and spatial occupancy models using Markov Chain Monte Carlo (MCMC). Models are fit using Polya-Gamma data augmentation detailed in Polson, Scott, and Windle (2013) doi:10.1080/01621459.2013.829001. Spatial models are fit using either Gaussian processes or Nearest Neighbor Gaussian Processes (NNGP) for large spatial datasets. Details on NNGP models are given in Datta, Banerjee, Finley, and Gelfand (2016) doi:10.1080/01621459.2015.1044091 and Finley, Datta, and Banerjee (2020) . Provides functionality for data integration of multiple single-species occupancy data sets using a joint likelihood framework. Details on data integration are given in Miller, Pacifici, Sanderlin, and Reich (2019) doi:10.1111/2041-210X.13110. Details on single-species and multi-species models are found in MacKenzie, Nichols, Lachman, Droege, Royle, and Langtimm (2002) doi:10.1890/0012-9658(2002)083[2248:ESORWD]2.0.CO;2 and Dorazio and Royle doi:10.1198/016214505000000015, respectively.