PLOT.LRA {easyCODA}

Plot the results of a logratio analysis

Description

Various maps and biplots of the results of a logratio analysis using function LRA.

Usage

PLOT.LRA(obj, map="symmetric", rescale=1, dim=c(1,2), axes.inv = c(1,1), main="",
         cols=c("blue","red"), colarrows = "pink", cexs=c(0.8,0.8), fonts=c(2,4))

Arguments

obj

An LRA object created using function LRA

map

Choice of scaling of rows and columns: "symmetric" (default), "asymmetric" or "contribution"

rescale

A rescaling factor applied to column coordinates (default is 1 for no rescaling)

dim

Dimensions selected for horizontal and vertical axes of the plot (default is c(1,2))

axes.inv

Option for reversing directions of horizontal and vertical axes (default is c(1,1) for no reversing, change one or both to -1 for reversing)

main

Title for plot

cols

Colours for row and column labels (default is c("blue","red"))

colarrows

Colour for arrows in asymmetric and contribution biplots (default is "pink")

cexs

Character expansion factors for row and column labels (default is c(0.8,0.8))

fonts

Fonts for row and column labels (default is c(2,4))

Details

The function PLOT.LRA makes a scatterplot of the results of a logratio analysis (computed using function LRA), with various options for scaling the results and changing the direction of the axes. By default, dimensions 1 and 2 are plotted on the horizontal and vertical axes, and it is assumed that row points refer to samples and columns to compositional parts.

By default, the symmetric scaling is used, where both rows and columns are in principal coordinates and have the same amount of weighted variance along the two dimensions. The other options are the asymmetric option, when columns are in standard coordinates, and the contribution option, when columns are in contribution coordinates. In cases where the row and column displays occupy widely different extents, the column coordinates can be rescaled using the rescale option.

Author(s)

Michael Greenacre

References

Greenacre, M. (2013), Contribution biplots, Journal of Computational and Graphical Statistics, 22, 107-122.

See Also

plot.ca

Examples

# perform LRA on the Roman cups data set and plot the results
data(cups)
cups.LRA <- LRA(cups)
PLOT.LRA(cups.LRA, map="contribution", rescale=0.2)