pContrib_locus | R Documentation |
Compute a matrix of posterior probabilties P(m|n0) where m
ranges from 1 to \code{m.max}, and n0 is
0,…,\code{m.max}. This is done by evaluating
Pr(m|n0)=Pr(n0|m)Pr(m)/Pr(n), where
Pr(n0|m) is evaluated by pNoA
.
pContrib_locus( prob = NULL, m.prior = NULL, m.max = 8, pnoa.locus = NULL, theta = 0 )
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 P(N(m)=n), n=1,…,2m. |
theta |
The coancestery coefficient |
Computes a matrix of Pr(m|n0) values for a specific locus.
Returns a matrix [Pr(m|n0)] for m = 1,…,m.max and n_0 = 1,…,2m.max.
Torben Tvedebrink, James Curran
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>
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.