likeLTD (version 6.3.0)

locus.likes: likeLTD::locus.likes

Description

Vector with individual locus likelihoods.

Usage

locus.likes(hypothesis,results,…)

Arguments

hypothesis

The hypothesis generated by either prosecution.hypothesis or defence.hypothesis.

results

Either prosecution or defence results from evaluate e.g. results$Pros or results$Def.

…

Any extra parameter that was passed to evaluate or optimisation.params.

Value

Vector.

Details

Convert the overall likelihood returned by DEoptim into locus specific likelihoods.

See Also

evaluate

Examples

# NOT RUN {
# datapath to example files
datapath = file.path(system.file("extdata", package="likeLTD"),"hammer")

# File paths and case name for allele report
admin = pack.admin.input(
            cspFile = file.path(datapath, 'hammer-CSP.csv'),
            refFile = file.path(datapath, 'hammer-reference.csv'),
            caseName = "hammer",
	    kit= "SGMplus"
             )

# Enter arguments
args = list(
        nUnknowns = 1,
        doDropin = FALSE,
        ethnic = "EA1",
        adj = 1,
        fst = 0.02,
        relatedness = c(0,0)
        )

# Create hypotheses
hypP = do.call(prosecution.hypothesis, append(admin,args))
hypD = do.call(defence.hypothesis, append(admin,args))

# Get parameters for optimisation
paramsP = optimisation.params(hypP)
paramsD = optimisation.params(hypD)

# Run optimisation
# n.steps set for demonstration purposes
results = evaluate(paramsP, paramsD, n.steps=1)

# get locus likelihoods under prosecution
locus.likes(hypP,results$Pros)

# get locus LRs
locus.likes(hypP,results$Pros)/locus.likes.peaks(hypD,results$Def)
# }