forensim (version 2.0)

PE: The random man exclusion probability

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.

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)