楊宇雷,吳榮奇,費 耿,林中圣,施 群,肖 碧,王黎揚,葛延昌,孟航,馬開軍,陳新
(上海市公安局物證鑒定中心 上海市現場物證重點實驗室,上海200083)
尸溫相關因素推斷死亡時間的研究
楊宇雷,吳榮奇,費 耿,林中圣,施 群,肖 碧,王黎揚,葛延昌,孟航,馬開軍,陳新
(上海市公安局物證鑒定中心 上海市現場物證重點實驗室,上海200083)
目的 統計分析影響尸體溫度下降多因素和死亡時間的相關性,探索尸溫推斷死亡時間的創新應用。方法 建立影響尸體溫度(直腸溫度)下降相關因素的采集標準,收集實踐工作中的157例明確死亡時間真實案例的相關數據,根據實際工作經驗結合統計學的基本要求對各因素進行量化評分,利用EXCEL和SPSS軟件對數據進行處理,采用多元線性回歸的方法統計分析各影響系數同死亡時間的相關性。結果 獲得了具有統計學意義的回歸方程,Y=25.993+0.04X1+0.172X2+0.88X3+0.047X4+0.373X5+0.347X6-0.766X7,決定系數R2=0.876。結論 該方法為尸溫推斷死亡時間的創新應用,經測試可用于實際工作。
尸體溫度;死亡時間;影響因素;量化評分;多元線性回歸
Abstract:Objective To statistically analyze the correlation between factors related to body temperature drop and postmortem interval,and to establish an innovative method for estimating postmortem interval by body temperature.Methods The collection standard of the factors influencing body temperature drop was established.Data of 157 cases which had exact postmortem intervals were collected.The related factors were quantitatively scored according to practical experience and basic statistic requirements.The data were processed by EXCEL and SPSS software,and the correlations between the factors and the postmortem intervals were analyzed with multiple linear regression.Results A statistically significant regression equation was obtained:Y=25.993+0.04X1+0.172X2+0.88X3+0.047X4+0.373X5+0.347X6-0.766X7.The correlation coefficient R2was 0.876.Conclusion The study presented an innovation application of body temperature for estimating postmortem interval,which had been tested in the practical work.
Key words:body temperature;postmortem interval;related factor;quantitative score;multiple linear regression
死亡時間是指法醫檢查尸體時刻距死者死亡時刻的間隔時間,死亡時間推斷是法醫實際工作中需要首先解決的重要任務之一,也一直是法醫學鑒定的難點和研究的熱點[1]。由于死亡時間推斷的影響因素繁多,還尚無一種能完全克服各種影響因素或者忽略影響因素的死亡時間推斷方法[2],而在公安法醫的實際工作中,死亡時間的推斷往往要滿足于案件偵破等需要,快速、準確是一切技術手段的首要要求。
本研究就是以“簡便實用,快速準確”為目的,根據大量的文獻資料和工作經驗,選擇最有可能精確判斷死亡時間的尸體溫度變化規律,統籌考慮環境溫度、個體情況、衣著情況、載體情況、死亡原因、環境通風情況、濕度等影響因素,根據實際工作經驗創新地量化重要的影響因素,進而通過一系列的統計分析尋求回歸方程并加以驗證,尋找出影響尸溫下降相關因素的系數,以實現利用尸溫推斷死亡時間這個技術手段的應用創新。
2007~2011年,上海市公安局物證鑒定中心法醫室根據統一的規范收集的157例已知死亡時間的案例,并根據實際情況,列舉了實踐工作中常見的各種情況,按照統計分析的要求結合工作經驗進行量化評分。
(1)主要收集檢驗時測得的直腸溫度大于環境溫度的案例;
(2)詳細記錄案發當時的天氣情況;
(3)詳細記錄死者所處環境的實時溫度、載體情況、環境通風情況等。載體需要詳細描述材質;環境需要區分室內室外,空間大小,有無空調,有無可感風速等;
(4)尸體檢驗時詳細記錄個體情況、衣著情況、載體情況,死亡原因。個體情況包括死者年齡,性別,胖瘦等情況,其中胖瘦情況需記錄腹部脂肪最厚處的脂肪厚度;衣著情況需要詳細記錄衣服的長短袖情況,衣服的件數,布料的材質,有無填充物及其材料等;載體情況需要記錄詳細的載體名稱,主要分成泥石類、液體類、木質類、布棉類;死亡原因記錄大失血、顱腦損傷、機械性窒息等情況,大失血要區分內出血或者外出血;
(5)準確測量、詳細記錄尸體溫度,本研究采用直腸溫度表示。
(6)實際死亡時間的追蹤。根據案情反映、醫院病史、監控信息、證人證言、犯罪嫌疑人供述等綜合判斷分析確切的死亡時間。
(1)個體情況:著重考慮對尸體溫度下降影響最大的個體胖瘦情況,據有關研究報道,人體的胖瘦與腹壁的脂肪厚度成正相關性[3],本研究即采用腹壁脂肪的厚度來體現個體情況中人體的胖瘦。單位為:cm。
(2)環境溫度:以現場實測的溫度表示。單位為:℃。
(3)衣著情況、載體情況、死亡原因和通風情況的列舉如下(見表1):本研究所使用的各個案例情況均符合相應的法律法規情況和相關的保密條款。

表1 影響尸體溫度下降因素的列舉表
(1)環境溫度的測量。經計量單位校準過的水銀溫度計和煤油溫度計各一根,水銀溫度計放置于尸體邊,水銀球距周圍物體15cm以上,煤油溫度計置于尸體所處空間的中心位置,取兩次測量值的平均數。
(2)直腸溫度的測量。經計量單位校準過的水銀溫度計和煤油溫度計各一根,先由水銀溫度計經肛門向上插入10cm,放置5min,然后讀取數值,緊接著用煤油溫度計再測一次,取兩次測量值的平均數。
(1)個體情況:腹壁脂肪最厚處厚度的實際測量值。單位為:cm。
(2)環境溫度:以現場實測的溫度表示。單位為:℃。
根據實際工作經驗判斷,將具有近似保暖功能和熱傳導功能的衣著和載體歸類并同通風情況和死亡時間量化評分(見表2)。

表2 影響尸體溫度下降因素的評分表
將表2的評分數值代入157例實際案例中,X1=個體情況、X2=環境溫度、X3=衣著情況、X4=載體情況、X5=死亡原因、X6=環境通風情況、X7=直腸溫度,得出Y=死亡時間。
采用SPSS17.0統計軟件包對數據進行多元線性回歸分析。檢驗水準α=0.05。
復相關系數R=0.936,決定系數R2=0.876。
F=150.060,F 檢驗結果概率 p=0.000<0.001,差異有統計學意義,回歸模擬的擬合是有效的。
自相關問題的診斷 :DW值一般要求1.5~2.5h,殘差與自變量互為獨立。回歸模型的DW值為1.888,說明該模型無自相關的問題,此模型可以被使用。
根據以上回歸系數的檢驗及結果,建立回歸方程如下:

復相關系數R=0.936,決定系數R2=0.876。
Y為死亡時間,X1為個體胖瘦情況 (脂肪厚度表示),X2為環境溫度,X3為衣著情況,X4為載體情況,X5為環境通風情況,X6為死亡原因,X7為直腸溫度。
組內考核,誤差在1h以內的為60.6%,誤差在2h以內的為87.9%,誤差在3h以內的為95.5%;組外驗證3個案例,顯示2個案例誤差在1h以內,1個案例誤差在2h以內。
根據尸溫的下降速度來推測死亡時間是最早且較準確的方法之一,國內外法醫學家對其的研究一直沒有停止過,也找出了一定的規律性,但在實際檢案中往往做不到簡便而準確,要么是沒有充分考慮多種因素的影響而得出簡單的線性公式,要么是各種條件下多個多元回歸方程式[4-5],不利于實戰中快速準確地得到結果。尸體溫度下降產生影響的這六個因素里,各個因素在尸體溫度下降的過程中占有的比重究竟有多大,一直是困擾著我們的一個問題。
事實證明,尸體溫度下降的規律和影響因素之間的關系研究,僅僅通過實驗室研究的手段是完全不夠的。規律的掌握或者模擬不僅僅只能依靠實驗室工作,實踐證明有些規律雖然人們不知道他確切的聯系性,但是可以通過大量數據的統計研究,盡量接近或者模擬這種規律性。公安法醫工作的性質決定了我們有得天獨厚的收集各類數據及獲得確切死亡時間的優勢,于是就有了該課題的設想和實施。

表3 模型效果檢驗匯總2)

表4 回歸系數1)的檢驗及結果
根據統計分析的基本要求,影響尸體溫度下降的六個指標,即個體條件、環境溫度、衣著與履蓋物、載體、環境通風情況及死亡原因,除去個體中的胖瘦情況、環境溫度本身即以數值表示,其他因素需要進行數字化評分。每一類影響因素本身包含有多種情況,也許不勝枚舉,但在實際工作中我們發現每類影響因素中均有常見的情況,同時根據經驗判斷將有些對尸體溫度下降有較類似作用的常見情況進行歸類,然后再加以列舉,最后根據各種常見情況的影響力大小,依據工作中的經驗判斷對于每項進行賦值量化。
經采用SPSS軟件,多元逐步回歸死亡時間Y和個體情況 X1、環境溫度X2、衣著及覆蓋物情況X3、載體情況X4、環境通風情況X5、死亡原因X6、直腸溫度X7之間的相關性,得到回歸方程:

復相關系數R=0.936,決定系數R2=0.876。
從回歸方程的復相關系數和方差分析結果表明,本課題采用的死亡時間同影響因素之間的多元回歸方程是有統計學意義的。從組內考核和組外驗證情況來看,該回歸方程計算出的死亡時間能夠較準確反映出真實的死亡時間,且相應的誤差基本能夠在公安工作需要的允許范圍之內,具有一定的實際應用價值。
目前根據157例案例的統計分析,得出各影響因素對應的影響系數分別為:個體胖瘦情況0.040,環境溫度0.172,衣著情況0.88,載體情況0.047,環境通風情況0.373,死亡原因0.347,直腸溫度-0.766。其中個體胖瘦情況和載體情況的影響系數相對較小,僅僅從回歸方程本身來講,其他因素保持不變的情況下,個體胖瘦和載體情況的變動對結果的影響較小,該兩項影響系數似乎與實際情況和經驗判斷有一定誤差。分析其中的原因,我們考慮如下:(1)由于樣本量的限制,目前僅有157例,各影響因素的系數準確度還有待進一步提高;(2)根據實際情況,觀察我們收集的數據本身,在個體胖瘦情況和載體情況中這兩項評分中,數據均較集中,也即大部分個體胖瘦均較一致,且上海地區主要以泥石類為主要現場;(3)在評分環節,特別是載體類因素中,由于缺乏對熱傳導相關數據的證實,評分中未體現出不同類的具體差異。綜上,受制于樣本量等因素,本課題所得出的多元回歸方程,不可否認,還是存在著一定的誤差,在實際應用中需要借助其他推斷死亡時間的方法綜合判斷,能做出更接近真實死亡時間的推斷。
本課題的創新應用研究,不僅僅是為了獲得一個具有統計學意義且可為實踐所用的多元回歸方程,更重要的是結合公安法醫在收集數據方面得天獨厚的優勢,從課題設計思路上進行創新,而方程的獲得和檢驗結果的分析,更加證明本課題的設計是可以滿足尸體溫度下降因素同死亡時間之間的相關性研究,且該統計方法是有效的,能發揮作用的,是利用尸體溫度下降規律來推算死亡時間的創新應用。
為了使本課題所研究的推斷死亡時間的新方法能夠得到普遍應用,在繼續收集各種場所中各種已知明確死亡時間的案例,用于不斷補充統計研究的樣本量的同時,還需要對研究思路、自變量的選擇、自變量的評估方面不斷地完善。
[1]黃光照,麻永昌.中國刑事科學技術大全法醫病理學[M].北京:中國人民公安大學出版社,2002:116-133.
[2]C.Henssge,B,Madea.Estimation of the time since death.Forensic Sci Int[J].2007,(165):182-184.
[3]劉力劍,趙文海,榮梅.成人腹壁脂肪厚度與體重指數、腰圍的相關性分析[J].中國誤診學雜志,2009,(31):7577-7578.
[4]八省市(區)公安廳(局)尸溫研究小組.581例非正常死亡尸體直腸溫度與時間關系的研究[J].刑事技術,1984,(4):1-10.
[5]徐長苗,陳希生.推定死亡時間的新方法[J].刑事技術,1994,(3):23-25.
(本文編輯:秦志強)
Study on Estimating Postmortem Interval with Factors Related to Body Temperature
YANG Yu-lei,WU Rong-qi,FEI Geng,LIN Zhong-sheng,SHI Qun,XIAO Bi,WANG Li-yang,GE Yan-chang,MENG Hang,MA Kai-jun,CHEN Xin
(Shanghai Key Laboratory of Crime Scene Evidence,Institute of Forensic Science,Shanghai Public Security Bureau,Shanghai 200083,China)
DF795.4
A
10.3969/j.issn.1671-2072.2012.05.009
1671-2072-(2012)05-0075-04
2012-05-27
公安部應用創新計劃項目(2008YYCXSHSJ032)
楊宇雷(1978—)男,碩士研究生,主要從事法醫學研究工作。 E-mail:wangnianjiao2007@yahoo.com.cn。
馬開軍(1975—),男,副主任法醫師,主要從事法醫病理學和法醫臨床學研究工作。E-mail:makaijun@sina.cn。