Bayesian inference and model choice for Holling’s disc equation: a case study on an insect predator-prey system

Papanikolaou, N.E. and Williams, H. and Demiris, Nikolaos and Preston, S.P. and Milonas, P.G. and Kypraios, Theodore (2016) Bayesian inference and model choice for Holling’s disc equation: a case study on an insect predator-prey system. Community Ecology, 17 (1). pp. 71-78. ISSN 1588-2756

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
Download (428kB) | Preview


The dynamics of predator-prey systems relate strongly to the density (in)dependent attributes of the predator’s feeding rate, i.e., its functional response. The outcome of functional response models is often used in theoretical or applied ecology in order to extract information about the mechanisms associated with the feeding behavior of predators. The focus of this study centres upon Holling’s type II functional response model, commonly known as the disc equation, which describes an inverse-density dependent mortality caused by a single predator to its prey. A common method to provide inference on functional response data involves nonlinear least squares optimization, assuming independent Gaussian errors, an assumption often violated in practice due to the heteroscedasticity which is typically present in the data. Moreover, as prey depletion is common in functional response experiments, the differential form of disc equation ought to be used in principle. We introduce a related statistical model and adopt a Bayesian approach for estimating parameters in ordinary differential equation models. In addition, we explore model uncertainty via Bayes factors. Our approach is illustrated via the analysis of several data sets concerning the functional response of a widespread ladybird beetle (Propylea quatuordecimpunctata) to its prey (Aphis fabae), predicting the efficiency of this predator on a common and important aphid species. The results showed that the approach developed in this study is towards a direction for accurate estimation of the parameters that determine the shape of the functional response of a predator without having to make unnecessary assumptions. The R ( code for fitting the proposed model to experimental data is made freely available.

Item Type: Article
Additional Information: This is not the final published version of the paper
Keywords: Functional response, Markov Chain Monte Carlo, Bayes Factors, Ordinary differential equation models, Coccinellidae
Schools/Departments: University of Nottingham, UK > Faculty of Science > School of Mathematical Sciences
Identification Number: 10.1556/168.2016.17.1.9
Depositing User: Eprints, Support
Date Deposited: 06 Jul 2017 10:06
Last Modified: 06 Jul 2017 12:17

Actions (Archive Staff Only)

Edit View Edit View