 Technical Advance
 Open Access
 Published:
Computationally approximated solution for the equation for Henssge’s time of death estimation
BMC Medical Informatics and Decision Making volume 19, Article number: 201 (2019)
Abstract
Background
Time of death estimation in humans for the benefit of forensic medicine has been successfully approached by Henssge, who modelled body cooling based on measurements of Marshall and Hoare. Thereby, body and ambient temperatures are measured at the death scene to estimate a time of death based on a number of assumptions, such as initial body temperature and stable ambient temperature. While so far, practical use of the method resorted to paper print outs or copies of a nomogram using a ruler, increasingly, users are interested in computer or mobile device applications. We developed a computational solution that has been available online as a web accessible PHP program since 2005. From that, we have received numerous requests not so much to detail our code but to explain how to efficiently approximate the solution to the Henssge equation.
Methods
To solve Henssge’s double exponential equation that models physical cooling of a body, it is sufficient to determine a difference term of the equation that will be close to zero for the correct time of death using a discrete set of all sensible possible solutions given that the modelled time frame has practical upper limits. Best postmortem interval approximation yields minimal difference between equation terms
Results
The solution is approximated by solving the equation term difference for a discrete set of all possible time of death intervals that are sensibly found, and by then determining the particular time of death where equation term difference is minimal.
Conclusions
The advantage of a computational model over the nomogram is that the user is also able to model hypothermia and hyperthermia. While mathematically impossible to solve in a straightforward way, solutions to the Henssge equation can be approximated computationally.
Background
The Henssge method for the estimation of time of death [1] a is popular tool for an early postmortem period up to around two days.
It is predominantly used in forensic medicine, where a time of death estimate may have reconstructive and legal relevance [2, 3]. Thereby, body and ambient temperatures are measured at the death scene to estimate a time of death based on a number of assumptions, such as assuming an initial body temperature and assuming a stable ambient temperature. While the computational aspect can be addressed, it is important to point out that the method itself can only estimate a postmortem interval, but it is not a determination or calculation of “the” postmortem interval. Much rather, it is an estimate based on the input parameters that will be approximated by calculation.
The equation that defines this method was simplified by Henssge using the equation developed by Marshall and Hoare (1962) [4]. It mathematically models the temperature curves of actual measured physical cooling experiments. As that, it contains the relevant term that is to be determined by investigating authorities, the time of death estimate t_{i}, twice, each time as part of an exponent.
To allow any user to practically obtain valid and fast solutions, nomograms were created by Henssge for easy replication through photocopies, and published [1]. These contain the relevant parameters in a form where the solution can be rapidly found by using a pencil and a ruler; however, one restriction that these nomograms contain is a fixed value for an assumed initial body temperature of 37,2^{∘}C.
With the advent of increasingly more affordable computation power, calculating the solution of the Henssge equation seems to increasingly become an reality. Users of the conventional nomograms now increasingly require the same method to be made available on mobile devices.
Mathematically, there is no plain and easy solution for this double exponential equation. Since one cannot provide an analytical mathematical solution, approximating the result is the method of choice. We had programmed a fast approximation to practically correct, internally validated solutions using PHP (Hypertext PreProcessor) in 2005 ^{Footnote 1}. Internal validation was performed by crossreferencing arbitrary artificial numerical data across the paper nomograms and by performing stepwise equation term evaluation. After it had been extensively validated by our team then, it was externally validated in 2013 by a South African team [5].
Over the years, we have received numerous requests not so much to detail our code but to explain how to efficiently approximate the solution to the Henssge equation. This will allow users of various programming environments to efficiently program software for this there.
Methods
The method chosen for this solution was that of approximation [6] rather than choosing more complex approaches [7] such as a stepwise iteration originally suggested [1].
For that, the problem’s ideal continuous solution space t is replaced with a finite discrete subset that contains discrete values across a practically sensible subspace A of all possible values of time, t_{A} [6, 8].
Results
The equation to calculate the time of death estimate or interval t=t_{m}−t_{0} for any given point in time t_{m} (time of examination, time of measurement of actual rectal temperature of the body) after time of death t_{0} bases on ambient temperature T_{U}=T_{U}(t_{m}) that is implicitly assumed to be the same as the initial ambient temperature (so T_{U}=T_{U}(t_{m})=T_{U}(t_{0}) is set), an assumed initial body temperature at the time of death T_{D}=T_{body}(t_{0}) (which is usually but not necessarily 37,2^{∘}C), measured actual body core temperature T_{B}=T_{body}(t_{m}) obtained from the environment adjacent to the body, outside of any present cover or clothing, an estimated or actual body weight M and an empirically determined corrective factor f that adapts the calculation to both body weight [9] and ambient conditions (basically modelling environmental aspects of cooling such as humidity and body covers) [10, 11].
For the correct solution for the time of death interval estimate t where factor Z is given in Eq. 3, the following term (up to an ambient temperature T_{U}≤23.2^{∘}C: Eq. 1; for any ambient temperature T_{U}≥23.3^{∘}C: Eq. 2) models the physical cooling according to Henssge [12, 13]:
As the variable t cannot be isolated from these equations, they cannot be analytically solved for t.
The result can be approximated, however. Simply performing a forward calculation of all terms of that equation using a large vector containing a discretized set of N possible time of death estimates t_{A}(i) with i={1,N} will yield N different equations.
The difference between both terms will contain the best approximation for t_{A}(i) once that difference is minimal.
So, a vector t_{A} of discrete times of death n={1,N} is used to determine all N values for the difference vector D with Eq. 4 (T_{U}≤23.2^{∘}C) or Eq. 5 (T_{U}≥23.3^{∘}C):
For the presented best approximation, the one value t_{A}(i) for which D(i)=min(D) is identified.
Discretization and computational error
Practically, N≈300 discrete values are enough to cover the extent of the method between time of death estimates ranging from around an hour to approximately two days. The upper limit of around two days is a technical limitation in that body cooling then has mostly yielded such small differences between ambient and body temperature that estimation of time of death based on cooling or temperature becomes practically irrelevant.
A relatively small error will systematically result from the grain of the values: using decimal time values for t_{A}, increasing t_{A}(i) by tenths of an hour (6 min) may cause time of death estimate errors within maximally ∼1.5% of the usual 95% confidence interval the method reportedly gives (the smallest 95% confidence interval is 2,8 h [12, 14, 15]).
This computational error seems to be of no practical relevance. In the report the forensic pathologist will issue, the notation of results from time of death estimates are usually not given in decimal notation with a 1/10th h step. Issuing results with a maximal precision of a quarter of an hour, but normally with a maximal precision of 1/2 h, will even out the much smaller systematic result deviations incurred by discretization, rendering these effectively insignificant in relation to the application. Any systematic approximation derived error had neither been reported by users since the first installation of our webbased solution in 2005^{Footnote 2}, nor by the study validating our electronic solution in 2013 [5].
Algorithm complexity
With regard to the discretization to a number of N≈300 discrete values, the algorithm will scale with Θ(N). Once defined, this algorithm conforms to a simple complexity of Θ(1). With that, it is laid out to require the same amount of time regardless of the value of the user’s input variables.
Discussion
A particular advantage of a computational model over the nomogram is that the user is also able to model hypothermia and hyperthermia with practical ease. That is possible because the mathematical model of Henssge bases on Marshall and Hoare with a variable initial body temperature [4].
Otherwise, there is no practically relevant difference in terms of the results than when a correctly formatted [5] paper nomogram is used.
Errors of the method and its application
At this point, it may be relevant to consider that this paper only addresses the efficient computational approximation (and the error related to it) of the Henssge method for time of death estimation given its published technical aspects, as it is typically used in the early postmortem period in forensic medicine.
Here, we do not address any potential errors related to the application of the Henssge method, These may, for example, be rooted in temperature measurement, body weight estimation. We also do not address any potential methodological errors inherent in the design or description of the method. These may, for example, stem from issuing a confidence interval based on two standard deviations as containing 95% (implied: 95,0 rather than ∼95,45%) [12, 14, 15], or from numerically simplifying the initially published 23,3^{∘}Ccutoff for application of the correct equation [12, 15] to 23^{∘}C (implied: 23,0) [16]. One study had examined the validity of the Henssge method and reported a violation of the predicted 95%confidence intervals in 57,1% of their 84 cases. Also, they reported larger standard deviations from their own experimental data and a trend to overestimate the Henssge method based postmortem interval in cases with high body mass and large body surface area [17].
Nevertheless, Henssge defined a scientifically based method for the estimation of postmortem intervals, with a number of case specific variables to carefully consider, that is popular and accepted in the forensic medicine community. For the correct application of the method in any given case, original literature should be consulted [1–3, 12, 15, 16, 18, 19].
Conclusions
While mathematically impossible to solve in a straightforward way, solutions to the Henssge equation can be approximated computationally by discretization with negligible error.
Availability of data and materials
Not applicable.
Notes
 1.
http://www.swisswuff.ch/calculators/todeszeit.php
 2.
http://www.swisswuff.ch/calculators/todeszeit.php
Abbreviations
 PHP:

Hypertext preprocessor
References
 1
Henssge C, Madea B. Estimation of the time since death in the early postmortem period. Forensic Sci Int. 2004; 144(2):167–75.
 2
Madea B, Brinkmann B. Handbuch Gerichtliche Medizin. 2 Bde. Berlin, Heidelberg: Springer; 2003.
 3
Saukko P, Knight B. Knight’s Forensic Pathology Fourth Edition. Boca Raton: CRC press; 2015.
 4
Marshall T, Hoare F. Estimating the time of death: the use of the cooling formula in the study of postmortem body cooling. J Forensic Sci. 1962; 7:189–210.
 5
Burger E, Dempers J, Steiner S, Shepherd R. Henssge nomogram typesetting error. Forensic Sci Med Pathol. 2013; 9(4):615–7.
 6
Smith E. Principles and analysis of approximation techniques. PhD thesis. 2016.
 7
Schöning U. Algorithmics in exponential time. In: Annual Symposium on Theoretical Aspects of Computer Science. Springer: 2005. p. 36–43.
 8
Cheney W. Discretization In: Cheney W, editor. Analysis for Applied Mathematics. New York: Springer: 2001. p. 170–6.
 9
Henssge C. Rectal temperature time of death nomogram: dependence of corrective factors on the body weight under stronger thermic insulation conditions. Forensic Sci Int. 1992; 54(1):51–66.
 10
Henssge C. Rectal temperature time of death nomogram: dependence of corrective factors on the body weight under stronger thermic insulation conditions. Forensic Sci Int. 1992; 54(1):51–66.
 11
Henssge C, Althaus L, Bolt J, Freislederer A, Haffner HT, Henssge C, Hoppe B, Schneider V. Experiences with a compound method for estimating the time since death. Int J Legal Med. 2000; 113(6):303–19.
 12
Henßge C. Todeszeitschätzungen durch die mathematische beschreibung der rektalen leichenabkühlung unter verschiedenen abkühlungsbedingungen. Zeitschrift für Rechtsmedizin. 1981; 87(3):147–78.
 13
Henßge C. Todeszeitbestimmung an leichen. Rechtsmedizin. 2002; 12(2):112–32.
 14
Stipanits E, Henßge C. Präzisionsvergleich von todeszeitrückrechnungen aus der rektaltemperatur ohne und mit berücksichtigung von einflussfaktoren. Beitr Gerichtl Med. 1985; 43:323–9.
 15
Henßge C. Die präzision von todeszeitschätzungen durch die mathematische beschreibung der rektalen leichenabkühlung. Zeitschrift für Rechtsmedizin. 1979; 83(1):49–67.
 16
Henssge C. Death time estimation in case work. i. the rectal temperature time of death nomogram. Forensic Sci Int. 1988; 38(34):209–36.
 17
Hubig M, Muggenthaler H, Sinicina I, Mall G. Temperature based forensic death time estimation: the standard model in experimental test. Legal Med. 2015; 17(5):381–7.
 18
Henssge C, Madea B, Gallenkemper E. Death time estimation in case work. ii. integration of different methods. Forensic Sci Int. 1988; 39(1):77–87.
 19
Henssge C, Madea B. Estimation of the time since death. Forensic Sci Int. 2007; 165(2):182–4.
Acknowledgements
The authors express their gratitude to Emma Louise Kessler, MD for her generous donation to the Zurich Institute of Forensic Medicine, University of Zurich, Switzerland.
Funding
No external funding was used for this.
Author information
Affiliations
Contributions
WS drafted the paper. MT contributed the institutional setting with the data for validation. Both authors read and agreed to the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Schweitzer, W., Thali, M.J. Computationally approximated solution for the equation for Henssge’s time of death estimation. BMC Med Inform Decis Mak 19, 201 (2019). https://doi.org/10.1186/s129110190920y
Received:
Accepted:
Published:
Keywords
 Computational approximation
 Estimation of early postmortem interval
 Forensic pathology
 PHP