Conditional probability distribution (CPD) method in temperature based death time estimation: Error propagation analysis
Introduction
In their recent paper [2] Biermann and Potente applied Bayesian estimation to death time determination by the nomogram method (NM) and the compound method (CM) introducing this application as conditional probability distribution (CPD) method. The authors describe the computation of conditional probabilities and demonstrate the power of the method in cases from a field study by Henssge [3]. This study was originally intended to emphasize the value of the CM in reducing the 95%-confidence interval width from the NM. By integrating external information when the deceased was last seen alive and found dead the CPD-method was able to further narrow the width of the 95%-confidence intervals. Biermann and Potente [2] conclude that “the CPD-method … will always lead to a correct 95.45% interval for use in court and investigations. This interval will always be smaller than the initial interval determined by the NM algorithm and in some cases smaller than the interval after application of CM.”
While CPD is in principle a promising tool for a more accurate death time estimation, Biermann and Potente [2] do not discuss potential error sources of the method. But the CPD-method does not only allow integration of external information that sets reliable limits to the possible death time interval but also computation of probabilities for small time intervals of interest (e.g. no-alibi time intervals of suspects). In the light of the importance of the CPD-method as court evidence we identify and complement potential error sources.
The knowledge of the correct probability distribution of the death time estimator t^ from temperature based death time estimation (TDE) represents the central presupposition of CPD. CPD initially estimates this probability distribution. If this estimated probability distribution is biased – e.g. by a systematic model error in TDE or by random variation of the estimated TDE value from the true value – the results of CPD may be severely distorted. We deduce a formula for the errors in CPD results caused by deviations of the estimator t^ from the true value and provide asymptotic approximations of these errors for large t^-deviations. The study also highlights a paradox of CPD in certain case scenarios, in which an increasing TDE-induced bias produces increasing CPD-probability values (up to 100%) for small time intervals. We demonstrate the results using a synthetic homicide case example.
Section snippets
Death time estimation and inherent errors
TDE tries to reconstruct the time difference between true death time t and time tM of rectal temperature measurement TM. The measured rectal temperature TM = T(tM, ?) monotonously depends on time tM; its (physiological) initial value T0 = T(t, ?) at death time t is known. The parameter ? refers to any vector of measurable quantities, which in case of TDE by Marshall and Hoare with parameter values by Henssge (NM) [4], [5], [6], [7], [8], [9], [10] are: rectal temperature T0 at death time t,
Bayesian probability and external information
We will now describe the abstract framework for CPD-application: Let a < b be two points in time and let further r be a value of the death time estimator t^. The starting point of the Bayesian approach (see e.g. [15], [16]) is to ask for the probability P(t?[a,b]|t^ = r) of the true death time t lying in the interval [a,b] under the condition that the estimator t^ assumes the value r. The symbol P(t|t^ = r) represents the a posteriori distribution or posterior, whereas P(t) stands for the a priori
General influence of TDE-induced bias on Bayesian estimation
The core problem of CPD is, that it requires the correct probability distribution of an (in the mathematical sense) unbiased estimator as input. In case of CPD-application in TDE this very distribution is estimated since CPD takes the estimated death time as the expected value and therefore as the ‘center of mass’ of the probability distribution. Any deviation of the estimator t^s value from the true value t will be transformed to a bias of the CPD-estimated distribution. If the TDE value lies
Synthetic case example
For the sake of clarity we present the results in the form of a hypothetical application example. The hypothetical example (E) was constructed to show the power and the risks of Bayesian estimation and CPD.
The scenario for the example (E) is a homicide investigation. There is non-temperature related information from testimonies that the deceased was still alive at time c = 10:00 and found dead at time d = 16:00 the same day, which attributes a 100% probability a priori to the time of death in the
Discussion
The present study deals with errors caused by estimator deviations in the conditional probability distribution (CPD) method [1], [2]. The general mathematical framework is the well-known Bayesian estimation. Biermann and Potente apply Bayesian estimation in their CPD-method to integrate external information (e.g. from testimonies) in death time estimation and to calculate probabilities for time intervals of special interest (e.g. no-alibi-time intervals of suspects).
All TDE approaches use
References (19)
- et al.
The deployment of conditional probability distributions for death time estimation
Forensic Sci. Int.
(2011) Death time estimation in case work I. The rectal temperature time of death nomogram
Forensic Sci. Int.
(1988)- et al.
Integration qualitativer Daten in die temperaturbasierte Todeszeitbestimmung nach Henßge durch Verwendung bedingter Wahrscheinlichkeiten [Abstract]
Rechtsmedizin
(2008) - et al.
Experiences with a compound method for estimating the time since death II. Integration of non-temperature-based methods
Int. J. Legal Med.
(2000) - et al.
Präzisionsvergleich von Todeszeitrückrechnungen aus der Rektaltemperatur ohne und mit Berücksichtigung von Einflussfaktoren
Beitr. Gerichtl Med.
(1985) - et al.
Estimating the time of death: the rectal cooling after death and its mathematical expression
J. Forensic Sci.
(1962) - et al.
Experiences with a compound method for estimating the time since death I. Rectal temperature nomogram for time since death
Int. J. Legal Med.
(2000) Die Präzision von Todeszeitschätzungen durch die mathematische Beschreibung der rektalen Leichenabkühlung
Z. Rechtsmed
(1979)- et al.
Methoden zur Bestimmung der Todeszeit an Leichen
(1988)
Cited by (7)
Prism – A novel approach to dead body cooling and its parameters
2021, Forensic Science InternationalCitation Excerpt :After calculation of ‘t’, a TsD interval of 2 SD (95.45%) width is produced. It may be limited further using the so called compound method or non-medical limitations in conditional probability distributions [10,11]. MHH is based on single rectal temperature readings only.1
Minimum time since death when the body has either reached or closely approximated equilibrium with ambient temperature
2017, Forensic Science InternationalCitation Excerpt :The NM provides a death time interval of 2 SD (95.45%) width in the form of a Gaussian distribution, which may be altered for certain ranges of cooling (see below) and Ta and for non standard cooling conditions (CF ≠ 1.0). Resulting death time intervals may be further limited due to findings in the so-called compound method [2] or due to non-medical limitations [13,14], taking confidence interval limits into account. The originally published confidence intervals [15,16] have been subject of debate recently [17–20].
A quantitative method for risk assessment of agriculture due to climate change
2018, Theoretical and Applied ClimatologyThermology 2014 - A computer-assisted literature survey
2015, Thermology International