Last chance! 50% off unlimited learning
Sale ends in
Biplot AMMI.
# S3 method for AMMI
plot(x,first=1,second=2,third=3,type=1,number=FALSE,gcol=NULL,ecol=NULL,
icol=NULL,angle=25,lwd=1.8,length=0.1,xlab=NULL,ylab=NULL,xlim=NULL,ylim=NULL,...)
object AMMI
position axis x, 0=Y-dependent, 1=PC1, 2=PC2, 3=PC3
position axis y,0=Y-dependent, 1=PC1, 2=PC2, 3=PC3
position axis z,0=Y-dependent, 1=PC1, 2=PC2, 3=PC3
1=biplot, 2= triplot 3=influence genotype
TRUE or FALSE names or number genotypes
genotype color
environment color
influence color
angle from the shaft of the arrow to the edge of the arrow head
parameter line width in function arrow
parameter length in function arrow
x labels
y labels
x limites
y limites
other parameters of plot
type=1 produce graphs biplot. type=2 produce graphs triplot, the components are normalizad in scale 0-1. type=3 produce graphs on a 2d point set that are all subgraphs of the Delaunay triangulation. The Gabriel graph is a subgraph of the delaunay triangulation and has the relative neighbor graph as a sub-graph.
help(graphneigh) package=spdep
# NOT RUN {
library(agricolae)
data(plrv)
model<- with(plrv,AMMI(Locality, Genotype, Rep, Yield))
# biplot PC2 vs PC1
plot(model)
## plot PC1 vs Yield
plot(model,0,1,gcol="blue",ecol="green")
## triplot PC 2,3,4
if (requireNamespace("klaR", quietly = TRUE)) {
plot(model,first=2,second=3,third=4, type=2,number=TRUE)
}
# biplot with influence genotype in pc3 vs pc2
if (requireNamespace("spdep", quietly = TRUE)) {
plot(model,first=2,second=3, type=3,number=TRUE,icol="green")
}
# }
Run the code above in your browser using DataLab