findUndisputed: Undisputed identifications in a DVI problem

Description

This function uses the pairwise LR matrix to find undisputed matches between victims and missing individuals. An identification \(V_i = M_j\) is called undisputed, relative to a threshold T, if the corresponding likelihood ratio \(LR_{i,j} \geq T\) AND \(LR_{i,j}\) is at least T times greater than all other pairwise LRs involving \(V_i\) or \(M_j\).

Usage

findUndisputed(
  dvi,
  pairings = NULL,
  ignoreSex = FALSE,
  threshold = 10000,
  strict = FALSE,
  relax = !strict,
  limit = 0,
  nkeep = NULL,
  numCores = 1,
  verbose = TRUE
)

Value

A list with the following entries:

dviReduced: A reduced version of dvi, where undisputed victims/missing persons are removed, and data from undisputed victims inserted into the reference data.
summary: A data frame summarising the undisputed matches.
LRmatrix: Output from pairwiseLR() applied to the reduced problem.

Arguments

dvi: A dviData object, typically created with dviData().
pairings: A list of possible pairings for each victim. If NULL, all sex-consistent pairings are used.
ignoreSex: A logical.
threshold: A non-negative number. If no pairwise LR exceed this, the iteration stops.
strict: A logical affecting the definition of being undisputed (see Details). Default: FALSE.
relax: Deprecated; use strict = FALSE instead.
limit: A positive number. Only pairwise LR values above this are considered.
nkeep: An integer, or NULL. If given, only the nkeep most likely pairings are kept for each victim.
numCores: An integer; the number of cores used in parallelisation. Default: 1.
verbose: A logical. Default: TRUE.

Details

If the parameter strict is set to TRUE, the last criterion is replaced with the stronger requirement that all other pairwise LRs involving \(V_i\) or \(M_j\) must be at most 1.

Examples

Run this code


# \donttest{
u1 = findUndisputed(planecrash, verbose = FALSE)
u1$summary 

# With `strict = TRUE`, the match M3 = V2 goes away
u2 = findUndisputed(planecrash, strict = TRUE, verbose = FALSE)
u2$summary

# Reason: M3 has LR > 1 also against V7
u2$LRmatrix[, "M3"] |> round(2)
# }