Bayesian Analysis

Intrinsic Bayesian Analysis for Occupancy Models

Daniel Taylor-Rodríguez, Andrew J. Womack, Claudio Fuentes, and Nikolay Bliznyuk

Full-text: Open access


Occupancy models are typically used to determine the probability of a species being present at a given site while accounting for imperfect detection. The survey data underlying these models often include information on several predictors that could potentially characterize habitat suitability and species detectability. Because these variables might not all be relevant, model selection techniques are necessary in this context. In practice, model selection is performed using the Akaike Information Criterion (AIC), as few other alternatives are available. This paper builds an objective Bayesian variable selection framework for occupancy models through the intrinsic prior methodology. The procedure incorporates priors on the model space that account for test multiplicity and respect the polynomial hierarchy of the predictors when higher-order terms are considered. The methodology is implemented using a stochastic search algorithm that is able to thoroughly explore large spaces of occupancy models. The proposed strategy is entirely automatic and provides control of false positives without sacrificing the discovery of truly meaningful covariates. The performance of the method is evaluated and compared to AIC through a simulation study. The method is illustrated on two datasets previously studied in the literature.

Article information

Bayesian Anal., Volume 12, Number 3 (2017), 855-877.

First available in Project Euclid: 9 September 2016

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

imperfect detection intrinsic priors model priors strong heredity Bayesian variable selection AIC

Creative Commons Attribution 4.0 International License.


Taylor-Rodríguez, Daniel; Womack, Andrew J.; Fuentes, Claudio; Bliznyuk, Nikolay. Intrinsic Bayesian Analysis for Occupancy Models. Bayesian Anal. 12 (2017), no. 3, 855--877. doi:10.1214/16-BA1014.

Export citation


  • Akaike, H. (1983). Information measures and model selection. Bulletin of the International Statistical Institute, 50:277–290.
  • Albert, J. H. and Chib, S. (1993). Bayesian-analysis of binary and polychotomous response data. Journal of the American Statistical Association, 88(422):669–679.
  • Berger, J. and Pericchi, L. (1996). The intrinsic Bayes factor for model selection and prediction. Journal of the American Statistical Association, 91(433):109–122.
  • Burnham, K. and Anderson, D. (2003). Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. Springer New York.
  • Burnham, K. P. (2004). Multimodel Inference: Understanding AIC and BIC in Model Selection. Sociological Methods & Research, 33(2):261–304.
  • Casella, G. and Moreno, E. (2006). Objective Bayesian variable selection. Journal of the American Statistical Association, 101(473):157–167.
  • Chib, S. (1995). Marginal likelihood from the Gibbs output. Journal of the American Statistical Association, 90(432):1313–1321.
  • Chib, S. and Jeliazkov, I. (2001). Marginal likelihood from the metropolis–hastings output. Journal of the American Statistical Association, 96(453):270–281.
  • Chipman, H. (1996). Bayesian variable selection with related predictors. Canadian Journal of Statistics, 24(1):17–36.
  • Dorazio, R. M. and Rodríguez, D. T. (2012). A Gibbs sampler for Bayesian analysis of site-occupancy data. Methods in Ecology and Evolution, 3(6):1093–1098.
  • Fiske, I. and Chandler, R. (2011). unmarked: An R package for fitting hierarchical models of wildlife occurrence and abundance. Journal of Statistical Software, 43(10).
  • Flegal, J. M. and Gong, L. (2013). Relative fixed-width stopping rules for markov chain monte carlo simulations. arXiv:1303.0238.
  • Green, P. J. (1995). Reversible Jump Markov Chain Monte Carlo Computation and Bayesian Model Determination. Biometrika, 82(4):711–732.
  • Guillera-Arroita, G., Lahoz-Monfort, J. J., MacKenzie, D. I., Wintle, B. A., and McCarthy, M. A. (2014). Ignoring imperfect detection in biological surveys is dangerous: A response to: Fitting and interpreting occupancy models’. PLoS ONE, 9(7):e99571.
  • Hooten, M. B. and Hobbs, N. T. (2015). A guide to Bayesian model selection for ecologists. Ecological Monographs, 85(1):3–28.
  • Hurvich, C. M. and Tsai, C.-L. (1989). Regression and time series model selection in small samples. Biometrika, 76:297–307.
  • Kass, R. E. and Wasserman, L. (1996). The selection of prior distributions by formal rules. Journal of the American Statistical Association, 91(435):1343.
  • Kery, M., Gardner, B., and Monnerat, C. (2010). Predicting species distributions from checklist data using site-occupancy models. Journal of Biogeography, 37(10):1851–1862. Kery, Marc Gardner, Beth Monnerat, Christian.
  • Kéry, M., Royle, J. A., and Schmid, H. (2005). Modeling avian abundance from replicated counts using binomial mixture models. Ecological Applications, 15(4):1450–1461.
  • Kéry, M. and Schmid, H. (2004). Monitoring programs need to take into account imperfect species detectability. Basic and Applied Ecology, 5(1):65–73.
  • Leon-Novelo, L., Moreno, E., and Casella, G. (2012). Objective Bayes model selection in probit models. Statistics in Medicine, 31(4):353–65.
  • Liu, J. S. and Wu, Y. N. (1999). Parameter expansion for data augmentation. Journal of the American Statistical Association, 94(448):1264–1274.
  • MacKenzie, D. I., Nichols, J. D., Lachman, G. B., Droege, S., Royle, J. A., and Langtimm, C. A. (2002). Estimating site occupancy rates when detection probabilities are less than one. Ecology, 83(8):2248–2255.
  • Mazerolle, M. and Mazerolle, M. (2013). Package ‘AICcmodavg’. (c).
  • McQuarrie, A., Shumway, R., and Tsai, C.-L. (1997). The model selection criterion AICu. Statistics & Probability Letters, 34(3):285–292.
  • Moreno, E., Bertolino, F., and Racugno, W. (1998). An intrinsic limiting procedure for model selection and hypotheses testing. Journal of the American Statistical Association, 93(444):1451–1460.
  • Peixoto, J. L. (1987). Hierarchical variable selection in polynomial regression models. American Statistician, 41(4):311–313.
  • Peixoto, J. L. (1990). A property of well-formulated polynomial regression-models. American Statistician, 44(1):26–30.
  • Pérez, J. and Berger, J. (2002). Expected-posterior prior distributions for model selection. Biometrika, 491–511.
  • Rao, C. R. and Wu, Y. (2001). On model selection, volume 38 of Lecture Notes–Monograph Series, pages 1–57. Institute of Mathematical Statistics, Beachwood, OH.
  • Roy, V. and Hobert, J. P. (2007). Convergence rates and asymptotic standard errors for Markov chain Monte Carlo algorithms for Bayesian probit regression. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 69(4):607–623.
  • Taylor-Rodriguez, D., Womack, A., and Bliznyuk, N. (2016). Bayesian Variable selection on model spaces constrained by heredity conditions. Journal of Computational and Graphical Statistics 25(2):515–535.
  • Taylor-Rodríguez, D., Womack, A. J., Fuentes, C., and Bliznyuk, N. (2016). “Supplementary Appendices of “Intrinsic Bayesian Analysis for Occupancy Models”.” Bayesian Analysis.
  • Wasserman, L. (2000). Bayesian model selection and model averaging. Journal of Mathematical Psychology, 44(1):92–107.
  • Yang, Y. (2005). Can the strengths of AIC and BIC be shared? Biometrika, 92(4):937–950.

Supplemental materials