mixMLX {KinMixLite}

Maximises mixture likelihood when extra findings present

Description

Replacement for mixML in DNAmixtures that calls logLX instead of logL.

Usage

mixMLX(mixture, expr.extra.findings, pars, constraints = NULL, phi.eq = FALSE, 
    val = NULL, trace = FALSE, order.unknowns = TRUE, initialize = FALSE, ...)

Arguments

mixture

Compiled DNAmixture object.

expr.extra.findings

expression containing the extra findings

pars

Parameters, in mixpar format.

constraints

as in mixML

phi.eq

as in mixML

val

as in mixML

trace

as in mixML

order.unknowns

as in mixML

initialize

should all entered evidence be removed from the networks in mixture

...

as in mixML

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 ([e-mail])

See Also

See also mixML.

Examples

data(test2data)

# set threshold C
C<-0.001

mixD<-DNAmixture(list(epg),k=2,C=list(C),database=db)

# find MLE's and maximised likelihood
# adding evidence individual 1 is Male 

expr.extra.findings<-expr.make.findings(list(list('Male',ind=1)))

startpar<-mixpar(rho=list(60),eta=list(24),xi=list(0.16),phi=list(c(U1=0.75,U2=0.25)))
mlDM<-mixMLX(mixD,expr.extra.findings,startpar,trace=FALSE) 
pars<-mlDM$mle
cat('\nBaseline model maximised likelihood:',mlDM$lik,'\n')
cat('and MLEs:\n')
print(mlDM$mle)