# NOT RUN {
data(australia.soybean)
dat <- australia.soybean
require(reshape2)
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)], main="australia.soybean", horiz=FALSE)
# }
Run the code above in your browser using DataLab