A latent variable model for identifying and characterizing food adulteration

Recently, growing consumer awareness of food quality and sustainability has led to a rising demand for effective food authentication methods. Vibrational spectroscopy techniques have emerged as a promising tool for collecting large volumes of data to detect food adulteration. However, spectroscopic data pose significant challenges from a statistical viewpoint, highlighting the need for more sophisticated modeling strategies. To address these challenges, in this work we propose a latent variable model specifically tailored for food adulterant detection, while accommodating the features of spectral data. Our proposal offers greater granularity with respect to existing approaches, since it does not only identify adulterated samples but also estimates the level of adulteration, and detects the spectral regions most affected by the adulterant. Consequently, the methodology offers deeper insights, and could facilitate the development of portable and faster instruments for efficient data collection in food authenticity studies. The method is applied to both synthetic and real honey mid-infrared spectroscopy data, delivering precise estimates of the adulteration level and accurately identifying which portions of the spectra are most impacted by the adulterant.

💡 Research Summary

The paper addresses the growing demand for reliable food authentication methods driven by heightened consumer awareness of quality and sustainability. While vibrational spectroscopy techniques such as near‑infrared (NIR) and mid‑infrared (MIR) provide rapid, non‑destructive acquisition of high‑dimensional spectral data, the statistical analysis of such data remains challenging due to the large number of wavelengths and strong inter‑wavelength correlations. Existing chemometric approaches (e.g., PCA‑LDA, PLS‑DA, or standard mixture models) often yield coarse binary decisions and struggle to quantify low levels of adulteration.

To overcome these limitations, the authors formulate the problem within the framework of individual‑level mixture models. They model each observed spectrum yᵢ (a p‑dimensional vector) as a Gaussian random vector with mean μ_P + gᵢ δ, where μ_P is the mean spectrum of the pure food, δ is the mean‑shift vector representing the difference between pure food and pure adulterant spectra, and **gᵢ ∈