(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.")**

To

**leave a comment**for the author, please follow the link and comment on his blog:**NIR-Quimiometría**.R-bloggers.com 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...