Supplement 1. R-code used to simulate occupancy data and to fit the IFM naïve, IFM missing, and IFM robust models.

File List 1_MCMC_FUNCTIONS.R 2_IFM_NO_MISSING_MCMC_FUNCTION.R 3_IFM_MISSING_MCMC_FUNCTION.R 4_IFM_ROBUST_MCMC_FUNCTION.R 5_CHAIN_DIAGNOSTICS_FUNCTIONS.R. 6_IFM_SIM_DATA_PREP.R 7_AUDIT_INM.R 8_AUDIT_IFM_MISSING.R 9_AUDIT_IFM_ROBUST.R Description The files in this Supplement are used to fit the three models described in this study (IFM Naïve, IFM Missing, and IFM Robust) and include a program that simulates occupancy data following the Incidence Function Model. The file 1_MCMC_FUNCTIONS.R contains four auxiliary functions used in 2_IFM_NO_MISSING_MCMC_FUNCTION.R, 3_IFM_MISSING_MCMC_FUNCTION.R, and 4_IFM_ROBUST_MCMC_FUNCTION.R. It includes a function that accepts or rejects proposed values in the Metropolis-Hastings algorithm, a function that counts the acceptance rates, a function that counts the number of missing values, and a function that is used when proposing values from a bivariate normal distribution for the parameters gamma and beta. The file 2_IFM_NO_MISSING_MCMC_FUNCTION.R contains the function that uses MCMC to estimate the IFM Naïve. The methods are described in Appendix A. The file 3_IFM_MISSING_MCMC_FUNCTION.R contains the function that uses MCMC to estimate the IFM Missing. The methods are described in Appendix A. The file 4_IFM_ROBUST_MCMC_FUNCTION.R contains the function that uses MCMC to estimate the IFM Robust. It implements the IFM Robust, and the methods are described in Appendix A. The file 5_CHAIN_DIAGNOSTICS_FUNCTIONS.R contains a number of functions used to diagnose the convergence of MCMC chains. In particular, the function coda.create converts the files created by the MCMC functions to a format that can be read by the “coda” R-package. The file 6_IFM_SIM_DATA_PREP.R simulates occupancy data from the Incidence Function Model. These data sets are then used in the subsequent programs. The file 7_AUDIT_INM.R estimates the parameters of the IFM Naïve using data created in 6_IFM_SIM_DATA_PREP.R. The file 8_AUDIT_IFM_MISSING.R estimates the parameters of the IFM Missing using data created in 6_IFM_SIM_DATA_PREP.R. The file 9_AUDIT_IFM_ROBUST.R estimates the parameters of the IFM Robust using data created in 6_IFM_SIM_DATA_PREP.R.