Second-order modeling of variability and uncertainty in microbial hazard characterization

This study describes an analytical framework that permits quantitative consideration of variability and uncertainty in microbial hazard characterization. Second-order modeling that used two-dimensional Monte Carlo simulation and stratification into homogeneous population subgroups was applied to integrate uncertainty and variability. Specifically, the bootstrap method was used to simulate sampling error due to the limited sample size in microbial dose-response modeling. A data set from human feeding trials with Campylobacter jejuni was fitted to the log-logistic dose-response model, and results from the analysis of FoodNet surveillance data provided further information on variability and uncertainty in Campylobacter susceptibility due to the effect of age. Results of our analyses indicate that uncertainty associated with dose-response modeling has a dominating influence on the analytical outcome. In contrast, inclusion of the age factor has a limited impact. While the advocacy of more closely modeling variability in hazard characterization is warranted, the characterization of key sources of uncertainties and their consistent propagation throughout a microbial risk assessment actually appear of greater importance.