Identifies Mendelian inconsistencies between father-child pairs.

check.mendelian.inconsistencies(genotypes, father, child)

Arguments

genotypes

Genotype matrix.

father

Vector of father sample IDs.

child

Vector of child sample IDs.

Value

Data frame summarizing inconsistencies per pair.

Examples

set.seed(1)
genotypes <- matrix(sample(1:3, 30, TRUE), nrow = 3,
                    dimnames = list(c("F1", "C1", "C2"), NULL))
check.mendelian.inconsistencies(genotypes,
                                father = "F1",
                                child  = c("C1", "C2"))
#> F1 - C1 = 2 5 0.4
#> F1 - C2 = 1 5 0.2
#>   Father Child N Total Rate
#> 1     F1    C1 2     5  0.4
#> 2     F1    C2 1     5  0.2