R: Set tables for auxiliary variables 'D'

setCPT.D {DNAmixturesLite}

R Documentation

Set tables for auxiliary variables `D`

Description

Function for setting the tables on all nodes D_r_a for mixture r in the network corresponding to one marker, using peak heights and a set of parameters.

IMPORTANT: This is the DNAmixturesLite package, which is intended as a service to enable users to try DNAmixtures without purchasing a commercial licence for Hugin. When at all possible, we strongly recommend the use of DNAmixtures rather than this lite-version. See https://dnamixtures.r-forge.r-project.org/ for details on both packages.

While the lite-version seeks to provide the full functionality of DNAmixtures, note that computations are much less efficient and that there are some differences in available functionality. Be aware that the present documentation is copied from DNAmixtures and thus may not accurately describe the implementation of this lite-version.

Usage

setCPT.D(
  domain,
  rho,
  xi,
  eta,
  phi,
  C,
  n.unknown,
  n_K,
  D,
  gets_stutter,
  can_stutter,
  stutter.from
)

Arguments

`domain`	A `hugin.domain` modelling one marker
`rho`	rho
`xi`	xi
`eta`	eta
`phi`	phi
`C`	Detection threshold
`n.unknown`	Number of unknown contributors
`n_K`	Allelecounts for known contributors
`D`	Names for binary nodes
`gets_stutter`	Does the allele receive stutter?
`can_stutter`	Can the allele stutter?
`stutter.from`	From which allele does the stutter come?

Details

The function sets conditional probabilities for auxiliary variables D_r_a, which are used for obtaining probabilities of a peak falling above resp. below the detection threshold under the current model. The conditional probability tables are defined using the gamma c.d.f. as

P(\texttt{D[a]} = \texttt{TRUE} | n_{ia}, n_{i,a+1}, i = 1,\ldots, k) = G(\texttt{C}).

Value

invisibly NULL