In this study, we introduce enzymatic perturbation combined with Fourier transform infrared (FTIR) spectroscopy as a concept for quantifying casein in subcritical heated skim milk using chemometric multiway analysis. Chymosin is a protease that cleaves specifically caseins. As a result of hydrolysis, all casein proteins clot to form a creamy precipitate, and whey proteins remain in the supernatant. We monitored the cheese-clotting reaction in real time using FTIR and analyzed the resulting evolution profiles to establish calibration models using parallel factor analysis and multiway partial least squares regression. Because we observed casein-specific kinetic changes, the retrieved models were independent of the chemical background matrix and were therefore robust against possible covariance effects. We tested the robustness of the models by spiking the milk solutions with whey, calcium, and cream. This method can be used at different stages in the dairy production chain to ensure the quality of the delivered milk. In particular, the cheese-making industry can benefit from such methods to optimize production control.