Objects of class
"pedigree" represent a set of
individuals that can have two parents including their parent-child
relations. The terminology has been taken from cattle breeding.
Objects from the Class
Objects can be created by calls of the form
...), or more conveniently,
pedigree(sire= ., dam = ., label = .).
R. A. Mrode, Linear Models for the Prediction of Animal Breeding Values, 2nd ed, CABI Publishing, 2005.
## Rather use, pedigree()! The following is "raw code": ## Example from chapter 2 of Mrode (2005) p1 <- new("pedigree", sire = as.integer(c(NA,NA,1, 1,4,5)), dam = as.integer(c(NA,NA,2,NA,3,2)), label = as.character(1:6)) p1 (dtc <- as(p1, "dtCMatrix")) # T-inverse in Mrode's notation solve(dtc) chol(p1) # (left) Cholesky factor of A, the relationship matrix