# NOT RUN {
data(CPdata)
head(CPpheno)
CPgeno[1:4,1:4]
#### create the variance-covariance matrix
A <- A.mat(CPgeno)
#### look at the data and fit the model
head(CPpheno)
mix1 <- mmer2(color~1,
random=~ eig(g(id)),
rcov=~units,
G=list(id=A),
data=CPpheno)
summary(mix1)
# mix2 <- mmer2(cbind(Yield,color)~1,
# random=~ us(trait):eig(g(id)),
# rcov=~us(trait):units,
# G=list(id=A),
# data=CPpheno)
# summary(mix2)
# }
Run the code above in your browser using DataLab