Conditional probability distribution (CPD) method in temperature based death time estimation: Error propagation analysis

https://doi.org/10.1016/j.forsciint.2014.02.016Get rights and content

Abstract

Bayesian estimation applied to temperature based death time estimation was recently introduced as conditional probability distribution or CPD-method by Biermann and Potente. The CPD-method is useful, if there is external information that sets the boundaries of the true death time interval (victim last seen alive and found dead). CPD allows computation of probabilities for small time intervals of interest (e.g. no-alibi intervals of suspects) within the large true death time interval. In the light of the importance of the CPD for conviction or acquittal of suspects the present study identifies a potential error source. Deviations in death time estimates will cause errors in the CPD-computed probabilities. We derive formulae to quantify the CPD error as a function of input error. Moreover we observed the paradox, that in cases, in which the small no-alibi time interval is located at the boundary of the true death time interval, adjacent to the erroneous death time estimate, CPD-computed probabilities for that small no-alibi interval will increase with increasing input deviation, else the CPD-computed probabilities will decrease. We therefore advise not to use CPD if there is an indication of an error or a contra-empirical deviation in the death time estimates, that is especially, if the death time estimates fall out of the true death time interval, even if the 95%-confidence intervals of the estimate still overlap the true death time interval.

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)

  • F.M. Biermann et al.

    The deployment of conditional probability distributions for death time estimation

    Forensic Sci. Int.

    (2011)
  • C. Henssge

    Death time estimation in case work I. The rectal temperature time of death nomogram

    Forensic Sci. Int.

    (1988)
  • S. Potente et al.

    Integration qualitativer Daten in die temperaturbasierte Todeszeitbestimmung nach Henßge durch Verwendung bedingter Wahrscheinlichkeiten [Abstract]

    Rechtsmedizin

    (2008)
  • C. Henssge 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)
  • E. Stipanits et al.

    Präzisionsvergleich von Todeszeitrückrechnungen aus der Rektaltemperatur ohne und mit Berücksichtigung von Einflussfaktoren

    Beitr. Gerichtl Med.

    (1985)
  • T.K. Marshall et al.

    Estimating the time of death: the rectal cooling after death and its mathematical expression

    J. Forensic Sci.

    (1962)
  • C. Henssge 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)
  • C. Henssge

    Die Präzision von Todeszeitschätzungen durch die mathematische Beschreibung der rektalen Leichenabkühlung

    Z. Rechtsmed

    (1979)
  • C. Henßge et al.

    Methoden zur Bestimmung der Todeszeit an Leichen

    (1988)
There are more references available in the full text version of this article.

Cited by (7)

  • Prism – A novel approach to dead body cooling and its parameters

    2021, Forensic Science International
    Citation 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 International
    Citation 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].

View all citing articles on Scopus
View full text