Preservation likelihood ratio from measured evidence

Convenience wrapper: evaluate measured non-invasive evidence under a natural preservation model (\(H_1\)) and an alternative model (\(H_2\)) and return the forensic likelihood ratio. Each model is a list with numeric mean and sd (recycled over the evidence vector).

Usage

morie_taphonomy_preservation_lr(evidence, natural, alternative)

Arguments

evidence

Numeric vector of measured values (e.g. pXRF calcium).

natural

List list(mean=, sd=) for the natural-preservation hypothesis \(H_1\).

alternative

List list(mean=, sd=) for the alternative hypothesis \(H_2\).

Value

The list from morie_taphonomy_likelihood_ratio, with loglik_h1 and loglik_h2 attached.

Examples

morie_taphonomy_preservation_lr(
  evidence = c(1200, 1310, 1180),
  natural = list(mean = 1250, sd = 90),      # lime-processed signature
  alternative = list(mean = 300, sd = 120)   # untreated/decayed signature
)$verbal
#> [1] "extremely strong support for H1"