| plot.anacor {anacor} | R Documentation |
These functions produce various 2- and 3-dimensional plots for objects of class "anacor"
## S3 method for class 'anacor': plot(x, plot.type, plot.dim = c(1,2), legpos = "top", arrows = FALSE, conf = 0.95, wlines = 0, xlab, ylab, main, type, xlim, ylim, ...) ## S3 method for class 'anacor': plot3d(x, plot.type, plot.dim = c(1,2,3), col.r = "RED", col.c = "BLUE", arrows = TRUE, xlab, ylab, zlab, main, ...) ## S3 method for class 'anacor': plot3dstatic(x, plot.type, plot.dim = c(1,2,3), col.r = "RED", col.c = "BLUE", arrows = TRUE, main, xlab, ylab, zlab, xlim, ylim, zlim, ...)
x |
Object of class "anacor". |
plot.type |
Type of plot to be produced (details see below): 2-D and 3-D for "jointplot", "rowplot", and "colplot"; 2-D for "regplot", "graphplot", "benzplot", and "transplot". |
plot.dim |
Vector of length 2 with Dimensions to be plotted. For "regplot" a single value should be provided, for "transplot" more than two dimensions are allowed. |
legpos |
Position of the legend (for "transplot" only) |
conf |
Ellipsoid confidence level for "jointplot", "rowplot", and "colplot". If NULL, no ellipsoids are drawn. |
arrows |
If TRUE the points in "jointplot", "rowplot", and "colplot" are connected with the origin. |
wlines |
For "graphplot" only: If 0, all lines are of the same thickness. For values > 0 line thickness indicates the strength of the pull |
col.r |
Color for row score plots |
col.c |
Color for column score plots |
xlab |
Label x-axis. |
ylab |
Label y-axis. |
zlab |
Label z-axis. |
xlim |
Scale x-axis. |
ylim |
Scale y-axis. |
zlim |
Scale z-axis. |
main |
Plot title. |
type |
Whether points, lines or both should be plotted; for "regplot" and "transplot" only. |
... |
Additional graphical parameters. |
The following plot types are provided: "jointplot" plots row and column scores into the same device, "rowplot" and "colplot" plot the row scores and column scores, respectively, in separate devices. For these types of plots 3-dimensional versions are provided. The graph plot is an unlabeled version of the joint plot where the points are connected by lines. Options are provided (wlines) to steer the line thickness indicating the connection strength.
The regression plot ("regplot") provides two plots.
First, the unscaled solution is plotted. A frequency grid for the row
categories (x-axis) and column categories (y-axis) is produced.
The regression line is based on the category weighted means of the relative frequencies:
the blue line on the column-wise means on the x-axis and the column category on the y-axis,
the red line is based on the row categories on the x-axis and the row-wise means on the y-axis.
In a second device the scaled solution is plotted. The frequency grid is determined by
the row scores (x-axis) and the column scores(y-axis).
Now, instead of the row/column categories, the column scores (black line y-axis)
and the row scores (red line x-axis) are used.
The transformation plot ("transplot") plots the row/column categories against the row/column scores. The Benzecri plot ("benzplot") plots the observed distances against the fitted distances. It is assumed that the CA result is Benzecri scaled.
Jan de Leeuw, Patrick Mair
de Leeuw, J., & Mair, P. (2007). Simple and canonical correspondence analysis using the R package anacor. Preprint available at http://gifi.stat.ucla.edu/anacor.pdf
##Regression plots using Glass data
data(glass)
res <- anacor(glass)
plot(res, plot.type = "regplot", xlab = "fathers occupation", ylab = "sons occupation")
## Benzecri Plots for bitterling data
data(bitterling)
res1 <- anacor(bitterling, ndim = 2, scaling = c("Benzecri", "Benzecri"))
res2 <- anacor(bitterling, ndim = 5, scaling = c("Benzecri", "Benzecri"))
res2
plot(res1, plot.type = "benzplot", main = "Benzecri Distances (2D)")
plot(res2, plot.type = "benzplot", main = "Benzecri Distances (5D)")
## Column score plot and transformation plot for canonical CA
data(maxwell)
res <- anacor(maxwell$table, row.covariates = maxwell$row.covariates, scaling = c("Goodman", "Goodman"))
res
plot(res, plot.type = "colplot", xlim = c(-1.5,1), conf = NULL)
plot(res, plot.type = "transplot", legpos = "topright")