A Generalized QMRA Beta-Poisson Dose-Response Model

Abstract

Quantitative microbial risk assessment (QMRA) is widely accepted for characterizing the microbial risks associated with food, water, and wastewater. Single-hit dose-response models are the most commonly used dose-response models in QMRA. Denoting (Formula presented.) as the probability of infection at a given mean dose d, a three-parameter generalized QMRA beta-Poisson dose-response model, (Formula presented.), is proposed in which the minimum number of organisms required for causing infection, Kmin, is not fixed, but a random variable following a geometric distribution with parameter (Formula presented.). The single-hit beta-Poisson model, (Formula presented.), is a special case of the generalized model with Kmin = 1 (which implies (Formula presented.)). The generalized beta-Poisson model is based on a conceptual model with greater detail in the dose-response mechanism. Since a maximum likelihood solution is not easily available, a likelihood-free approximate Bayesian computation (ABC) algorithm is employed for parameter estimation. By fitting the generalized model to four experimental data sets from the literature, this study reveals that the posterior median (Formula presented.) estimates produced fall short of meeting the required condition of (Formula presented.) = 1 for single-hit assumption. However, three out of four data sets fitted by the generalized models could not achieve an improvement in goodness of fit. These combined results imply that, at least in some cases, a single-hit assumption for characterizing the dose-response process may not be appropriate, but that the more complex models may be difficult to support especially if the sample size is small. The three-parameter generalized model provides a possibility to investigate the mechanism of a dose-response process in greater detail than is possible under a single-hit model.