tabAinv: Inverse of the relationship matrix *A* in a tabular format

View source: R/tabAinv.R

tabAinvR Documentation

Inverse of the relationship matrix A in a tabular format

Description

Creates the data.frame of the inverse of the pedigree-based genetic relationship matrix.

Usage

tabAinv(ped, inbr)

Arguments

ped

: data.frame with integer columns corresponding to ID, SIRE, DAM. Missing value is 0.

inbr

: Vector of inbreeding coefficients in the order of individuals in the relationship matrix.

Value

data.frame of the inverse of the genetic relationship matrix

Examples

ped = data.frame(ID=1:6, SIRE=c(0,0,1,3,1,4), DAM=c(0,0,2,2,2,5))
inbr = c(0, 0, 0, 0.25, 0, 0.25)
# or
(inbr = diag(buildA(ped)) - 1)
# or
inbr = tabA(ped); (inbr = inbr[inbr[,1]==inbr[,2],]$a - 1)
# or
# For individual inbreeding values, use function inb.
tabAinv(ped, inbr)


ggroups documentation built on March 28, 2022, 1:06 a.m.