I have been working on creating a multivariate predictive model (most likely PLS) using FTIR spectra to predict the amount of a component in complex mixtures. In order to calibrate the model I have considered two pathways - creating synthetic mixtures of the component and confounders or using an established technique to give me the true value of the component in real mixtures. Both pathways have big limitations. I am now considering spiking real mixtures with several known (by weight) amount of the component. How can I implement this in a X-Y chemometric model ?