####=========================================####
#### For CRAN time limitations most lines in the
#### examples are silenced with one '#' mark,
#### remove them and run the examples
####=========================================####
data(DT_fulldiallel)
DT <- DT_fulldiallel
head(DT)
mix <- mmer(stems~1, random=~female+male, data=DT)
summary(mix)
mixb <- mmec(stems~1, random=~female+male, data=DT)
summary(mixb)$varcomp
####=========================================####
####=========================================####
#### Multivariate model example
####=========================================####
####=========================================####
data(DT_fulldiallel)
DT <- DT_fulldiallel
head(DT)
mix <- mmer(cbind(stems,pods,seeds)~1,
random=~vsr(female) + vsr(male),
rcov=~vsr(units),
data=DT)
summary(mix)
#### genetic variance covariance
cov2cor(mix$sigma$`u:female`)
cov2cor(mix$sigma$`u:male`)
cov2cor(mix$sigma$`u:units`)
Run the code above in your browser using DataLab