| relfactor {pedigreemm} | R Documentation |
Relationship factor from a pedigree
relfactor(ped, labs)
ped |
a pedigree that includes the individuals who occur in labs |
labs |
a character vector or a factor giving the labels to
which to restrict the relationship matrix. If labs is a
factor then the levels of the factor are used as the labels.
Default is the complete set of labels in the pedigree. |
Determine the right Cholesky factor of the relationship matrix for the
pedigree ped, possibly restricted to the specific labels that
occur in labs.
an upper triangular, sparse (right) Cholesky factor of the relationship matrix
## Example from chapter 2 of Mrode (2005)
ped <- pedigree(sire = c(NA,NA,1, 1,4,5),
dam = c(NA,NA,2,NA,3,2), label= 1:6)
(fac <- relfactor(ped))
crossprod(fac) # the relationship matrix