# NOT RUN {
if (require("ChemometricsWithRData")) {
data(bdata, package = "ChemometricsWithRData")
D1.efa <- efa(bdata$d1, 3)
matplot(D1.efa$forward, type = "l")
matplot(D1.efa$backward, type = "l")
matplot(D1.efa$pure.comp, type = "l")
D1.opa <- opa(bdata$d1, 3)
matplot(D1.opa$pure.comp, type = "l")
D1.mcr.efa <- mcr(bdata$d1, D1.efa$pure.comp, what = "col")
matplot(D1.mcr.efa$C, type = "l", main = "Concentration profiles")
matplot(t(D1.mcr.efa$S), type = "l", main = "Pure spectra")
}
# }
Run the code above in your browser using DataLab