## S3 method for class '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,...)
AMMI
library(agricolae)
data(plrv)
attach(plrv)
model<- AMMI(Locality, Genotype, Rep, Yield)
detach(plrv)
# biplot
# plot PC2 vs PC1
plot(model)
# plot PC1 vs Yield
plot(model,0,1,gcol="blue",ecol="green")
# triplot PC 2,3,4
require(klaR) # type=2
plot(model,first=2,second=3,third=4, type=2,number=TRUE)
# biplot with influence genotype in pc3 vs pc2
require(spdep) # type=3
plot(model,first=2,second=3, type=3,number=TRUE,icol="green")
Run the code above in your browser using DataLab