Description Usage Arguments Details Value Author(s) References See Also Examples
Calculates conditional probability of DNA evidence, given proposition about who the known and unknown contributors to the mixture were. All the individuals involved in a case are assumed to come from the same subpopulation with a given coancestry coefficient. Independence of alleles in the subpopulation (i.e., relatives are excluded) and dependence in the whole population is assumed.
1 |
alleles |
vector of distinct alleles (from one specific locus) found in the crime sample. |
prob |
vector of corresponding allele proportions in a population. The allele proportions are estimated from the whole population. |
x |
nonnegative integer. The number of unknown contributors to the mixture. |
T |
object of class genotype (package genetics), or a vector of strings where each string contains two alleles seperated by '/', corresponding to one known contributor. The length of the vector equals the number of known contributors. Default is NULL. |
V |
object of class genotype (package genetics), or a vector of strings where each string contains two alleles seperated by '/', corresponding to one known noncontributor. The length of the vector equals the number of known noncontributors. Default is NULL. |
theta |
number from the interval [0,1). Coancestry coefficient |
The general formula for the evaluation of the probability of DNA evidence (and its derivation) can be found in Zoubkova and Zvarova (2004).
If theta = 0
, Pevid.gen
returns the same numerical result as Pevid.ind
.
probability of the DNA evidence
Miriam Marusiakova maruskay@gmail.com
The work was supported by the project 1M06014 of the Ministry of Education, Youth and Sports of the Czech Republic.
Curran JM, Triggs CM, Buckleton J, Weir BS (1999), Interpreting DNA mixtures in structured populations, Journal of Forensic Sciences 44, 987-995.
Evett IW, Weir BS (1998), Interpreting DNA evidence; Statistical genetics for forensic scientists. Sinauer, Sunderland, MA
Fung WK, Hu YQ (2000), Interpreting forensic DNA mixtures: allowing for uncertainty in population substructure and dependence, Journal of the Royal Statistical Society A 163, 241-254.
National Research Council (1996), The evaluation of forensic DNA evidence National Academy Press, Washington, DC.
Zoubkova K, Zvarova J (2004), Statisticke metody ve forenzni genetice, Master's thesis, Charles University, Prague.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | ## Rape case
## The evidence profile:
m <- c(13, 14, 15)
## the victim's genotype:
victim <- "13/14"
## the suspect's genotype
suspect <- "15/15"
## frequencies of alleles {13, 14, 15}:
p <- c(0.042, 0.166, 0.106)
## consensual partner of the victim
partner <- "15/16"
## Prosecution proposition:
## Contributors were the victim and the suspect.
Pevid.gen( alleles = m, prob = p, T = c(victim, suspect),
V = partner, x = 0, theta = 0.03)
## Defence proposition:
## Contributors were the victim and one unknown person.
##
## Likelihood ratio for DNA evidence:
## structured population
1/Pevid.gen( alleles = m, prob = p, T = victim,
V = c(suspect, partner), x = 1, theta = 0.03)
## Note: a person carrying both alleles different from the alleles
## in the crime sample (e.g., with genotype "16/16") has no effect
## on the value of LR:
1/Pevid.gen( alleles = m, prob = p, T = victim,
V = c(suspect, partner, "16/16"), x = 1, theta = 0.03)
## But the consensual partner of the victim having genotype "15/16"
## influences the value of LR, compare:
1/Pevid.gen( alleles = m, prob = p, T = victim, V = suspect,
x = 1, theta = 0.03)
##
## population in Hardy - Weinberg equilibrium
1/Pevid.gen( alleles = m, prob = p, T = victim,
V = c(suspect, partner), x = 1)
1/Pevid.gen( alleles = m, prob = p, T = victim, x = 1)
## compare
1/Pevid.ind( alleles = m, prob = p, u = 15, x = 1)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.