Plotting the “Mean Spectrum”

February 17, 2012

(This article was first published on NIR-Quimiometría, and kindly contributed to R-bloggers)

Mean spectrum calculation is important:
To center a matrix of spectra, we subtract the mean spectrum, from every spectrum in the matrix.
There are also many options to use the mean spectrum, like average subsamples.
Let´s calculate and plot the mean spectra for the Yarn NIR Data:
> yarn_mean<-colMeans(yarn$NIR)
> wavelength<-c(1:268)
> matplot(wavelength,yarn_mean,lty=1,pch=21,col="red",
  + xlab="data_points",ylab="mean spec.")
We can see in the following plot the X matrix centered:

To leave a comment for the author, please follow the link and comment on his blog: NIR-Quimiometría. offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...

If you got this far, why not subscribe for updates from the site? Choose your flavor: e-mail, twitter, RSS, or facebook...


Comments are closed.