Sunday, September 28, 2025

Risk Assessment Model for Foodborne Pathogens using R statistical programming environment

Dr. Kshitij Shrestha, SFRO, DFTQC

Abstract

This study replicates and extends the FDA-iRISK comparative risk assessment framework, originally described by Chen et al.(2013), within the R statistical programming environment. The primary objective was to reconstruct the core computational models for assessing the public health burden of foodborne hazards, specifically Salmonella in peanut butter and Listeria monocytogenes in soft ripened cheese and cantaloupe, while introducing advanced visualization capabilities not present in the original web-based tool. The R-based model faithfully implements the key components of iRISK: stochastic process models simulating pathogen fate from production to consumption, population-specific consumption patterns, dose-response relationships (Beta-Poisson for Salmonella, exponential for L. monocytogenes), and health impact quantification using Disability-Adjusted Life Years (DALYs). Monte Carlo simulation techniques (mc2d package) were employed to propagate variability and uncertainty. The model successfully reproduced the central findings of the original case studies, confirming the significant annual DALY burden of Salmonella in peanut butter (61.9 DALYs) and the variable risk of listeriosis across perinatal (9.8 DALYs), elderly (5.6 DALYs), and intermediate-age (1.0 DALYs) populations. Beyond replication, this implementation enhances the original iRISK system by generating a suite of dynamic visualizations, including distributional analyses of concentrations, doses, and risks, comparative population risk matrices, and detailed intervention effectiveness plots. These visualizations provide deeper insights into the underlying risk dynamics and the impact of potential interventions, such as pathogen log-reduction and temperature control. The R script offers a transparent, modifiable, and extensible platform for quantitative microbiological risk assessment (QMRA), making it a valuable tool for researchers and risk managers seeking to evaluate and communicate food safety risks. The complete R code is provided for full reproducibility.

Keywords: Quantitative Microbiological Risk Assessment (QMRA), FDA-iRISK, SalmonellaListeria monocytogenes, Monte Carlo Simulation, R Programming, DALYs, Food Safety, Intervention Analysis.



Popular Posts