《生存分析醫(yī)學(xué)生》

引例現(xiàn)用新舊兩種方法治療某癌病患者各100人，結(jié)果如下表，試評價其療效表3

新舊兩種方法治療某癌病

死亡人數(shù)新法組舊法組RR一年16151.01二年2535

0.70五年40800.50*2023/4/91編輯ppt生存分析2023/4/92編輯ppt一、基本概念生存時間（survivaltime)：疾病治療的預(yù)后情況，一方面看結(jié)局好壞，另一方面還要看出現(xiàn)這種結(jié)局所經(jīng)歷的時間長短。所經(jīng)歷的時間稱為生存時間。完全與不完全數(shù)據(jù)一部分研究對象可觀察到死亡，從而得到準確的生存時間，所提供的信息是完全的，稱為完全數(shù)據(jù)；另一部分病人由于失訪、意外事故、或到觀察結(jié)束時仍存活等原因，無法知道確切的生存時間，它提供了不完全的信息，稱為不完全數(shù)據(jù)（截尾數(shù)據(jù)、刪失數(shù)據(jù):censordata）。2023/4/93編輯ppt始點終點2023/4/94編輯ppt始點終點2023/4/95編輯ppt生存分析(survivalanalysis)：生存時間一般是通過隨訪收集。不完全數(shù)據(jù)提供了部分信息。須要用專門的方法進行統(tǒng)計處理，這類統(tǒng)計方法起源于對壽命資料的統(tǒng)計分析，故稱為生存分析。

2023/4/96編輯ppt

“生存”的概念

生物生存與死亡

儀器始使正常與出現(xiàn)故障疾病產(chǎn)生與治愈疾病治愈與復(fù)發(fā)

陰性與陽性 2023/4/97編輯ppt

截尾數(shù)據(jù)的處理

因為不太好處理截尾數(shù)據(jù)，很多臨床研究工作者常常將失訪或中止等原因造成的截尾數(shù)據(jù)在分析時拋棄。截尾數(shù)據(jù)提供的信息雖然是不完全的，但也很有價值，不應(yīng)隨便刪掉它。2023/4/98編輯ppt

生存分析的主要內(nèi)容第一，描述生存過程

研究生存時間的分布特點，估計生存率，生存曲線；第二，比較生存過程（假設(shè)檢驗）對兩組或多組生存率進行比較；第三，影響生存時間的因素分析了解影響生存過程的主要因素為改善預(yù)后提供指導(dǎo)。2023/4/99編輯ppt例在對資料進行描述時：

5名癌癥患者存活時間（月）

610142020n=5平均生存時間,mean=18，median=１４78+2535

+50

當(dāng)有截尾數(shù)據(jù)時，?2023/4/910編輯ppt

死亡概率

(mortalityprobability)是指死于某時段內(nèi)的可能性大小.記為q，年死亡概率計算公式為：某年內(nèi)死亡數(shù)

q＝───────某年年初人口數(shù)

若年內(nèi)有刪失，則分母用校正人口數(shù)：校正人口數(shù)＝年初人口數(shù)－1/2刪失數(shù)

二、生存資料常計算的幾種率2023/4/911編輯ppt死亡概率（2）死亡概率是指已活滿t時刻的個體，在此后一段時期內(nèi)（t至）死亡的可能性。死亡概率的計算公式為2023/4/912編輯ppt若t至內(nèi)有截尾值，上式的死亡概率也偏低，可對分母加以調(diào)整，為：如果很短，即時，死亡概率為t時刻的瞬時死亡概率，或稱為危險函數(shù)（hazardfunction,h(t)），其意義為剛剛活過t時刻的個體在t時刻死亡的概率。2023/4/913編輯ppt

生存概率

(survivalprobability)表示在某單位時段開始時存活的個體到該時段結(jié)束時仍存活的可能性大小,記為p。年生存概率計算公式為：活滿某一年人數(shù)

p＝1－q＝────────某年年初人口數(shù)若年內(nèi)有刪失，則分母用校正人口數(shù)。2023/4/914編輯ppt若在t至?xí)r段內(nèi)有截尾值，也應(yīng)計算調(diào)整生存概率生存概率與死亡概率之間有以下關(guān)系:p=1-q2023/4/915編輯ppt生存率與生存函數(shù)累積生存概率，簡稱生存率，是指病人經(jīng)過t個時間單位之后仍存活的概率。若令s(t)為任意時刻t的生存率，T為生存期，得s(t)=p(T≥t)0≤t<∞，稱s(t)為生存率函數(shù)，簡稱生存函數(shù)2023/4/916編輯ppt三、生存率估計（描述）

醫(yī)學(xué)研究中的生存資料分布常常不規(guī)則、不確定或未知分布，因而常用非參數(shù)法估計生存率。根據(jù)樣本含量的大小可分別選擇乘積極限估計法（小樣本時）壽命表法（大樣本時）2023/4/917編輯ppt1.乘積限估計法

Product-limitmethod簡稱積限法或PL法,直接用概率乘法原理估計生存率，它是由統(tǒng)計學(xué)家Kaplan和Meier于1958年首先提出的,因此又稱為Kaplan-Meier法.基本思想：將生存時間（包括截尾數(shù)據(jù)）逐個由小到大依次排列，并對其中的每個死亡點進行死亡概率、生存概率和生存率估計。應(yīng)用條件：數(shù)據(jù)個數(shù)較少時。2023/4/918編輯ppt2.壽命表法

（Lifetablemethod）

簡稱LT法。

基本思想：隨訪時間劃分成若干個時間區(qū)間，t時刻的生存率為t時刻前各時間區(qū)間生存概率的乘積。

應(yīng)用條件：數(shù)據(jù)個數(shù)多，經(jīng)分組時應(yīng)用。2023/4/919編輯ppt3.Kaplan-meier法計算生存率此法計算生存率時，先將每個個體的生存時間按照由小到大的順序排列，排序時若截尾值與非截尾值的觀察時間相同，則規(guī)定非截尾值小于截尾值，排在截尾值之前。然后依次計算出各時段的死亡概率、生存概率，進而計算出從觀察開始至各時刻的生存率及其標(biāo)準誤。2023/4/920編輯ppt小樣本生存分析Kaplan-Meier方法在SPSS中的實現(xiàn)生存狀態(tài)：1死亡，0截尾2023/4/921編輯pptSPSS軟件實現(xiàn)方法

File→Open→相應(yīng)數(shù)據(jù)文件→Analyze→Survival→Kaplan-Meier→Time(時間)→Status→Defineevent→singlevalue(1)→Continue→OK2023/4/922編輯ppt2023/4/923編輯ppt選擇選項2023/4/924編輯ppt2023/4/925編輯ppt例1：某臨床試驗對20名第III或第IV期黑色素瘤的患者進行隨訪研究，截至研究期結(jié)束，記錄的生存資料見下表，試計算100周的生存率。4.大樣本壽命表法2023/4/926編輯ppt表120名第III或第IV期黑色素瘤的患者的治療后的生存時間（周）12.8077.2015.6082.4024.00+87.20+26.4094.40+29.2097.20+30.80+106.00+39.20114.80+42.00117.20+58.40+140.00+72.00+168.00+2023/4/927編輯ppt建立數(shù)據(jù)文件“生存分析.sav”定義兩個變量：反應(yīng)變量，即生存時間變量，“time”分類變量，即生存狀態(tài)變量，“status”2023/4/928編輯pptAnalyze—Survival—LifeTables….2023/4/929編輯ppt2023/4/930編輯ppt本例100周生存率為53%。2023/4/931編輯ppt2023/4/932編輯ppt

三、生存過程的比較（假設(shè)檢驗）非參數(shù)方法:

對數(shù)秩檢驗(log-ranktest)

廣義Wilcoxon檢驗似然比（Likelihoodratio)檢驗

生存率組間比較實際上是對兩條或多條生存曲線的分布情況比較（齊性檢驗）

特點：檢驗過程一般不顧及生存率，而是利用死亡數(shù)和死亡率函數(shù)作統(tǒng)計推斷2023/4/933編輯ppt應(yīng)用范圍均為兩組或多組生存率比較的方法：對數(shù)秩檢驗Log－ranktest

：且近似服從Weibull分布，數(shù)值較小的失效時間所起的作用較大，（是SAS默認的選項檢驗方法）；廣義Wilcoxon檢驗：最適合數(shù)值較大的失效時間所起的作用較大。似然比（Likelihoodratio)檢驗：最適合生存資料近似服從指數(shù)分布。2023/4/934編輯ppt例題：22例非小細胞肺癌患者在不同日期經(jīng)隨機化分配到放療組和放化療聯(lián)合組，從緩解出院日開始隨訪，隨訪時間（月）見下表，試比較放療和放化療聯(lián)合兩種治療方案的生存率曲線有無差別。甲（放療組）：1，2，3，5，6，9+，11，13，16