PE: The random man exclusion probability

In forensim: Statistical tools for the interpretation of forensic DNA mixtures

Description

Computes the random man exclusion probability of a mixture stored in a simumix object

Usage

PE(mix, freq, refpop = NULL, theta = 0, byloc = FALSE)

Arguments

mix	a simumix object
freq	a tabfreq object giving the allele frequencies from which to compute the exclusion probability
refpop	character giving the reference population, used only if freq contains allele frequencies for multiple populations
theta	a float from [0,1[ giving Wright's Fst coefficient. theta accounts for population subdivision while computing the likelihood of the data.
byloc	logical, if TRUE, than the exclusion probability is computed per locus, if FALSE (default), the calculations are done for all loci simultaneously

Details

PE gives the exclusion probability at a locus, or at several loci when conditions for Hardy Weinberg are met. If this condition is not met in the population, than a value for theta must be supplied to take into account dependencies between alleles. The formula of the exclusion probability that allows taking into account departure from Hardy Weinberg proportions due to population subdivision was provided by Bruce Weir, please see the references section.

Author(s)

Hinda Haned <[e-mail]>

References

Clayton T, Buckleton JS. Mixtures. In: Buckleton JS, Triggs CM, Walsh SJ, editors. Forensic DNA Interpretation. CRC Press 2005;217-74

Examples

data(strusa)
geno1<-simugeno(strusa,n=c(0,0,100))
mix2 <-simumix(geno1,ncontri=c(0,0,2))
PE(mix2,strusa,"Hisp",byloc=TRUE)

forensim documentation built on May 2, 2019, 6:09 p.m.