### Apply MSC to spectra from simdata
library(mdatools)
data(simdata)
spectra = simdata$spectra.c
wavelength = simdata$wavelength
res = prep.msc(spectra)
cspectra = res$cspectra
par(mfrow = c(2, 1))
mdaplot(cbind(wavelength, t(spectra)), type = 'l', main = 'Before MSC')
mdaplot(cbind(wavelength, t(cspectra)), type = 'l', main = 'After MSC')
Run the code above in your browser using DataLab