生存分析PPT學習教案

1、會計學1 生存分析生存分析 疾病確診 死亡 治療開始 痊愈 手 術 復發(fā) 接觸毒物 出現(xiàn)毒性反映 吸 毒 HIV感染 HIV感染 自殺未遂 出 生 發(fā)病 第1頁/共27頁 第2頁/共27頁 圖 102名黑色素瘤患者的生存時間分布 time (Month) 051015202530354045 0 10 20 30 40 第3頁/共27頁 生存時間的分布為非正態(tài)分布 對數(shù)正態(tài)分布 指數(shù)分布， Weibull分布 Gamma分布，logistic分布， 含有刪失數(shù)據(jù)(censored data) 刪失數(shù)據(jù)提供的信息是不完全的 (uncompleted data) 數(shù)據(jù)特點 第4頁/共27頁 T39

2、1520202630 Y1011010 T：生存時間；：生存時間； Y：刪失指示量，：刪失指示量，Y1表示未刪失；表示未刪失；Y0表示刪失。表示刪失。 第5頁/共27頁 主要內(nèi)容 生存分析中的基本概念生存分析中的基本概念 生存率的估計和生存曲線生存率的估計和生存曲線 生存曲線的生存曲線的log-rank檢驗檢驗 Cox比例風險回歸模型比例風險回歸模型 第6頁/共27頁 生存率的估計 小樣本時，不需要對生存時間進行分組，小樣本時，不需要對生存時間進行分組， 常用常用 乘積極限法乘積極限法(product-limited，PL)，也稱也稱 Kaplan-Meier(K-M)法法估計生存率。估計生存

3、率。 大樣本時，需要對生存時間進行分組，常大樣本時，需要對生存時間進行分組，常 用壽命表用壽命表(lifetime table)法估計生存率。法估計生存率。 屬于非參數(shù)方法屬于非參數(shù)方法 第7頁/共27頁 第8頁/共27頁 01224364860 t (month) 0.0 0.2 0.4 0.6 0.8 1.0 S(t) treat control 圖29.5兒童橫紋肌肉瘤加化療組(treat)與對照組(control)緩解曲線比較 第9頁/共27頁 1 )( 2 2 2 v v ea j jj L 第10頁/共27頁 繪制生存曲 線 PL法估計生存曲 線 第11頁/共27頁 第12頁/共2

4、7頁 單因素分析方法單因素分析方法 Log-rank 檢驗 多因素分析方法多因素分析方法 Cox回歸分析回歸分析 第13頁/共27頁 Cox回歸分析(比例危險率回歸) 英國統(tǒng)計學家英國統(tǒng)計學家D.R.Cox于于1972年提出的年提出的比例比例 危險率回歸模型危險率回歸模型(proportional hazards regression model)，簡稱，簡稱Cox回歸?；貧w。 主要用于腫瘤與其他慢性病的預后分析，也可用主要用于腫瘤與其他慢性病的預后分析，也可用 于隊列研究的病因探索于隊列研究的病因探索 第14頁/共27頁 01122mm h(t,x)h (t)exp(xxx ) 半?yún)?shù)模型？

5、 第15頁/共27頁 1、建立、建立偏似然函數(shù)偏似然函數(shù)用用最大似然估計最大似然估計方法估計系數(shù)方法估計系數(shù) 2、回歸系數(shù)的檢驗方法、回歸系數(shù)的檢驗方法 似然比檢驗似然比檢驗 Score 檢驗檢驗 Wald 檢驗檢驗 回歸系數(shù)的估計與檢驗 第16頁/共27頁 少表示病人死亡的風險減取值越大，時， 無影響的取值對時， 大表示病人死亡的風險越取值越大，時， ,),(0 ),(0 ,),(0 XthX XthX XthX ii ii ii ，危險無關因素 危險因素危險因素 保護因素保護因素 exp() j RRb 回歸系數(shù) j的含義 第17頁/共27頁 第18頁/共27頁 第19頁/共27頁 調(diào)用調(diào)

6、用Cox模型分析模塊模型分析模塊 ，分析數(shù)據(jù)集為，分析數(shù)據(jù)集為cox 建立生存時間為建立生存時間為t ，截尾指示變，截尾指示變 量為量為d 的的Cox模型，模型，d 取值為取值為1時時 表示截尾表示截尾 選入和剔出水平均為選入和剔出水平均為 0.05 計算篩選因素的相對危險計算篩選因素的相對危險 度及其度及其95%的可信區(qū)間的可信區(qū)間 第20頁/共27頁 Summary of Stepwise Selection Variable Number Score Wald Step Entered Removed In Chi-Square Chi-Square Pr ChiSq 1 x4 1 13

7、.0399 . 0.0003 2 x5 2 4.7039 . 0.0301 篩選過程總結篩選過程總結 Analysis of Maximum Likelihood Estimat Parameter Standard Hazard 95%HR Variable DF Estimate Error Chi-Square PrChiSq Ratio Confidence Limits x4 1 1.76128 0.54785 10.3356 0.0013 5.820 1.989 17.031 x5 1 0.93133 0.44455 4.3890 0.0362 2.538 1.062 6.066

8、最大似然估計結果最大似然估計結果 治療方式（治療方式（x4）：相對危險度為：相對危險度為5.820，說明傳統(tǒng)的治療方式和新的治療方式相，說明傳統(tǒng)的治療方式和新的治療方式相 比，病人死亡的風險為比，病人死亡的風險為5.820倍，相對危險度的倍，相對危險度的95%的可信區(qū)間為的可信區(qū)間為1.98917.031； 淋巴結是否轉移（淋巴結是否轉移（x5）的）的RR值的含義與治療方式相同。值的含義與治療方式相同。 第21頁/共27頁 例例3 3 某醫(yī)藥公司與中山醫(yī)科大學協(xié)作，研究天花粉注射治某醫(yī)藥公司與中山醫(yī)科大學協(xié)作，研究天花粉注射治 療絨癌的療效問題。將療絨癌的療效問題。將1616只體表接種絨癌成功

9、的裸鼠隨機分只體表接種絨癌成功的裸鼠隨機分 為為4 4組，其中一組作空白對照，其余三組分別注射天花粉、甲組，其中一組作空白對照，其余三組分別注射天花粉、甲 藥和乙藥，經(jīng)整理，數(shù)據(jù)如表所示。試作藥和乙藥，經(jīng)整理，數(shù)據(jù)如表所示。試作CoxCox回歸分析?；貧w分析。 第22頁/共27頁 proc phreg; model day*d(0)=td v0 tr1tr2 tr3/selection=forward; run; 第23頁/共27頁 h(t,x)=h0(t)exp(0.42td-2.94tr1) 1、其他因素保持不變的情況下，帶瘤天數(shù)具有統(tǒng)計學意義（ x2=6.65，P0.05 ），帶瘤天數(shù)每增加一個單位，死亡風險增 加0.522倍。 2、其他因素保持不變的情況下，注射天花粉具有統(tǒng)計學意義 （x2=7.53，P0.05 ），注射天花粉者死亡風險是不注射者的 0.053倍。 第24頁/共27頁 例：某醫(yī)院泌尿外科醫(yī)師選擇例：某醫(yī)院泌尿外科醫(yī)師選擇1996-2000年間經(jīng)手術治療的年間經(jīng)手術治療的 膀胱腫瘤患者，對可能影響膀胱腫瘤術后生存的因素進行了膀胱腫瘤患者，對可能影響膀胱腫瘤術后生存的因素進行了 調(diào)查，研究者不僅要看是否出現(xiàn)了感興趣的終點事件或結局調(diào)查，研究者不僅要看是否出現(xiàn)了感興趣的終點事件或結局 （死亡）