## generation of the data
data(tea)
new.lambdas <- seq(260, 500, by = 2)
tea <- lapply(tea.raw,
preprocess,
dim2 = new.lambdas)
tea.split <- splitTimeWindow(tea, c(12, 14), overlap = 10)
tea.alslist <- lapply(tea.split,
function(Xl) {
Xl.opa <- opa(Xl, 4)
doALS(Xl, Xl.opa)
})
teaMerged <- mergeTimeWindows(tea.alslist)
## This is the object saved in teaMerged.RData
ncomp <- ncol(teaMerged$S)
myPalette <- colorRampPalette(c("black", "red", "blue", "green"))
mycols <- myPalette(ncomp)
plot(teaMerged, what = "spectra", col = mycols)
legend("top", col = mycols, lty = 1, bty = "n", ncol = 2,
legend = paste("C", 1:ncol(teaMerged$S)))
Run the code above in your browser using DataLab