findNonidentifiable {dvir}

Nonidentifiable missing persons

Description

A missing person in a DVI case is nonidentifiable if unrelated to all (genotyped) reference individuals and all other missing persons in the reference family. It is often wise to ignore such individuals in jointDVI() and other analyses, to relieve the computational burden.

Usage

findNonidentifiable(dvi)

Arguments

dvi

A dviData object, typically created with dviData().

Details

The implementation uses ribd::kinship() to identify individuals having kinship coefficient 0 with all relevant individuals.

Value

A list with the following entries:

• nonidentifiable: A character vector (possibly empty) with the names of the nonidentifiable missing persons.

• dviReduced: A reduced dviData object, where the nonidentifiable individuals are removed from the list of missing persons. If there are no nonidentifiable, this is just a copy of dvi.

• summary: A data frame summarising the findings.

Examples

# Example 1: No nonidentifiables in dataset `example1`
findNonidentifiable(example1)

# Example 2: Add nonidentifiable person "A"
amNew = example1$am[[1]] |>
  addSon(parents = c("NN", "A"))
missNew = c(example1$missing, "A")

dvi = dviData(pm = example1$pm, am = amNew, missing = missNew)
plotDVI(dvi, textAbove = c(A = "nonidentif."))

findNonidentifiable(dvi)