data(grey_cloud_2015)
if (FALSE) {
# takes too long for CRAN
modmat.sire <- model.matrix(~ 0 + fit:paternalID, redata)
modmat.dam <- model.matrix(~ 0 + fit:maternalID, redata)
modmat.siredam <- cbind(modmat.sire, modmat.dam)
rout <- reaster(resp ~ fit + varb,
list(parental = ~ 0 + modmat.siredam, block = ~ 0 + fit:block),
pred, fam, varb, id, root, data = redata)
}
Run the code above in your browser using DataLab