狂犬病潛伏期及病例暴露后生存與風險變化的Cox回歸分析

孫建偉,陳豪敏,許汴利

狂犬病潛伏期及病例暴露后生存與風險變化的Cox回歸分析

孫建偉,陳豪敏,許汴利

目的分析和探討狂犬病潛伏期及病例暴露后生存與風險隨時間變化的特點。方法對2004-2009年河南省狂犬病病例進行流行病學調(diào)查,用SPSS15.0建立數(shù)據(jù)庫,進行Cox回歸等統(tǒng)計分析。結果狂犬病潛伏期為正偏態(tài)分布,調(diào)查顯示平均潛伏期為60 d。病例暴露后總的生存率變化趨勢是,40 d時病例累計生存率約為50%,60 d時累計生存率約為35%,80 d時累計生存率約為10%,120 d時的累計生存率幾近為零。累計風險量隨暴露后時間逐漸增高,120 d時達到風險高峰,之后形成一風險平臺。暴露后60 d以內(nèi)男性病例累計生存率要高于女性病例,60 d以后則低于女性。個案調(diào)查資料分析結果顯示,在暴露后50 d以內(nèi),男性累計生存率均高于女性,且30 d以內(nèi)差異有顯著性(P＜0.05)。50 d之后,男性累積生存率逐漸低于女性,第70 d時差異有顯著性(P＜0.05)。結論狂犬病潛伏期為正偏態(tài)分布,病例暴露后120 d時發(fā)病累積風險達到最高,且60 d內(nèi)女性發(fā)病風險高于男性。

狂犬病;潛伏期;生存率;風險量;Cox回歸

狂犬病在全球廣泛流行,每年至少有5.5萬人發(fā)病死亡,是迄今為止人類病死率最高的疾病[1]。狂犬病潛伏期是病例暴露起至發(fā)病這段時間,其長短變化幅度很大,短者數(shù)天,長者數(shù)年。暴露后規(guī)范的預防處置可大大降低狂犬病暴露著的發(fā)病風險[2]。為深入研究狂犬病潛伏期,探索病例暴露后生存與風險隨時間變化的規(guī)律,為狂犬病理論研究和預防控制提供支持,本文對有關狂犬病個案調(diào)查資料進行了分析和研究。

1 材料與方法

1.1 資料 本文材料來自河南省2004-2009年狂犬病病例的流行病學調(diào)查,對其中468例符合研究目的的個案進行分析。調(diào)查表采用全國狂犬病監(jiān)測方案病例個案調(diào)查表,調(diào)查內(nèi)容主要包括潛伏期、性別、年齡、職業(yè)、暴露方式、暴露程度、暴露部位等18個項目。

1.2 研究方法 用SPSS15.0建立數(shù)據(jù)庫并進行統(tǒng)計分析。用Cox回歸對潛伏期數(shù)據(jù)及病例暴露后累計生存率和風險量變化進行分析。在Cox回歸分析中,以潛伏期(d)為“生存時間”,以“病死”為事件唯一結局,以個案調(diào)查中性別、年齡、職業(yè)、暴露方式、暴露程度、出血與否、暴露部位、傷口處理、處理地點、是否擠壓排血、是否沖洗傷口、是否傷口消毒、是否接種疫苗、肇事動物管理、肇事動物狀態(tài)、動物傷人原因16個因素為自變量,進行生存回歸分析,探索病例暴露后生存與發(fā)病風險的動態(tài)變化。

2 結 果

2.1 潛伏期及其分布 潛伏期最短為4d,最長為2465天(6年9個月),平均潛伏期(中位數(shù))為60d。第25及75百分位數(shù)分別為38d、120d,1年內(nèi)發(fā)病的占到90.5%。狂犬病潛伏期呈正偏態(tài)分布,其偏度為5.047、峰度為31.099,進行常用對數(shù)轉換后,分布近似于正態(tài)分布,其偏度為 0.727、峰度為0.792,雖已較轉換前更接近正態(tài)分布,但經(jīng)單樣本K-S正態(tài)性檢驗(Z=2.054,P=0.000),仍不服從統(tǒng)計學意義上的正態(tài)分布。潛伏期轉換前后頻數(shù)分布及正態(tài)性Q-Q分布見圖1、圖2。

圖1 狂犬病潛伏期頻數(shù)分布及正態(tài)性Q-Q圖Fig.1 Graph of frequency distribution amd normality Q-Q test on human rabies incubation period

2.2 病例暴露后生存與風險變化 利用Cox回歸生存分析,對病例暴露后生存與發(fā)病變化進行探索,可以估計或預測人群暴露后累計生存率、累計風險量變化趨勢。病例暴露后總的生存率變化趨勢是,40 d時病例累計生存率約為50%,60 d時累計生存率約為35%,80 d時累計生存率約為10%,120 d時的累計生存率幾近為零。累計風險量隨暴露后時間逐漸增高,120 d達到風險高峰,之后形成一風險平臺。暴露后60 d以內(nèi)男性病例累計生存率要高于女性病例,60 d以后則低于女性,即暴露后60 d內(nèi)女性病例累計風險量高于男性病例,60 d后則低于男性病例。見圖3、圖4。

圖2 常用對數(shù)轉換后的狂犬病潛伏期時間分布及正態(tài)性Q-Q圖Fig.2 Graph of frequency distribution amd normality Q-Q test on human rabies incubation period after transformation by common logarithm

對468份個案分性別進行暴露后累積生存率分析,結果顯示,在暴露后 50 d以內(nèi),男性累計生存率均高于女性,且30 d以內(nèi)差異有顯著性。50天之后,男性累積生存率逐漸低于女性,第70 d時差異有顯著性。和上述模型預測結果相比相差10 d左右。可見,該模型能較為準確地預測暴露后生存率與風險量變化的一般規(guī)律或趨勢。見表1。

3 討 論

表1 病例分性別暴露后不同時間生存情況分析Table 1 The survival rate dynamics after exposure to human rabies by sex

一種疾病潛伏期的數(shù)據(jù)分布,取決于人群暴露后發(fā)病風險的動態(tài)變化。通常,這里有三種情況,一是暴露后發(fā)病風險快速升高,隨后緩慢下降,這種情況下潛伏期呈正偏態(tài)分布;二是暴露后發(fā)病風險逐漸升高,到達高峰后又逐步下降,這種情況下潛伏期呈對稱分布;三是暴露后發(fā)病風險持續(xù)在一較低水平一段時間,然后快速升高,達到高峰后又較快下降,這種情況下潛伏期呈負偏態(tài)分布。根據(jù)本研究中Cox回歸模型所得暴露后狂犬病發(fā)病累積風險示意圖,暴露后其發(fā)病風險是快速增高,然后維持一平臺期,故狂犬病潛伏期應呈正偏態(tài)分布,這在本研究調(diào)查中得到證實。狂犬病潛伏期中位數(shù)相關報道略有差異,這與所用統(tǒng)計方法、數(shù)據(jù)質(zhì)量、地區(qū)差異等多種原因有關[3-5]。狂犬病潛伏期總的來說取決于三個因素,一是暴露后進入機體的病毒數(shù)量與病毒毒力,二是暴露部位神經(jīng)組織是否豐富及與中樞神經(jīng)系統(tǒng)的距離,三是暴露者的免疫狀況。第一個因素與肇事動物傷人時的狀態(tài)(傳染力)、暴露方式、暴露程度、傷口處理情況等密切相關。第二個因素與暴露部位、年齡等有關。第三個因素涉及兩個方面,一是機體自身的免疫狀態(tài),與性別、年齡等有關;二是暴露后人工自動免疫(接種狂犬疫苗)與人工被動免疫(注射抗狂犬病免疫球蛋白)實施情況。暴露后病例性別的風險差異客觀存在,且模型分析結果得到了很好的實際調(diào)查數(shù)據(jù)支持。發(fā)病風險性別差異的原因可能與所研究病例性別差異有關,也可能與對狂犬病毒抵抗力不同有關,這是需要進一步重點關注和研究的方向。Cox回歸是半?yún)?shù)統(tǒng)計分析模型,對所關注數(shù)據(jù)分布類型沒有嚴格限制和要求,通常較多地應用在慢性疾病生存狀況研究中[6],在傳染病及暴露后發(fā)病風險與生存狀況研究中少有運用。在本文狂犬病病例暴露后生存與風險變化研究中,Cox回歸得出了與實際調(diào)查分析基本一致的結論,表明Cox回歸可以用于傳染病暴露后發(fā)病風險變化的研究,且可以得出研究對象隨時間變化的累計生存率與累計風險量關系圖,直接、形象展示出研究者所關心的事件結局動態(tài)變化。

[1]WHO.WHO Expert consultation on rabies first report[R].Geneva,Switzerland,2004:12-15.

[2]中華人民共和國衛(wèi)生部.狂犬病暴露預防處置工作規(guī)范[M].北京:衛(wèi)生部,2009.

[3]唐青,李浩.中國狂犬病流行近況及相關因素分析[J].中華流行病學雜志,2005,26(3):223-224.

[4]孫立梅,吳文鋒,許桂鋒.138例狂犬病潛伏期特點及其影響因素研究[J].實用預防醫(yī)學,2006,13(3):563-564.

[5]劉福強,郭綬衡,高立冬.1059例狂犬病個案分析及潛伏期影響因素研究[J].中國自然醫(yī)學雜志,2008,10(4):263-268.

[6]孫振球.醫(yī)學統(tǒng)計學[M].2版.北京:人民衛(wèi)生出版社,2007,256-268.

Incubation period and survival-hazard tendency after exposure to human rabies by Cox regression analysis

SUN Jian-wei,CHEN Hao-min,XU Bian-li
(Henan Provincial Center for Disease Control and Prevention,Zhengzhou450016,China)

The data of 468 human rabies cases from 2004 to 2009 in Henan Province were investigated and analyzed with SPSS15.0 in order to explore and analyze the incubation period and survival-hazard tendency after exposure of human rabies.The results showed that the data of rabies incubation period was positively skew distribution,and the median was sixty-day.The results of Cox model analysis were as follows:the tendency of hazard after exposure of female human rabies cases was higher than male cases within sixty-day,and the accumulated survival rate after exposure in 40thday was 50%and 35%in 60thday,10%in 80thday,and about zero in 120thday.The results of survey on 468 cases were as follows:the cumulative survival rate of male cases was higher than female cases within fifty-day after exposure,and was significant within thirty-day(P＜0.05).However,the rate was lower than female after 50thday,and was significant after 70thday(P＜0.05).In conclusion,the incidence hazard of female rabies after exposure is higher than male cases within sixty-day.Cox regression model is the best method for the analysis on variables of“survival time” .

human rabies;incubation period;survival rate;hazard rate;Cox regression

R181.2

A

1002-2694(2011)08-0734-04

河南省疾病預防控制中心,鄭州 450016;

Email:sunruohan@163.com

2011-02-17;

2011-05-26