logLX: Evaluates mixture log likelihood when extra findings present

In KinMixLite: Inference About Relationships from DNA Mixtures

Evaluates mixture log likelihood when extra findings present

Description

Replacement for logL in DNAmixtures that calls calls LogL.UKX instead of logL.UK.

Usage

logLX(mixture, expr.extra.findings, presence.only = FALSE, initialize = FALSE)

Arguments

mixture	Compiled DNAmixture object.
expr.extra.findings	expression containing the extra findings
presence.only	Set to TRUE to ignore peak heights and evaluate the likelihood function considering peak presence and absence (heights above and below threshold) only. Defaults to FALSE
initialize	should all entered evidence be removed from the networks in mixture

Value

The log likelihood.

Author(s)

Peter Green ([e-mail])

See Also

See also logL.

Examples

data(test2data)

# set threshold C
C<-0.001

mixD<-DNAmixture(list(epg),k=2,C=list(0.001),database=db)
pars<-mixpar(rho=list(2),eta=list(100),xi=list(0.1),phi=list(c(U1=0.7,U2=0.3)))
baseline<-logL(mixD)(pars)

mixMBN<-DNAmixture(list(epg),k=2,C=list(C),database=db,triangulate=FALSE,compile=FALSE)
cgtcaca<-gt2aca(mixMBN,Cgt)
add.child.meiosis.nodes(mixMBN,cgtcaca,1)
log10LR<-(logLX(mixMBN,
	expr.make.findings(list(
	list('Male',ind=1),
	list('Caca',aca='cgtcaca')
	))
	)(pars)-attr(cgtcaca,'logGt')-baseline)/log(10)
cat('log10 LR',log10LR,'\n')

KinMixLite documentation built on Oct. 9, 2024, 9:07 a.m.