makeD: Creates the dominance genetic realationship matrix

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

Given a pedigree, the matrix of coefficients of fraternity are returned - the D matrix. Note, no inbreeding must be assumed. Will return the inverse of the D matrix by default, otherwise this operation can be skipped if desired.

Usage

1
makeD(pedigree, invertD = TRUE)

Arguments

pedigree

A pedigree with columns organized: ID, Dam, Sire

invertD

A logical indicating whether or not to invert the D matrix

Details

Missing parents (e.g., base population) should be denoted by either 'NA' or '0'.

There exists no convenient method of obtaining the inverse of the dominance genetic relatedness matrix (or the D matrix itself) directly from a pedigree (such as for the inverse of A, i.e., Quaas (1995)). Therefore, this function computes the coefficient of fraternity (Lynch and Walsh, 1998) for every individual in the pedigree with a non-zero additive genetic relatedness. Note, the construction of the D matrix is more computationally demanding (in time and space) than is the construction of A.

Value

A

the A matrix in sparse matrix form

D

the D matrix in sparse matrix form

logDet

the log determinant of the D matrix

Dinv

the inverse of the D matrix in sparse matrix form

listDinv

the three column form of the non-zero elements for the inverse of the D matrix

Author(s)

matthewwolak@gmail.com

References

Quaas, R.L. 1995. Fx algorithms. An unpublished note.

Lynch M., & Walsh, B. 1998. Genetics and Analysis of Quantitative Traits. Sinauer, Sunderland, Massachusetts.

See Also

makeDsim

Examples

1
  DinvMat <- makeD(Mrode9)$Dinv


Search within the nadiv package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.