estimatePD: Estimate the drop-out probability based on number of alleles
In DNAtools: Tools for Analysing Forensic Genetic DNA Data

estimatePD

R Documentation

Estimate the drop-out probability based on number of alleles

Description

An inferior may to estimate the drop-out probability compared to using the peak heights from the electropherogram. However, to compare the performance with Gill et al. (2007) this implements a theoretical approach based on their line of arguments.

Usage

estimatePD(n0, m, pnoa = NULL, probs = NULL, theta = 0, locuswise = FALSE)

Arguments

`n0`	Vector of observed allele counts - same length as the number of loci
`m`	The number of contributors
`pnoa`	The vector of Pr(N(m)=n) for n=1,…,2Lm, where L is the number of loci and m is the number of contributors OR
`probs`	List of vectors with allele probabilities for each locus
`theta`	The coancestery coefficient
`locuswise`	Logical. Indicating whether computations should be done locuswise.

Details

Computes the Pr(D) that maximises equation (10) in Tvedebrink (2014).

Value

Returns the MLE of Pr(D) based on equation (10) in Tvedebrink (2014)

Author(s)

Torben Tvedebrink

References

Gill, P., A. Kirkham, and J. Curran (2007). LoComatioN: A software tool for the analysis of low copy number DNA profiles. Forensic Science International 166(2-3): 128 - 138.

T. Tvedebrink (2014). 'On the exact distribution of the number of alleles in DNA mixtures', International Journal of Legal Medicine; 128(3):427–37. <https://doi.org/10.1007/s00414-013-0951-3>

Examples


  ## Simulate some allele frequencies:
  freqs <-  simAlleleFreqs()
  ## Assume 15 alleles are observed in a 2-person DNA mixture with 10 loci:
  estimatePD(n0 = 15, m = 2, probs = freqs)

DNAtools documentation built on March 18, 2022, 7:01 p.m.