DNAmixturesLite (version 0.0-1)

mixpar: 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 mixpar
print(x, digits = max(3L, getOption("digits") - 3L), scientific = FALSE, ...)

Value

An object of class "mixpar".

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 mixpar

x

An object of class mixpar.

digits

The number of digits to print

scientific

Should scientific notation be used?

...

arguments passed to print

Author

Therese Graversen

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.

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))