Forensic DNA random-match probability and likelihood ratio

Product rule over loci with the NRC II Recommendation 4.10 subpopulation correction (theta): homozygote [2t+(1-t)p][3t+(1-t)p]/((1+t)(1+2t)); heterozygote 2[t+(1-t)p1][t+(1-t)p2]/((1+t)(1+2t)); LR = 1/RMP. theta = 0 recovers the Hardy-Weinberg product rule.

Usage

morie_forsnp(genotype, freqs, theta = 0)

Arguments

genotype

List of length-2 allele vectors per locus.

freqs

List of per-locus named allele-frequency vectors.

theta

Coancestry coefficient.

Value

A list with elements rmp, lr, locus_rmp, n_loci, theta, method.

References

Buckleton, J., Triggs, C. M. and Walsh, S. J. (2005). Forensic DNA Evidence Interpretation. CRC Press.

Examples

genotype <- list(c("A", "a"), c("B", "B"))
freqs <- list(c(A = 0.6, a = 0.4), c(B = 0.7, b = 0.3))
morie_forsnp(genotype, freqs)
#> $rmp
#> [1] 0.2352
#> 
#> $lr
#> [1] 4.251701
#> 
#> $locus_rmp
#> [1] 0.48 0.49
#> 
#> $n_loci
#> [1] 2
#> 
#> $theta
#> [1] 0
#> 
#> $method
#> [1] "forensic RMP/LR, NRC II 4.10 (Buckleton 2005)"
#>