I read in an article that Ian Cowe said that what normally chemometricians do is to look to the graphics, of course interpret those graphics. So I still go on trying to develop a function can help me to understand the graphics and all the statistics there are behind.

I add some more lines to the monitor function:

plot(x~y,main=”X-Y plot”,xlab=”predicted”,ylab=”reference”)

abline(0,1,col=”blue”)

abline(intercept,slope,col=”red”)

abline(intercept+(2.5*sep),slope,col=”red”,lty=4)

abline(intercept-(2.5*sep),slope,col=”red”,lty=4)

I can do the same for the residual plot.

There are two warning lines which advice if the residual exceeds 2,5*SEP value. That is the T value warning limit.

Another line can be add, called the T value action limit (3*SEP).

Graphics show the 0-1 abline (blue) and the calculated slope-intercept abline (red). Limits are with dashed red lines.

We can see that almost both lines red and blue are almost over-plotted in this case.

*Related*

To

**leave a comment** for the author, please follow the link and comment on their blog:

** NIR-Quimiometría**.

R-bloggers.com offers

**daily e-mail updates** about

R news and

tutorials on topics such as:

Data science,

Big Data, R jobs, 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...

**Tags:** "R" Chemometrics