Gaussian log-likelihood of forensic evidence under a model

Sum of independent normal log-densities for a vector of measured evidence (e.g. pXRF residual calcium, CT tissue density) under a model with expected mean and sd. Building block for morie_taphonomy_likelihood_ratio: evaluate the evidence under each competing hypothesis' model, then take the ratio.

Usage

morie_taphonomy_evidence_loglik(evidence, mean, sd)

Arguments

evidence

Numeric vector of measured values.

mean

Numeric model mean(s), recycled to length(evidence).

sd

Numeric model standard deviation(s) (> 0), recycled.

Value

A single double: the total log-likelihood.

Examples

morie_taphonomy_evidence_loglik(c(1200, 1310), mean = 1250, sd = 80)
#> [1] -11.07849