likeLTD::locus.likes.peaks

Description

Vector with individual locus likelihoods for peak height data.

Usage

locus.likes.peaks(hypothesis,results,...)

Arguments

Details

Convert the overall likelihood returned by evaluate.peaks into locus specific likelihoods.

Value

Vector.

See Also

evaluate.peaks

Examples

## Not run: 
# datapath to example files
datapath = file.path(system.file("extdata", package="likeLTD"),"laboratory")

# File paths and case name for allele report
admin = pack.admin.input.peaks(
            peaksFile = file.path(datapath, 'laboratory-CSP.csv'),
            refFile = file.path(datapath, 'laboratory-reference.csv'),
            caseName = "Laboratory",
            detectionThresh = 20
             )

# Enter arguments
args = list(
        nUnknowns = 1
        )

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

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

# reduce number of iterations for demonstration purposes
paramsP$control$itermax=25
paramsD$control$itermax=25

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

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

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

## End(Not run)

Results