library(asreml)
ginvA <- asreml.Ainverse(warcolak)$ginv
ginvD <- makeD(warcolak[,1:3])$listDinv
warcolak$IDD <- warcolak$ID
warcolak.mod <- asreml(trait1 ~ 1, random = ~ped(ID) + giv(IDD),
ginverse = list(ID = ginvA, IDD = ginvD), data = warcolak)
summary(warcolak.mod)$varcomp
varTrans(warcolak.mod)
Run the code above in your browser using DataLab