protected.mixML {KinMixLite}R Documentation

Protect against numerical errors in maximum likelihood computation

Description

Attempts to catch numerical errors in maximum likelihood computation, by replacing logL values by a default value instead of NaN or other failures

Usage

protected.mixML(mixture, pars, constraints = NULL, phi.eq = FALSE, val = NULL,
  trace = FALSE, order.unknowns = TRUE, default=-999999, ...)

Arguments

mixture

A DNAmixture object.

pars

A mixpar parameter used as a starting value for the optimisation.

constraints

Equality constraint function as function of an array of parameters.

phi.eq

Should the mixture proportions be the same for all traces? Defaults to FALSE.

val

Vector of values to be satisfied for the equality constraints.

trace

Print the evaluations of the likelihood-function during optimisation?

order.unknowns

Should unknown contributors be ordered according to decreasing contributions? Defaults to TRUE.

...

Further arguments to be passed on to solnp.

default

value of logL to be used if numerical errors are encountered

Value

A list containing

mle

The (suggested) MLE.

lik

The log of the likelihood (log e).

as well as the output from the optimisation.

Author(s)

Peter Green (P.J.Green@bristol.ac.uk)


[Package KinMixLite version 2.1.0 Index]