IS: Use importance sampling to determine the probability of m...

In relSim: Relative Simulator

Use importance sampling to determine the probability of m peaks from n contributors to a mixture

Description

WARNING: This function is experimental.

Usage

IS(
  Freqs,
  numContributors = 4,
  maxPeaks = NULL,
  numIterations = 100,
  bTail = FALSE
)

Arguments

Freqs	a set of allele frequencies. The format can be found in readFreqs
numContributors	the number of contributors to each mixture. Must be >= 1.
maxPeaks	either the number of peaks observed in the mixture or such that 1 <=1 maxPeaks <= min(2 * numContribuors, numAlleles). That is, if used in this way maxPeaks must be between 1 and the smaller of twice the number of contributors or the number of possible alleles because you cannot see more peaks than there are possible alleles.
numIterations	the number of iterations to use in the importance sampling scheme.
bTail	if TRUE then the tail probability is calculated.

Value

a list with as many elements as loci. If tail probabilities are selected then each locus element will be a vector of probabilities

Examples

## Not run: 
data(USCaucs)
IS(USCaucs, numContributors = 4, maxPeaks = 3, numIterations = 1e4)

## End(Not run)

relSim documentation built on Aug. 29, 2023, 9:07 a.m.