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