# Fitting incomplete block, single tree plot and animal model
dbh.4 <- asreml(dbh ~ 1, random = ~ Rep + Rep/Block+ ped(Tree),
                data = dat, ginverse = list(Tree = pedinv))