Plot ‘signatures’ (loadings of the Y matrix) or (projected) observed and fitted values (Y and X scores) from an “gslcca” object.
1 2 3 4 5 6 7  ## S3 method for class 'gslcca'
plot(x, type = "signature", series = x$treatment,
mean = FALSE, overlay = length(agrep(type, "signature")),
ask = dev.interactive(), lattice = !length(agrep(type, "signature")),
main = NULL, xlab = NULL, ylab = NULL, col = NULL, lty = 1, lwd = 1.5,
pch = NULL, legend.x = "topright", space = "bottom", corner = NULL,
columns = 2, ...)

x 
an object of class 
type 
the type of plot: either 
series 
if 
mean 
logical: if 
overlay 
for 
ask 
logical: if 
lattice 
logical: if 
main 
an overall title for the plot. 
xlab 
a title for the x axis. 
ylab 
a title for the y axis. 
col 
a vector specifying a colour palette to use for the
lines/points. By default the palette is specified by

lty 
a vector of line types, see 
lwd 
a vector of line widths, see 
pch 
a vector of plotting characters, see 
legend.x 
a character vector to use in the legend. 
space 
when 
corner 
when 
columns 
when 
... 
arguments passed on to the plotting functions
( 
For type = "signature"
, the loadings of the Y
matrix from the GSLCCA analysis are displayed for each subject or
their average is displayed, as specified by individual
. When both
individual
and overlay
are TRUE
, the arguments
col
, lty
and legend.x
control the line style and
legend text for each subject.
For type = "fitted"
, the X scores are
displayed with lines. If there are multiple X scores
for each time point, then a line is displayed for each times series
specified by series
(usually the treatment factor used in
GSLCCA) and the arguments col
, lty
,
pch
and legend.x
control the display for each
series. By default, grey or black is used to display values
corresponding to the reference treatment (specified by x$ref
).
For type = "scores"
the X scores are displayed as
for type = "fitted"
and points are added for the Y
scores, again coloured by series
.
Heather Turner
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33  data(clonidine)
### Separate critical exponential model for each treatment
result < gslcca(spectra, "Critical Exponential",
time = Time, treatment = Treatment, subject = Rat,
subject.smooth = 4, data = clonidine)
## plot of individual signatures
plot(result, type = "signature", overlay = FALSE)
## plot of individual projected observed and fitted values
plot(result, type = "scores", lattice = FALSE)
## plot of individual projected fitted values
plot(result, type = "fitted", lattice = FALSE)
## plot of individual signatures  overlaid
plot(result, type = "signature")
## plot of individual signatures  lattice
plot(result, type = "signature", overlay = FALSE, lattice = TRUE)
## plot of individual projected observed and fitted values  lattice
plot(result, type = "scores")
## plot of individual projected fitted values  lattice
plot(result, type = "fitted")
## plot of mean signature
plot(result, type = "signature", mean = TRUE)
## plot of mean fitted
plot(result, type = "fitted", mean = TRUE)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.