data(HDdata)
head(HDdata)
### GCA matrix for half diallel using male and female columns
#Z1 <- hdm(HDdata[,c(3:4)])
#Z2 <- model.matrix(~as.factor(geno)-1, data=HDdata)# SCA matrix
#y <- HDdata$sugar # response
#ETA <- list(list(Z=Z1), list(Z=Z2)) # Zu component
#modHD <- mmer(y=y, Z=ETA, draw=FALSE, silent=TRUE)
#summary(modHD)
Run the code above in your browser using DataLab