dataL: Generic formula of the likelihood of the observed alleles in...

In hindantation/forensim: Statistical Tools for the Interpretation of Forensic DNA Mixtures

Generic formula of the likelihood of the observed alleles in a mixture conditional on the number of contributors for a specific locus

Description

The function dataL gives the likelihood of a set of alleles observed at a specific locus conditional on the number of contributors that gave these alleles. Calculation is based upon the frequencies of the observed alleles.

Usage

dataL(x = 1, p, theta = 0)

Arguments

x	an integer giving the number of contributors
p	a numeric vector giving the frequencies of the observed alleles in the mixture
theta	a float in [0,1[. theta is equivalent to Wright's Fst. In case of population subdivision, it allows a correction of the allele frequencies in the subpopulation of interest

Note

dataL function has several similarities with the Pevid.gen function of the forensic package which computes the probability of the DNA evidence, dataL implements a particular case of this probability. Please see https://CRAN.R-project.org/package=forensim

Author(s)

Hinda Haned [e-mail]

References

Haned H, Pene L, Lobry JR, Dufour AB, Pontier D. Estimating the number of contributors to forensic DNA mixtures: Does maximum likelihood perform better than maximum allele count? J Forensic Sci, accepted 2010.

Curran JM, Triggs CM, Buckleton J, Weir BS. Interpreting DNA Mixtures in Structured Populations. J Forensic Sci 1999;44(5): 987-995

See Also

lik.loc and lik for calculating the likelihood of a given simumix object

Examples

## Not run: 
#likelihood of observing two alleles at frequencies 0.1 and 0.01 when the number of 
#contributors is 2, in two cases:  theta=0 and theta=0.03
dataL(x=2,p=c(0.1,0.01), theta=0)
dataL(x=2,p=c(0.1,0.01), theta=0.03)

## End(Not run)

hindantation/forensim documentation built on Oct. 8, 2022, 4:42 a.m.