IS {relSim}

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)