DT_fulldiallel | R Documentation |
This dataset contains phenotpic data for 36 winter bean hybrids, coming from a full diallel design and evaluated for 9 traits. The column male and female origin columns are included as well.
data("DT_fulldiallel")
The format is: chr "DT_fulldiallel"
This data was generated by a winter bean study and originally included in the agridat package.
Covarrubias-Pazaran G (2016) Genome assisted prediction of quantitative traits using the R package sommer. PLoS ONE 11(6): doi:10.1371/journal.pone.0156744
The core functions of the package mmer
and mmec
####=========================================####
#### 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`)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.