if (FALSE) {
library(agridat)
data(bond.diallel)
dat <- bond.diallel
# Because these data are means, we will not be able to reproduce
# the anova table in Bond. More useful as a multivariate example.
libs(corrgram)
corrgram(dat[ , 3:11], main="bond.diallel",
lower=panel.pts)
# Multivariate example from sommer package
corrgram(dat[,c("stems","pods","seeds")],
lower=panel.pts, upper=panel.conf, main="bond.diallel")
libs(sommer)
m1 <- mmer(cbind(stems,pods,seeds) ~ 1,
random= ~ vs(female)+vs(male),
rcov= ~ vs(units),
dat)
#### genetic variance covariance
cov2cor(m1$sigma$`u:female`)
cov2cor(m1$sigma$`u:male`)
cov2cor(m1$sigma$`u:units`)
}
Run the code above in your browser using DataLab