Case studies in Bayesian microbial risk assessments
© Kennedy et al; licensee BioMed Central Ltd. 2009
Published: 21 December 2009
The quantification of uncertainty and variability is a key component of quantitative risk analysis. Recent advances in Bayesian statistics make it ideal for integrating multiple sources of information, of different types and quality, and providing a realistic estimate of the combined uncertainty in the final risk estimates.
We present two case studies related to foodborne microbial risks. In the first, we combine models to describe the sequence of events resulting in illness from consumption of milk contaminated with VTEC O157. We used Monte Carlo simulation to propagate uncertainty in some of the inputs to computer models describing the farm and pasteurisation process. Resulting simulated contamination levels were then assigned to consumption events from a dietary survey. Finally we accounted for uncertainty in the dose-response relationship and uncertainty due to limited incidence data to derive uncertainty about yearly incidences of illness in young children. Options for altering the risk were considered by running the model with different hypothetical policy-driven exposure scenarios. In the second case study we illustrate an efficient Bayesian sensitivity analysis for identifying the most important parameters of a complex computer code that simulated VTEC O157 prevalence within a managed dairy herd. This was carried out in 2 stages, first to screen out the unimportant inputs, then to perform a more detailed analysis on the remaining inputs. The method works by building a Bayesian statistical approximation to the computer code using a number of known code input/output pairs (training runs).
We estimated that the expected total number of children aged 1.5-4.5 who become ill due to VTEC O157 in milk is 8.6 per year, with 95% uncertainty interval (0,11.5). The most extreme policy we considered was banning on-farm pasteurisation of milk, which reduced the estimate to 6.4 with 95% interval (0,11). In the second case study the effective number of inputs was reduced from 30 to 7 in the screening stage, and just 2 inputs were found to explain 82.8% of the output variance. A combined total of 500 runs of the computer code were used.
These case studies illustrate the use of Bayesian statistics to perform detailed uncertainty and sensitivity analyses, integrating multiple information sources in a way that is both rigorous and efficient.
Quantitative Microbial Risk Assessment (QMRA) is a model-based tool for managing food chain risks and enhancing the science behind food safety regulations. QMRA models commonly describe complex farm-to-fork food production chains, involving uncertain parameters and processes. Interest concerns the impact that uncertainty in model inputs and structure can have on uncertainty in model outputs (estimates of "risk"). Bayesian statistics provides a unified approach for handling uncertainty and variability: by contrast with classical approaches it describes uncertainty probabilistically, permits strengthening of inference via incorporation of expert opinions, and synthesizes multiple types of uncertainty in a mathematical framework. Recent developments in Bayesian statistics, such as elicitation of expert opinions , Monte Carlo methods , and sensitivity analyses , have been applied in other fields and have great potential within QMRA.
We present two new case studies as examples. The first quantifies the risk from Vero-cytotoxigenic E. coli (VTEC) O157 in pasteurised milk. The second illustrates an efficient Bayesian sensitivity analysis of a complex simulation model that estimates the prevalence of VTEC O157 within a managed dairy herd.
Case study 1: dietary exposure to VTEC O157 from milk products
To complete the exposure assessment data was taken on daily consumption of milk products of 1858 UK toddlers (aged 1.5-4.5) from the National Diet and Nutrition Survey (NDNS)  to be representative of the UK toddler population. The consumption and contamination models were linked by randomly assigning consumption events in the NDNS to either on- or off-farm contamination levels (cfu/ml), in proportion to the estimated market share of both types. As the market share for on-farm pasteurised milk is very small (around 0.16%), we separately generated results for the entire UK toddler population and for a 'high-risk' subpopulation of toddlers consuming milk only from on-farm pasteurisers. Sufficient milk contamination levels were generated to match the number of consumed items by running the code repeatedly, to represent the distribution of between-item variability. The total daily exposure E i was then calculated for each toddler together with the illness probability p = 1 - (1 + E i /β)-α. Illness events were simulated using a Bernoulli distribution with probability p. Finally, approximate scaling factors were used to scale up from the NDNS-based incidences to the relevant populations.
Bayesian methods were used to quantify uncertainties due to limited information about some of the model inputs. For example, starting from uninformative beta(1,1) prior distributions for many of the rate and proportion inputs listed above, the posterior distributions were also beta distributed, after updating them using a sample of observed 'true' and 'false' outcomes. In each case, the amounts of data used in the updates determined the strength of information in the respective posterior distributions. The implied uncertainty in incidences was calculated via Monte Carlo simulation. The chain of models described above was repeated 1000 times with values for the uncertain inputs generated from their probability distributions. This included the posterior distribution of (α, β) derived in , so that a different relationship between dose and illness probability was used in each iteration.
Management options and how we modelled them
Status Quo. The existing level of risk was estimated, based on current consumption levels and farming practice;
Ban on-farm pasteurisation, so that all consumed milk is assumed to come from off-farm pasteurisers;
Test more post-pasteurisation batches leaving a small dairy, effectively reducing the amount of contaminated milk by an amount depending on the accuracy of the test and the proportion of batches tested;
More training for on-farm pasteurisers. We assume that the ratio of human errors to mechanical errors (ordinarily assumed to be around 4:1) can be reduced through training down to 2:1 or 1:1, altering the off-farm pasteurisation failure rate accordingly. Risk managers could then decide how much reduction is desirable and practical.
Case study 2: sensitivity analysis of a model for VTEC O157 transmission in a UK dairy herd
Our next case study used a deterministic SIR-type (Susceptible/Infective/Recovered) model to simulate the transmission of VTEC O157 within a typical UK dairy herd . To capture the flow of infection around this structured population, the model included 4 interacting sub-groups (unweaned, weaned, dry and lactating), demographic processes and direct, indirect and pseudovertical (occurring in the first 24-48 h after birth) transmission. The result was a complex computer model containing 39 independent parameters. The aim was to identify key parameters related to prevalence of infection and suggest strategies for its reduction. The output of interest was prevalence of VTEC O157 in the weaned group, 150 days after the introduction of infection, as in .
As with many Bayesian methods for analysing computer codes [3, 10] we assume the code output is an unknown function of its inputs (unknown in the Bayesian sense that before the code is run at a given input, the corresponding output is not available to the analyst). A Gaussian process prior distribution is assumed for this function , with a linear mean function that has unknown coefficients and a covariance between any pair of outputs that is a decreasing function of the distance between the corresponding inputs. This covariance represents a belief that the output is a smooth, continuous function of the inputs. After observing a series of selected code input-output (training) runs, the prior is updated to yield a posterior distribution for the function. This posterior is commonly referred to as the Bayesian emulator of the code. It includes the posterior mean, which is an approximation to the function, and a measure of residual uncertainty at any untried input. The emulation uncertainty is zero at the training inputs and increases smoothly at more distant points. Various uncertainty and sensitivity measures can be derived analytically from the emulator, if the input distributions are uniform or normal, without the need for additional code runs. These include variance based sensitivity measures that quantify the relative contribution to the total output variance due to variation in single inputs (main effects) and additional contributions to this variance due to variations in pairs of effects (interaction effects). A two-stage Bayesian sensitivity analysis  was carried out using the GEM-SA software  to investigate which of the uncertainties in the inputs contribute most to uncertainty in the output. The software is limited to 30 independent input parameters, so we first reduced the number of inputs to 30 by constraining a number of parameters to be common between groups (e.g. by making the recovery rate for the weaned group equal to the recovery rate for the unweaned group). In the first stage, 400 code runs were generated with inputs generated using the LP-Tau algorithm in GEM-SA. Since this was only a preliminary screening stage to identify those inputs that might possibly have a non-negligible influence on the output, uniform input distributions were assumed for all 30 inputs, with minimum and maximum values chosen to be conservative, based on the minimum and maximum values specified in . These ranges represent a mixture of uncertainty and between-farm variability. A more detailed second stage sensitivity analysis was then carried out only considering the 'important' screened inputs (those contributing more than 0.5% to the total output variance). At this stage, more realistic distributions for the key inputs should be obtained, as the outputs are most sensitive to these and we would like a realistic assessment of the impact of genuine uncertainty in these inputs. Data analysis or methods for eliciting expert opinions  could be employed to obtain more detailed distributions, but this was not available in our study.
Case study 1: dietary exposure to VTEC O157 from milk products
Risk estimates for Case Study 1
Toddler group & scenario
Estimated incidence (per year)
95% credible intervals
Current UK toddlers
Consumers of on-farm milk products only
Current UK toddlers
Testing (50% reduction)
Training (40% reduction)
Training (60% reduction)
Case study 2: sensitivity analysis of a model for VTEC O157 transmission in a UK dairy herd
Of the initial 30 inputs, only 7 were found to contribute more than 0.5% to the total output variance. This screening of important input parameters meant that a more detailed investigation could be carried out at the second stage. A smaller input design was created, using 100 code runs but now allowing the 'common' inputs to vary independently. For example, the restriction from Stage 1 was relaxed so that the weaned and unweaned group were assumed to have distinct recovery rates. The new emulator was built using only the 8 sensitive inputs with the remaining inputs fixed at their default values. The uniform distributions were unchanged from the first stage. The effects of just 2 inputs, including the interaction effect based on this pair, explained 82.8% of the output variance. The inputs are recovery rate (i.e. the inverse of the duration of infection), which contributed 54%, and concentration of VTEC O157 in faeces, which contributed 22.6%. Both inputs have an impact on the level of indirect transmission (i.e. transmission via the environment). The joint effect contribution was 6.2%. The implication is that efforts to reduce prevalence in the weaned group should focus on changing the values of these parameters. Measures such as vaccination  and dietary management [13, 14] can reduce the level and duration of shedding of VTEC O157 in cattle. Our results suggest that these should be more successful than alternative control measures (e.g. improving hygiene). A more conventional sensitivity analysis of this model, using partial rank correlation coefficients, is reported in , and was based on 1000 runs of the model. Our 2-stage Bayesian approach achieved similar results based on a total of 500 runs.
The case studies presented here have demonstrated some of the potential uses for Bayesian approaches in QMRA to integrate information sources in an efficient, rigorous way. In Case Study 1, multiple sources of uncertainty are propagated through a farm-to-fork model to produce uncertainty in the exposures for the population. This is coupled with uncertainty derived from a Bayesian dose-response model to produce uncertainty intervals for health outcomes that are relevant to policy-makers and other stakeholders (Table 1). Case Study 2 showed how Bayesian sensitivity analysis could be used to understand the influence of uncertain inputs, in a way that is far more efficient than conventional Monte Carlo simulation.
The peer review of this article can be found in Additional file 1.
This work was funded under the joint Environment and Human Health Programme (grant NE/E008526/1), which is supported by the Natural Environment Research Council (NERC), the Environment Agency, the Department for Environment, Food and Rural Affairs (Defra), the Ministry of Defence, Medical Research Council, The Wellcome Trust, Economic and Social Research Council (ESRC), Biotechnology and Biological Sciences Research Council (BBSRC), Engineering and Physical Sciences Research Council and Health Protection Agency. Norval Strachan and Peter Teunis provided information about the dose-response model in Case Study 1. We are also grateful to Clare Bailey, Simon French and Andy Hart from the RELU-RISK team, for their involvement in Case Study 1. RELU-RISK (RES 224-25-0090) is funded jointly by ESRC, BBSRC and NERC, with additional funding from Defra and the Scottish Executive Environment and Rural Affairs Department.
This article has been published as part of Environmental Health Volume 8 Supplement 1, 2009: Proceedings of the Joint Environment and Human Health Programme: Annual Science Day Conference and Workshop. The full contents of the supplement are available online at http://www.ehjournal.net/supplements/8/S1.
- O' Hagan A, Buck CE, Daneshkhah A, Eiser JE, Garthwaite PH, Jenkinson DJ, Oakley JE, Rakow T: Uncertain Judgements: Eliciting Expert Probabilities. 2006, Wiley, ChichesterView ArticleGoogle Scholar
- Robert CP, Casella G: Monte Carlo Statistical Methods. 2006, New York: Springer-Verlag, secondGoogle Scholar
- Oakley JE, O'Hagan A: Probabilistic sensitivity analysis of complex models: a Bayesian approach. Journal of the Royal Statistical Society (Series B). 2004, 66: 751-769. 10.1111/j.1467-9868.2004.05304.x.View ArticleGoogle Scholar
- RELU-RISK. [http://www.relurisk.org.uk]
- Clough HE, Clancy D, French NP: Quantifying exposure to Vero-cytotoxigenic Escherichia coli O157 in milk sold as pasteurized: a model-based approach. International Journal of Food Microbiology. 2009, 131: 95-105. 10.1016/j.ijfoodmicro.2008.12.036.View ArticleGoogle Scholar
- Gregory J, Collins DL, Davies PSW, Hughes JM, Clarke PC: National Diet and Nutrition Survey: children aged 1 1/2 to 4 1/2 years: Report of the diet and nutrition survey. 1995, London: HMSO, 1:Google Scholar
- Teunis PFM, Ogden ID, Strachan NJC: Hierarchical dose response of E. coli O157:H7 from human outbreaks incorporating heterogeneity in exposure. Epidemiol Infect. 2008, 136: 761-770. 10.1017/S0950268807008771.Google Scholar
- Turner J, Begon M, Bowers RG, French NP: A model appropriate to the transmission of a human food-borne pathogen in a multigroup managed herd. Preventative Veterinary Medicine. 2003, 57: 175-198. 10.1016/S0167-5877(03)00006-0.View ArticleGoogle Scholar
- Turner J, Bowers RG, Begon M, Robinson SE, French NP: A semi-stochastic model of the transmission of Escherichia coli O157 in a typical UK dairy herd: Dynamics, sensitivity analysis and intervention/prevention strategies. Journal of Theoretical Biology. 2006, 241: 806-822.View ArticleGoogle Scholar
- O'Hagan A: Bayesian analysis of computer code outputs: a tutorial. Reliability Engineering and System Safety. 2006, 91: 1290-1300. 10.1016/j.ress.2005.11.025.View ArticleGoogle Scholar
- GEM-SA, Version 1.1. Software: Gaussian Emulation Machine for Sensitivity Analysis. [http://www.ctcd.group.shef.ac.uk/gem.html]
- Potter AA, Klashinsky S, Li YL, Frey E, Townsend H, Rogan D, Erickson G, Hinkley S, Klopfenstein T, Moxley RA: Decreased shedding of Escherichia coli O157: H7 by cattle following vaccination with type III secreted proteins. Vaccine. 2004, 22: 362-369. 10.1016/j.vaccine.2003.08.007.View ArticleGoogle Scholar
- Van Baale MJ, Sargeant JM, Gnad DP, DeBey BM, Lechtenberg KF, Nagaraja TG: Effect of forage or grain diets with or without monensin on ruminal persistence and fecal Escherichia coli O157:H7 in cattle. Applied and Environmental Microbiology. 2004, 70: 5336-5342. 10.1128/AEM.70.9.5336-5342.2004.View ArticleGoogle Scholar
- Schamberger GP, Phillips RL, Jacobs JL, Diez-Gonzalez F: Reduction of Escherichia coli O157: H7 populations in cattle by addition of colicin E7-producing E-coli to feed. Applied and Environmental Microbiology. 2004, 70: 6053-6060. 10.1128/AEM.70.10.6053-6060.2004.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.