plot.AMMI {statgenGxE}

Plot function for class AMMI

Description

Two types of biplot can be made. A plot of genotype and environment means vs PC1 (AMMI1) or a biplot of genotypes and environment interaction with PC1 and PC2 (AMMI2).

If the AMMI analysis was done by year, a separate plot will be made for every year in the data. For some years the number of principal components may be lower than the number specified on the secondary axis. If this is the case this year is skipped when plotting. If this happens for all years the function returns an error.

Usage

## S3 method for class 'AMMI'
plot(
  x,
  ...,
  plotType = c("AMMI1", "AMMI2", "GGE2"),
  scale = 1,
  plotGeno = TRUE,
  colorGenoBy = NULL,
  colGeno = NULL,
  sizeGeno = 0,
  plotConvHull = FALSE,
  plotEnv = TRUE,
  colorEnvBy = NULL,
  colEnv = NULL,
  sizeEnv = 3,
  envFactor = 1,
  primAxis = "PC1",
  secAxis = "PC2",
  rotatePC = NULL,
  title = NULL,
  output = TRUE
)

Arguments

x

An object of class AMMI

...

Not used.

plotType

A character string indicating which plot should be made. Either "AMMI1" for an AMMI1 plot (genotype and environment means vs PC1) or "AMMI2" for an AMMI2 biplot (genotypes and environment interaction with PC1 and PC2) respectively. For results of a GGE analysis only an GGE2 biplot can be made and plotType may be ignored.

scale

A numerical value. The variables are scaled by lambda ^ scale and the observations by lambda ^ (1 - scale) where lambda are the singular values computed by princomp in gxeAmmi. Normally 0 <= scale <= 1, and a warning will be issued if the specified scale is outside this range.

plotGeno

Should genotypes be plotted?

colorGenoBy

A character string indicating a column in the TD used as input for the AMMI analysis by which the genotypes should be colored. If NULL all genotypes will be colored in black.

colGeno

A character vector with plot colors for the genotypes. A single color when colorGenoBy = NULL, a vector of colors otherwise.

sizeGeno

An numerical value indicating the text size for plotting the genotypes. Use sizeGeno = 0 for plotting genotypes as points instead of using their names.

plotConvHull

Should a convex hull be plotted around the genotypes. If TRUE a convex hull is plotted. For GGE2 biplots lines from the origin of the plot perpendicular to the edges of the hull are added. Only valid for AMMI2 and GGE2 biplots.

plotEnv

Should environments be plotted?

colorEnvBy

A character string indicating a column in the TD used as input for the AMMI analysis by which the environments should be colored. If NULL all genotypes will be colored in red.

colEnv

A character string with plot colors for the environments. A single color when colorEnvBy = NULL, a vector of colors otherwise.

sizeEnv

An integer indicating the text size for plotting the environments.

envFactor

A positive numerical value giving a factor by which to blow up the environmental scores. Providing a value between 0 and 1 will effectively blow up the genotypic scores.

primAxis

A character string indicating the principal component to be plotted on the primary axis of the AMMI2 plot. Has to be given as "PCn" where n is the number of the principal component.

secAxis

A character string indicating the principal component to be plotted on the secondary axis of the AMMI2 plot. Has to be given as "PCn" where n is the number of the principal component. n Has to differ from primAxis.

rotatePC

A character string indicating a genotype or environment that is to be aligned with the positive x-axis in the plot.

title

A character string used a title for the plot.

output

Should the plot be output to the current device? If FALSE only a list of ggplot objects is invisibly returned.

Value

A biplot depending on plotType. The ggplot object for the biplot is returned invisibly.

See Also

Other AMMI: fitted.AMMI(), gxeAmmi(), report.AMMI(), residuals.AMMI()

Examples

## Run AMMI analysis.
geAmmi <- gxeAmmi(TD = TDMaize, trait = "yld")

## Create an AMMI1 biplot.
plot(geAmmi)

## Create an AMMI2 biplot.
plot(geAmmi, plotType = "AMMI2", scale = 0.5)

## Create an AMMI2 biplot, with HN96b along the positive x-axis.
plot(geAmmi, plotType = "AMMI2", scale = 0.5, rotatePC = "HN96b")

## Run GGE analysis.
geGGE <- gxeGGE(TD = TDMaize, trait = "yld")

## Create an GGE2 biplot.
## Add a convex hull.
plot(geGGE, plotType = "GGE2", scale = 0.5, plotConvHull = TRUE)