The seasonal variation in fatty acids and minerals concentrations was investigated through the analysis of Pecorino Romano cheese samples collected in January, April, and June. A fraction of samples contained missing values in their fatty acid profile. Probabilistic Principal component analysis coupled with Linear Discriminant Analysis was employed to classify cheese samples on a production season basis while accounting for missing data and quantifying the missing Fatty acids concentration for the sample in which they were absent. The levels of rumenic acid, vac-cenic acid and omega-3 compounds were positively correlated with the spring season, while the length of the saturated fatty acids increased throughout the production seasons. Concerning the classification performances, the optimal number of principal components (i.e., 5) achieved an ac-curacy in cross-validation equal to 98 %. Then, when the model was tasked to impute the lacking Fatty acid concentration values, the optimal number of principal components resulted in an R2 value in cross-validation of 99.53%