calcLRmcmc: calcLRmcmc

In oyvble/euroformix: A Quantitative Model for Interpreting DNA Mixtures

calcLRmcmc

Description

A function for calculating LR from MCMC simulations

Usage

calcLRmcmc(
  mlefitHp,
  mlefitHd,
  niter = 2000,
  delta = 2,
  quantile = 0.05,
  seed = NULL,
  accRate = 0.25,
  accRateTol = 0.1,
  niterTune = 200,
  diffSeed = 999,
  verbose = TRUE,
  traceplot = TRUE,
  mcmcObjList = NULL
)

Arguments

mlefitHp	A fitted object returned from contLikMLE (under Hp)
mlefitHd	A fitted object returned from contLikMLE (under Hd)
niter	Number of samples in the MCMC-sampling.
delta	A numerical parameter to scale with the covariance function Sigma. Default is 2. Will be tuned automtically.
quantile	The quantile used to report conservative LR
seed	The user can set seed if wanted
accRate	aimed acceptance rate
accRateTol	Difference tolerance regarding acceptance rate
niterTune	Number of samples for tuning delta in the MCMC-sampling.
diffSeed	Seed difference
verbose	Whether progress should be printed
traceplot	Whether to add a traceplot with 95% CI
mcmcObjList	A list with hp/hd earlier run with contLikMCMC output

Details

Returns a conservative (quantile) LR or a Full Bayesian based LR

Value

A list with different LR statistics from the MCMC based simulations

Author(s)

Oyvind Bleka

oyvble/euroformix documentation built on Aug. 25, 2023, 11:14 a.m.