Q
Function for setting the tables on all nodes Q_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.
setCPT.Q(
domain,
rho,
xi,
eta,
phi,
heights,
n.unknown,
n_K,
Q,
gets_stutter,
can_stutter,
stutter.from
)
Scaling factors to be set as evidence when conditioning on the peak heights
A hugin.domain
modelling one marker
rho
xi
eta
phi ordered as phi[c(U,K)]
where U
and K
are the names of unknown and known contributors.
Observed peak heights
Number of unknown contributors
Allelecounts for known contributors
Names for binary nodes
Does the allele receive stutter?
Can the allele stutter?
From which allele does the stutter come?
Therese Graversen
The function sets conditional probabilities for auxiliary variables Q_r_a
,
which are used for finding probabilities of observing a smaller or
larger peak than the one observed for the DNA mixture.
The conditional probability tables are defined using the gamma c.d.f. as
$$P(\texttt{Q[a]} = \texttt{TRUE} | n_{ia}, n_{i,a+1}, i = 1,\ldots, k) = G(\texttt{heights[a]}).
$$
For further details see setCPT.O
.