Designing optimal food safety monitoring schemes using Bayesian network and integer programming: The case of monitoring dioxins and DL-PCBs.

Efficient food safety monitoring should achieve optimal resource allocation. In this article, a methodology is presented to optimize the use of resources for food safety monitoring aimed at identifying noncompliant samples and estimating background level of hazards in food products. A Bayesian network (BN) model and an optimization model were combined in a single framework. The framework was applied to monitoring dioxins and dioxin-like polychlorinated biphenyls (DL-PCBs) in primary animal-derived food products in the Netherlands. The BN model was built using a national dataset with monitoring results of dioxins and DL-PCBs in animal-derived food products over a 10-year period (2008-2017). These data were used to estimate the probability of detecting suspect samples with dioxins and DL-PCBs levels above preset thresholds, given certain sample conditions. The results of the BN model were then inserted into the optimization model to compute an optimal monitoring scheme. Model estimates showed that the probability of dioxins and DL-PCBs exceeding threshold limits was higher in laying hen eggs and sheep meat than in other animal-derived food (except deer meat). Compared with the monitoring scheme used in the Netherlands in 2018, the optimal monitoring scheme would save around 10,000 EUR per year. This could be obtained by reallocating monitoring resources from products with lower probability of dioxin and DL-PCBs exceeding threshold limits (e.g., pig meat) to products with higher probability (e.g., bovine animal meat), and by shifting sample collection from the last quarter of the year toward the first three quarters of the year.