mixpar {DNAmixturesLite} | R Documentation |
Parameters for DNA mixture models
Description
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
mixpar(rho = NULL, eta = NULL, xi = NULL, phi = NULL, parlist = NULL)
## S3 method for class 'mixpar'
print(x, digits = max(3L, getOption("digits") - 3L), scientific = FALSE, ...)
Arguments
rho |
Amplification factor. |
eta |
Scale parameter in the gamma distribution. |
xi |
Stutter parameter. |
phi |
Named vector of the fraction of DNA contributed by each contributor. |
parlist |
A list of parameters of class |
x |
An object of class |
digits |
The number of digits to print |
scientific |
Should scientific notation be used? |
... |
arguments passed to print |
Details
The mixture parameter is a two-way array of lists;
columns correspond to the four model parameters rho
,
eta
, xi
, and phi
, and rows correspond to the
mixtures included in the model.
The print method is currently somewhat specialised, in that it
assumes that rho
, eta
, and xi
are merely real
numbers. phi
is assumed to consist of a named vector per
mixture; the names, or order of names, can differ between mixtures.
Value
An object of class "mixpar".
Author(s)
Therese Graversen
Examples
## A parameter for two mixtures
q <- mixpar(rho = list(30, 30), eta = list(30, 30), xi = list(0.08, 0.08),
phi = list(c(Anna = 0.5, Peter = 0.2, U1 = 0.3),
c(U1 = 0.5, Anna = 0.2, Peter = 0.3)))
## Equivalent to specifying the parameter for each mixture and then combining.
p1 <- mixpar(rho = list(30), eta = list(30), xi = list(0.08),
phi = list(c(Anna = 0.5, Peter = 0.2, U1 = 0.3)))
p2 <- mixpar(rho = list(30), eta = list(30), xi = list(0.08),
phi = list(c(U1 = 0.5, Anna = 0.2, Peter = 0.3)))
p <- mixpar(parlist = list(p1, p2))