dat <- australia.soybean
dm <- melt(dat, id.var=c('env', 'year','loc','gen'))
# Joint plot of genotypes & traits. Similar to Figure 1 of Kroonenberg 1989
dmat <- acast(dm, gen~variable, fun=mean)
dmat <- scale(dmat)
biplot(princomp(dmat), main="australia.soybean trait x gen biplot", cex=.75)
# Figure 1 of Kozak 2010, lines 44-58
dmat2 <- dmat[44:58,]
require("lattice")
parallelplot(dmat2[,c(2:6,1)], horiz=FALSE)Run the code above in your browser using DataLab