sommer (version 1.1)

hdm: Half Diallel Matrix

Description

This function uses a dataframe with 2 columns named "female" and "male" in numeric format and creates an incidence matrix for a single explanatory variable corresponding to the GCA effect. The resulting incidence matrix can be used in the mmer function as a 'Z' argument in the 2-level list argument for random effects.

Usage

hdm(data)

Arguments

data

a dataframe with 2 columns named 'female' and 'male' with numeric or factor values indicating the male or female used to produce such hybrid.

Value

• [object Object]

References

Fikret Isik. 2009. Analysis of Diallel Mating Designs. North Carolina State University, Raleigh, USA.

Examples

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)