pContrib_locus {DNAtools} | R Documentation |
Compute the posterior probabilities for
for a given prior
.
Description
Compute a matrix of posterior probabilties where
ranges from 1 to
, and
is
. This is done by evaluating
, where
is evaluated by
pNoA
.
Usage
pContrib_locus(
prob = NULL,
m.prior = NULL,
m.max = 8,
pnoa.locus = NULL,
theta = 0
)
Arguments
prob |
Vectors with allele probabilities for the specific locus |
m.prior |
A vector with prior probabilities (summing to 1), where the
length of |
m.max |
Derived from the length of |
pnoa.locus |
A named vector of locus specific probabilities
|
theta |
The coancestery coefficient |
Details
Computes a matrix of values for a specific locus.
Value
Returns a matrix for
and
.
Author(s)
Torben Tvedebrink, James Curran
References
T. Tvedebrink (2014). 'On the exact distribution of the number of alleles in DNA mixtures', International Journal of Legal Medicine; 128(3):427–37. <https://doi.org/10.1007/s00414-013-0951-3>
Examples
## Simulate some allele frequencies:
freqs <- simAlleleFreqs()
## Compute Pr(m|n0) for m = 1, ..., 5 and n0 = 1, ..., 10 for the first locus:
pContrib_locus(prob = freqs[[1]], m.max = 5)