Produces a scatter plot of the correlation of the variables against their covariance for a chosen principal component. It allows visual identification of variables driving the separation and thus is a useful adjunct to traditional loading plots.
sPlotSpectra(spectra, pca, pc = 1, tol = 0.05, ...)
An object of S3 class
An integer specifying the desired pc plot.
A number describing the fraction of points to be labeled.
tol = 1.0 labels all the points;
tol = 0.05 labels
approximately the most extreme 5 percent. Set to
completely suppress labels. Note that a simple approach
based upon quantiles is used, assumes that both x and y are each normally
distributed, and treats x and y separately. Thus, this is not a formal
treatment of outliers, just a means of labeling points. Groups are lumped
Parameters to be passed to the plotting routines. Applies to base graphics only.
The returned value depends on the graphics option selected (see
None. Side effect is a plot.
The plot is displayed, and a
ggplot2 plot object is returned if the
value is assigned. The plot can be modified in the usual
Wiklund, Johansson, Sjostrom, Mellerowicz, Edlund, Shockcor, Gottfries, Moritz, and Trygg. "Visualization of GC/TOF-MS-Based Metabololomics Data for Identification of Biochemically Interesting Compounds Usings OPLS Class Models" Analytical Chemistry Vol.80 no.1 pgs. 115-122 (2008).
Additional documentation at https://bryanhanson.github.io/ChemoSpec/
# This example assumes the graphics output is set to ggplot2 (see ?GraphicsOptions). library("ggplot2") data(SrE.IR) pca <- c_pcaSpectra(SrE.IR) myt <- expression(bolditalic(Serenoa) ~ bolditalic(repens) ~ bold(IR ~ Spectra)) p <- sPlotSpectra(spectra = SrE.IR, pca = pca, pc = 1, tol = 0.001) p <- p + ggtitle(myt) p