####=========================================####
#### For CRAN time limitations most lines in the
#### examples are silenced with one '#' mark,
#### remove them and run the examples using
#### command + shift + C |OR| control + shift + C
####=========================================####
data(DT_rice)
# DT <- DT_rice
# GT <- GT_rice
# GTn <- GTn_rice
# head(DT)
# M <- atcg1234(GT)
# A <- A.mat(M$M)
# mix <- mmer(Protein.content~1,
# random = ~vsr(geno, Gu=A) + geno,
# rcov=~units,
# data=DT)
# summary(mix)$varcomp
#
# Ai <- as(solve(A + diag(1e-6,ncol(A),ncol(A))), Class="dgCMatrix")
# mixb <- mmec(Protein.content~1,
# random = ~vsc(isc(geno), Gu=Ai) + geno,
# rcov=~units,
# data=DT)
# summary(mixb)$varcomp
Run the code above in your browser using DataLab