狂犬病潛伏期及病例暴露后生存與風(fēng)險(xiǎn)變化的Cox回歸分析

孫建偉,陳豪敏,許汴利

狂犬病潛伏期及病例暴露后生存與風(fēng)險(xiǎn)變化的Cox回歸分析

孫建偉,陳豪敏,許汴利

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

狂犬病;潛伏期;生存率;風(fēng)險(xiǎn)量;Cox回歸

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

1 材料與方法

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

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

2 結(jié) 果

2.1 潛伏期及其分布 潛伏期最短為4d,最長為2465天(6年9個(gè)月),平均潛伏期(中位數(shù))為60d。第25及75百分位數(shù)分別為38d、120d,1年內(nèi)發(fā)病的占到90.5%?？袢摲诔收珣B(tài)分布,其偏度為5.047、峰度為31.099,進(jìn)行常用對(duì)數(shù)轉(zhuǎn)換后,分布近似于正態(tài)分布,其偏度為 0.727、峰度為0.792,雖已較轉(zhuǎn)換前更接近正態(tài)分布,但經(jīng)單樣本K-S正態(tài)性檢驗(yàn)(Z=2.054,P=0.000),仍不服從統(tǒng)計(jì)學(xué)意義上的正態(tài)分布。潛伏期轉(zhuǎn)換前后頻數(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 病例暴露后生存與風(fēng)險(xiǎn)變化 利用Cox回歸生存分析,對(duì)病例暴露后生存與發(fā)病變化進(jìn)行探索,可以估計(jì)或預(yù)測人群暴露后累計(jì)生存率、累計(jì)風(fēng)險(xiǎn)量變化趨勢(shì)。病例暴露后總的生存率變化趨勢(shì)是,40 d時(shí)病例累計(jì)生存率約為50%,60 d時(shí)累計(jì)生存率約為35%,80 d時(shí)累計(jì)生存率約為10%,120 d時(shí)的累計(jì)生存率幾近為零。累計(jì)風(fēng)險(xiǎn)量隨暴露后時(shí)間逐漸增高,120 d達(dá)到風(fēng)險(xiǎn)高峰,之后形成一風(fēng)險(xiǎn)平臺(tái)。暴露后60 d以內(nèi)男性病例累計(jì)生存率要高于女性病例,60 d以后則低于女性,即暴露后60 d內(nèi)女性病例累計(jì)風(fēng)險(xiǎn)量高于男性病例,60 d后則低于男性病例。見圖3、圖4。

圖2 常用對(duì)數(shù)轉(zhuǎn)換后的狂犬病潛伏期時(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

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

3 討 論

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

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

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

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

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

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

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

[6]孫振球.醫(yī)學(xué)統(tǒng)計(jì)學(xué)[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

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

Email:sunruohan@163.com

2011-02-17;

2011-05-26