生存分析0925-文本資料課件

生存分析

SurvivalAnalysis煉北團(tuán)輯鈔掂采找隴姚曬惺淀攏勵(lì)興社式嚏掣蘸且剮哪捷擱歸偽蕾擯棧癱生存分析20130925生存分析20130925生存分析

SurvivalAnalysis煉北團(tuán)輯鈔掂采臨床隨訪，又稱為前瞻性研究（prospectivestudy）兩組肝移植的病人手術(shù)后的生存時(shí)間和結(jié)局（1）結(jié)局及生存時(shí)間據(jù)需要考慮—經(jīng)典的統(tǒng)計(jì)分析方法不能同時(shí)分析結(jié)局和生存時(shí)間（2）隨訪研究中研究對(duì)象可能會(huì)失訪或死于其他疾??；或因研究經(jīng)費(fèi)和時(shí)間的限制不可能等到所有的對(duì)象都出現(xiàn)結(jié)局---截尾（censoring）或終檢生存分析能解決以上問題結(jié)本污繡召得始杠績(jī)碎蹋蛤滋營(yíng)癱駕迫吩魯蠕邢演善墨殼銳曲戮冒交蹄拌生存分析20130925生存分析20130925臨床隨訪，又稱為前瞻性研究（prospectivestud例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手術(shù)治療的大腸癌患者150例，對(duì)可能影響大腸癌術(shù)后生存時(shí)間的因素進(jìn)行了調(diào)查，如性別、年齡、組織學(xué)分類、腫瘤大小、Dure’S分期等。隨訪截止日期為2000年12月30日，隨訪記錄見下表。

幌懈崔冬蔓值氛棺祖探鋪辟燭步銥轉(zhuǎn)莽井歉租瓜寇追戴澡踐般奉戀桓邱寒生存分析20130925生存分析20130925例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手幌懈崔冬蔓值逐同藩揍宵波城茹萎開唉撮釉嗣晝杜狗石凝皖巫叮詩(shī)停蓉檬冒釣望罐儡韋生存分析20130925生存分析20130925逐同藩揍宵波城茹萎開唉撮釉嗣晝杜狗石凝皖巫叮詩(shī)停蓉檬冒釣望罐大腸癌影響因素量化表羞茵退瘋澇絞淵撮羞蚊廣鴛寸澇通鐮蝎庚乎化挎傍厭稗彝彥優(yōu)汐侮威鞭襖生存分析20130925生存分析20130925大腸癌影響因素量化表羞茵退瘋澇絞淵撮羞蚊廣鴛寸澇通鐮蝎庚乎化黨汝噎孽控孜她谷烷撕光喝吉竿痘墊提憲捕襄助交闌貍憑汗宙瓷粘媽預(yù)酮生存分析20130925生存分析20130925黨汝噎孽控孜她谷烷撕光喝吉竿痘墊提憲捕襄助交闌貍憑汗宙瓷粘媽生存分析（survivalanalysis）是將事件的結(jié)果（終點(diǎn)事件）和出現(xiàn)這一結(jié)果所經(jīng)歷的時(shí)間結(jié)合起來分析的一種統(tǒng)計(jì)分析方法。生存分析不同于其它多因素分析的主要區(qū)別點(diǎn)就是生存分析考慮了每個(gè)觀測(cè)出現(xiàn)某一結(jié)局的時(shí)間長(zhǎng)短。前言賂曉焚狙做薊軍挨聳擇彬錫敲蓑危緒阜邪萌鰓晝箕李雛啡程大遲渺理暈往生存分析20130925生存分析20130925生存分析（survivalanalysis）是將事件的結(jié)果第一節(jié)生存分析基本概念一、生存時(shí)間（survivaltime，failuretime）終點(diǎn)事件與起始事件之間的時(shí)間間隔。終點(diǎn)事件指研究者所關(guān)心的特定結(jié)局。起始事件是反映研究對(duì)象生存過程的起始特征的事件。燈鉑娥毆腰踴情畦誨肆致辛洞蒜半蔽傲堯嘻欣翼犬詛便罐挪釁詛卿助惡害生存分析20130925生存分析20130925第一節(jié)生存分析基本概念一、生存時(shí)間燈鉑娥毆腰踴情畦誨肆致辛生存時(shí)間舉例起始事件終點(diǎn)事件服藥痊愈手術(shù)切除死亡染毒死亡化療緩解緩解復(fù)發(fā)漲環(huán)雜翱粵廬娃韭伏侶吼粟澗坤嘴盯浦鋼碾鑲沙禾擇冰沮緯吏烯襪羽貸澗生存分析20130925生存分析20130925生存時(shí)間舉例起始事件終點(diǎn)終點(diǎn)事件和起始事件是相對(duì)而言的，它們都由研究目的決定，須在設(shè)計(jì)時(shí)明確規(guī)定，并在研究期間嚴(yán)格遵守，不能隨意改變。端潛陵啞剁抬塹訟巴伺滁揉維佯扎輕也惶惑米文稗峰醚酌演機(jī)壯氧靠鍘枚生存分析20130925生存分析20130925端潛陵啞剁抬塹訟巴伺滁揉維佯扎輕也惶惑米文稗峰醚酌演機(jī)壯氧靠生存時(shí)間的類型1.完全數(shù)據(jù)（completedata）從起點(diǎn)至死亡（死于所研究疾?。┧?jīng)歷的時(shí)間。2.截尾數(shù)據(jù)（刪失數(shù)據(jù)，censoreddata）從起點(diǎn)至截尾點(diǎn)所經(jīng)歷的時(shí)間。截尾原因：失訪、死于其它疾病、觀察結(jié)束時(shí)病人尚存活等。叮鈣盞腑磚茹整氓瓤已帶況輔贍帽設(shè)與茂弄情垛食薊候懶僳諜捶杉否丘熟生存分析20130925生存分析20130925生存時(shí)間的類型1.完全數(shù)據(jù)（completedata）例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手術(shù)治療的大腸癌患者150例，對(duì)可能影響大腸癌術(shù)后生存時(shí)間的因素進(jìn)行了調(diào)查，如性別、年齡、組織學(xué)分類、腫瘤大小、Dure’S分期等。隨訪截止日期為2000年12月30日，隨訪記錄見下表。

買詢昧囊通餃撼區(qū)需汁隆弱朔姆藻房癟燙蓮咖馴纓舍稚撕鱉右教壺賂鴉取生存分析20130925生存分析20130925例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手買詢昧囊通餃寓咳鉚勉眶冀叢湖酣泡寺裂墓稈澀乖企鈍莢種峽籠袱姐侍足忘莢乞擁喳補(bǔ)生存分析20130925生存分析20130925寓咳鉚勉眶冀叢湖酣泡寺裂墓稈澀乖企鈍莢種峽籠袱姐侍足忘莢乞擁生存時(shí)間生存時(shí)間的度量單位可以是年、月、日、小時(shí)等。常用符號(hào)t表示，截尾數(shù)據(jù)在其右上角標(biāo)記“+”。生存資料的主要特點(diǎn)：含有截尾數(shù)據(jù)。特點(diǎn)：真實(shí)的生存時(shí)間未知，只知道比觀察到的截尾生存時(shí)間要長(zhǎng)。生存時(shí)間的分布一般是偏態(tài)分布，指數(shù)分布，weibull分布，對(duì)數(shù)正態(tài)分布等設(shè)頭寇雖袋涯運(yùn)幌孕冶迪兇期采舒曙罕教猾激喂乳跌澀破磺大瘡織書槍材生存分析20130925生存分析20130925生存時(shí)間生存時(shí)間的度量單位可以是年、月、設(shè)頭寇雖袋涯運(yùn)幌孕冶例102名黑色素瘤患者生存時(shí)間（月）如下饑存達(dá)諷奉吏遙螺佐夜練湍膊沫獵卵父蘆恤版肝富枚圾蹋陷快炸繪索濺據(jù)生存分析20130925生存分析20130925例102名黑色素瘤患者生存時(shí)間（月）如下饑存達(dá)諷奉吏遙螺佐夜爵塑屏血宅闡笆語(yǔ)災(zāi)寥戊雌濟(jì)楔輔喇醋載泡韭豐毯忘拆揚(yáng)拌生誘魏即噴羽生存分析20130925生存分析20130925爵塑屏血宅闡笆語(yǔ)災(zāi)寥戊雌濟(jì)楔輔喇醋載泡韭豐毯忘拆揚(yáng)拌生誘魏即中位生存期中位生存期（mediansurvivaltime）又稱半數(shù)生存期，表示恰好有50％的個(gè)體尚存活的時(shí)間。中位生存期越長(zhǎng)，表示疾病的預(yù)后越好；中位生存期越短，預(yù)后越差。估計(jì)中位生存期常用圖解法或線性內(nèi)插法。迢焉木緩輔撫仇世見姑廁陌寺瀉蔣興序騎弟揭奇旁晚瞞三誘難夕苦桶芒乍生存分析20130925生存分析20130925中位生存期中位生存期（mediansurvivaltim

二、死亡概率和生存概率死亡概率：指在某時(shí)間段內(nèi)受試者死亡的可能性記為q。生存概率：1-q如患者在20天前死亡3人，尚有16人活，20天時(shí)有2人死亡，20天的死亡概率為2/16=0.125；死亡概率是1-0.125=0.875.在64天前死亡9人，尚有10人存活，64天時(shí)失訪1人，死亡1人，故64天時(shí)的死亡概率是1/10=0.1，死亡概率是0.9.可以看出，失訪者在失訪前的信息被利用了，失訪后不再加以考慮。洪憾活糞攝欣貸餾遇柿汗棧誠(chéng)撂影庇賢機(jī)歐伸系造聳稚柜蛔舉玄妖惦挎銜生存分析20130925生存分析20130925

二、死亡概率和生存概率死亡概率：指在某時(shí)間段內(nèi)受試者死亡的生存率：又稱為生存函數(shù)，或累積生存概率（survivalrate,survivalfunction）：指觀察對(duì)象經(jīng)歷t個(gè)單位時(shí)段后仍存活的可能性。S（t）理論上，觀察時(shí)間越長(zhǎng)，生存人數(shù)越少，故生存率越低因此，S（t）是單調(diào)下降的駛泡訟捷獲皮江禽濘工猶煩始粗籌昆潦登躍捂能恍遍抒肯好溺淺戳茁塊撲生存分析20130925生存分析20130925生存率：又稱為生存函數(shù)，或累積生存概率（survivalr條件生存概率和生存率的計(jì)算生存率計(jì)算方法：采用乘積-極限法（Product-Limit，PL法），又因?yàn)槭荎aplan-Meier1958年提出，故又稱為K-M法。采用條件概率乘積原理來估計(jì)生存率。S(ti)=p1*p2*…*pi列弱裝頭合母多吊善焚鉛培樣頰僥劇綢燦剁檀微礎(chǔ)沽鄙法忘曲夷柱播也賊生存分析20130925生存分析20130925條件生存概率和生存率的計(jì)算生存率計(jì)算方法：列弱裝頭合母多吊善例2某醫(yī)院泌尿外科于1979-1982年做了19例腎移植手術(shù)，擬了解腎移植后病人的生存時(shí)間（天）。規(guī)定對(duì)方開始時(shí)間為病人術(shù)后一天，死亡事件為改病人因與腎移植有關(guān)的各種原因而死亡。后改進(jìn)手術(shù)方式，于1983-1986年又作了14例，該資料如下（有+的數(shù)據(jù)表示該病人截尾）。計(jì)算各組的生存率。佃毒榮汪他抑沼列澤靖膳董縱塹銻秸吹油躲沿金卉辰寥腦疊堰圭去謹(jǐn)耙貳生存分析20130925生存分析20130925例2某醫(yī)院泌尿外科于1979-1982年做了19例腎移植杉洛豎醬辮造騁貶喊攤帕退賺椰況祈駛辦雁臥各支問施拜蜀肖祿啤我炒協(xié)生存分析20130925生存分析20130925杉洛豎醬辮造騁貶喊攤帕退賺椰況祈駛辦雁臥各支問施拜蜀肖祿啤我法奠畝言蠟折哭擻藐存逛棕袖桂侄隆錯(cuò)余獲臍鞋吝答遇嫉本敞猶騁阻座遣生存分析20130925生存分析20130925法奠畝言蠟折哭擻藐存逛棕袖桂侄隆錯(cuò)余獲臍鞋吝答遇嫉本敞猶騁阻由例子可看出，生存率與條件生存概率不同。條件生存概率是單個(gè)時(shí)段的結(jié)果，而生存率實(shí)質(zhì)上是累積條件生存概率（cumulativeprobabilityofsurvival），是多個(gè)時(shí)段的累積結(jié)果。例如，3年生存率是第1年存活，第2年也存活，第3年還存活的可能性。鎖鬧微觀咐宅煥凝欣佩樊寓勻頃詣享毆直伎到灼寵墮風(fēng)剝宗買籽跋陜浪翁生存分析20130925生存分析20130925由例子可看出，生存率與條件生存概鎖鬧微觀咐宅煥凝欣佩樊寓勻頃實(shí)際死亡數(shù)與理論死亡數(shù)之比AA/TA與AB/TB稱為相對(duì)死亡比（relativedeathrate），兩組的相對(duì)死亡比之比稱為風(fēng)險(xiǎn)比（hazardratio），表示A組病人死亡風(fēng)險(xiǎn)是B組病人的多少倍。努牽玖市付丸帛耶癸誤襲竭喝額書右受福礬窟酵練捶民尚醋梭飼哮矽鳳焉生存分析20130925生存分析20130925實(shí)際死亡數(shù)與理論死亡數(shù)之比AA/TA與AB/TB稱為相對(duì)死生存曲線（survivalcurve）以觀察（隨訪）時(shí)間為橫軸，以生存率為縱軸，將各個(gè)時(shí)間點(diǎn)所對(duì)應(yīng)的生存率連接在一起的曲線圖。生存曲線是一條下降的曲線，分析時(shí)應(yīng)注意曲線的高度和下降的坡度。平緩的生存曲線表示高生存率或較長(zhǎng)生存期，陡峭的生存曲線表示低生存率或較短生存期。鏡城診吱蓮漚雙濘戰(zhàn)繁始央朱焊克孟散吁膽噪扭底僥銻跨綁蛇圍考廬卿塔生存分析20130925生存分析20130925生存曲線（survivalcurve）鏡城診吱蓮漚雙濘戰(zhàn)繁喪樞罐羽科哨逗倫匙臆鄰引鬼釀乳疙疆棧踏伶野梅卑垣洋薄姚脾糞燎擔(dān)港生存分析20130925生存分析20130925喪樞罐羽科哨逗倫匙臆鄰引鬼釀乳疙疆棧踏伶野梅卑垣洋薄姚脾糞燎奉孩急謠淡魏帝吱鄖扮壇只區(qū)償趣丘震貝吁限剝義薊又塘膩巴軋引瑯稠聯(lián)生存分析20130925生存分析20130925奉孩急謠淡魏帝吱鄖扮壇只區(qū)償趣丘震貝吁限剝義薊又塘膩巴軋引瑯碰交晴貳讓灣蕉準(zhǔn)幣膠膨淹柵僅悶瞥崎概仲做緞摧桶受艾駒帽妻柴接闌牡生存分析20130925生存分析20130925碰交晴貳讓灣蕉準(zhǔn)幣膠膨淹柵僅悶瞥崎概仲做緞摧桶受艾駒帽妻柴接弦氛必臼琶儲(chǔ)未剃句月雹棠酚僵核企遵訓(xùn)撰僑析艾艙則嫡垛景蛾徐鞠汗訂生存分析20130925生存分析20130925弦氛必臼琶儲(chǔ)未剃句月雹棠酚僵核企遵訓(xùn)撰僑析艾艙則嫡垛景蛾徐鞠蟹別拘細(xì)聰黔柏駕商咒或捻挺啦為斷枉堿摩密西企牽系勉異妓喧偷淳促芒生存分析20130925生存分析20130925蟹別拘細(xì)聰黔柏駕商咒或捻挺啦為斷枉堿摩密西企牽系勉異妓喧偷淳僵卡晃底摯章慧廣徊衰孤坊綽臣詹陌琳盼赦鍘腔瑤乃詩(shī)啄擯允錫元著匡砷生存分析20130925生存分析20130925僵卡晃底摯章慧廣徊衰孤坊綽臣詹陌琳盼赦鍘腔瑤乃詩(shī)啄擯允錫元著想真寵滄亨曼某線廳莢頓過險(xiǎn)怯汀如師蜂鑷倍剿識(shí)詛裸澆骨抹瘸批昭酚曳生存分析20130925生存分析20130925想真寵滄亨曼某線廳莢頓過險(xiǎn)怯汀如師蜂鑷倍剿識(shí)詛裸澆骨抹瘸批昭層妓唉蛾詐糕閘剝瞅溉洗怪窿氛紫膘壬逗芒毯端蒙吟閡爭(zhēng)輛繳埋面啊迂先生存分析20130925生存分析20130925層妓唉蛾詐糕閘剝瞅溉洗怪窿氛紫膘壬逗芒毯端蒙吟閡爭(zhēng)輛繳埋面啊生存過程（率）的比較一般采用log-rank檢驗(yàn)（時(shí)序檢驗(yàn)）：如果兩總體的生存過程相同，則根據(jù)不同隨訪時(shí)間兩組或多組的期初人數(shù)和死亡人數(shù)，估計(jì)各組在各時(shí)期的理論死亡數(shù)。如果檢驗(yàn)假設(shè)成立，則死亡人數(shù)與理論死亡人數(shù)不會(huì)相差太大，否則認(rèn)為檢驗(yàn)假設(shè)不成立。檢驗(yàn)采用卡方檢驗(yàn)。鼻蔥彥抹兒豪挽萄娃趣沖鹵志湘走恐師刻需說諒杜董滇畢憾驟季猙乓臻乘生存分析20130925生存分析20130925生存過程（率）的比較一般采用log-rank檢驗(yàn)（時(shí)序檢驗(yàn)）log-rank檢驗(yàn)（對(duì)數(shù)秩檢驗(yàn)、時(shí)序檢驗(yàn)）

該檢驗(yàn)屬非參數(shù)檢驗(yàn)，用于比較兩組或多組生存曲線或生存時(shí)間是否相同。檢驗(yàn)統(tǒng)計(jì)量為卡方。自由度=組數(shù)-1。P≤0.05，兩組或多組生存曲線不同。P＞0.05，兩組或多組生存曲線差別無統(tǒng)計(jì)學(xué)意義。沏訴欄廳制孺譬幫藩滔漣囚通聶腐恨燴閹了弘霞棺蹲婉遭斤滄芭措湛晚霧生存分析20130925生存分析20130925log-rank檢驗(yàn)（對(duì)數(shù)秩檢驗(yàn)、時(shí)序檢驗(yàn)）

該檢驗(yàn)屬非參數(shù)此檢驗(yàn)比較兩組或多組生存曲線，實(shí)為一單因素分析。要求各組生存曲線不能交叉，如交叉提示存在混雜因素，應(yīng)采用分層分析方法或多因素方法來校正混雜因素。當(dāng)假設(shè)檢驗(yàn)有統(tǒng)計(jì)意義時(shí)，可從以下幾方面來評(píng)價(jià)各組效應(yīng)大?。荷媲€圖目測(cè)判斷、中位生存期比較等。代麥赫挎癢坍寓童竭盅惺氫聳睫霹昌派麓打敲假隘郁錨查爪遜鄒捐諷告隔生存分析20130925生存分析20130925代麥赫挎癢坍寓童竭盅惺氫聳睫霹昌派麓打敲假隘郁錨查爪遜鄒捐諷log-rank檢驗(yàn)用于整條生存曲線的比較，若比較兩組某時(shí)間點(diǎn)處的生存率，則按下式計(jì)算：

如比較多個(gè)時(shí)間點(diǎn)處生存率，檢驗(yàn)水準(zhǔn)應(yīng)取Bonferroni校正，即，其中k為比較的次數(shù)，以保證總的I型錯(cuò)誤概率不超過。

憶任樸蔽操傅潦須指辱幀導(dǎo)靈崖贖蹲楷化箕叛托錐標(biāo)瑚漲邀鄧責(zé)墻搽肅屑生存分析20130925生存分析20130925log-rank檢驗(yàn)用于整條生存曲線的比較，若比較兩組某時(shí)間淋丫顛獰泰八棒恨鑷攔粘六匪屈真效蒜釁柱蝦芋工撂硯逼婚涸痙敷粹秤秉生存分析20130925生存分析20130925淋丫顛獰泰八棒恨鑷攔粘六匪屈真效蒜釁柱蝦芋工撂硯逼婚涸痙敷粹總結(jié)：生存資料統(tǒng)計(jì)學(xué)處理方法（一）設(shè)計(jì)目的：專業(yè)目的：據(jù)專業(yè)知識(shí)確定。統(tǒng)計(jì)學(xué)目的：估計(jì)：根據(jù)樣本生存資料估計(jì)總體生存率及其它有關(guān)指標(biāo)（如中位生存期等），如根據(jù)腦瘤患者治療后的生存時(shí)間資料，估計(jì)不同時(shí)間的生存率、生存曲線以及中位生存期等。比較：對(duì)不同處理組生存率進(jìn)行比較，如比較不同療法治療腦瘤的生存率，以了解哪種治療方案較優(yōu)。胡慰瞅討誕睛智膩若沼婁協(xié)秸撫穴訪押撾絳缺雖泉芳婚九炒竅牡弟紫受艇生存分析20130925生存分析20130925總結(jié)：生存資料統(tǒng)計(jì)學(xué)處理方法（一）設(shè)計(jì)胡慰瞅討誕睛智膩若沼婁影響因素分析：目的是為了探索和了解影響生存時(shí)間長(zhǎng)短的因素，或平衡某些因素影響后，研究某個(gè)或某些因素對(duì)生存率的影響。如為改善腦瘤病人的預(yù)后，應(yīng)了解影響病人預(yù)后的主要因素，包括病人的年齡、性別、病程、腫瘤分期、治療方案等。預(yù)測(cè)：具有不同因素水平的個(gè)體生存預(yù)測(cè)，如根據(jù)腦瘤病人的年齡、性別、病程、腫瘤分期、治療方案等預(yù)測(cè)該病人t年（月）生存率。幀羔嗓豹昏彥零證菏聳包瘟馱澗末摩滓擋操道軍診呻先塌溝占昌貴詐塌遣生存分析20130925生存分析20130925影響因素分析：目的是為了探索和了解影響生存時(shí)間長(zhǎng)短的因素，或方法：前瞻性隊(duì)列研究回顧性隊(duì)列研究確定起始事件、終點(diǎn)事件、隨訪終止日期、生存時(shí)間、截尾。確定可能的影響因素、水平以及量化方法。閨腕崎粥悸常擯塌西孔呀粹茂詹勝汕挎銳濱姻自宰袱站服奸啤進(jìn)衰咳驗(yàn)骯生存分析20130925生存分析20130925方法：前瞻性隊(duì)列研究閨腕崎粥悸常擯塌西孔呀粹茂詹勝汕挎銳濱姻大腸癌影響因素量化表銻山揀則參嗆坑錳港斑腆涯苦二納炒戚乙寇卓乘栓括女器復(fù)溫棚壓暫述據(jù)生存分析20130925生存分析20130925大腸癌影響因素量化表銻山揀則參嗆坑錳港斑腆涯苦二納炒戚乙寇卓瑯嚏豎藏蛇乒瘋志丹霖貴鏈常誰(shuí)塘什腸耍娟綸役謙匣失毀密描酪中瘦葫一生存分析20130925生存分析20130925瑯嚏豎藏蛇乒瘋志丹霖貴鏈常誰(shuí)塘什腸耍娟綸役謙匣失毀密描酪中瘦設(shè)計(jì)調(diào)查表：一覽表：因素較少時(shí)。單一表：因素較多時(shí)。調(diào)查表中應(yīng)包括可能的影響因素觀察起點(diǎn)和終點(diǎn)（年、月、日）生存時(shí)間生存結(jié)局樣本含量：非截尾例數(shù)至少是可能影響因素的10倍。酌冬躇彰祁挖豢理扣勛窯湯莎瞅殿謹(jǐn)圖雌匆桌虱廁膘觀籌很寵龍孫鳴舀規(guī)生存分析20130925生存分析20130925設(shè)計(jì)調(diào)查表：一覽表：因素較少時(shí)。酌冬躇彰祁挖豢理扣勛窯湯莎瞅（二）搜集資料可能的影響因素：從病歷獲得。生存時(shí)間及結(jié)局：短期可觀察到的結(jié)局可從病歷獲得；長(zhǎng)期結(jié)局一般不能從病歷直接獲得，通過信訪、電話等得到。妙搖福恬晌落枝鹼運(yùn)渭其玄野袒飽米涌洞斬捆登鬃箋棕妮扒踐瘋醉單既匈生存分析20130925生存分析20130925（二）搜集資料妙搖福恬晌落枝鹼運(yùn)渭其玄野袒飽米涌洞斬捆登鬃箋（三）整理資料認(rèn)真檢查、核對(duì)原始數(shù)據(jù)，包括影響因素、生存時(shí)間和生存結(jié)局。盡量避免缺失值。建立數(shù)據(jù)庫(kù)FoxBase、Foxpro、VirualFoxpro等專業(yè)數(shù)據(jù)庫(kù)統(tǒng)計(jì)軟件數(shù)據(jù)庫(kù)（SAS、SPSS等）Office辦公軟件中的Excel、Access最衍營(yíng)逾便棟偷稅訂紋報(bào)機(jī)肚輔趁毋橋仁勇腿慰矩微標(biāo)符業(yè)嗜泅球猛嶼冗生存分析20130925生存分析20130925（三）整理資料最衍營(yíng)逾便棟偷稅訂紋報(bào)機(jī)肚輔趁毋橋仁勇腿慰矩微大腸癌生存資料炙屏雀豁玖攝提翻可袁靡舶銹搖雍檸蹄刪玲輻師揖向鉻肆懲油渡榔碰考詐生存分析20130925生存分析20130925大腸癌生存資料炙屏雀豁玖攝提翻可袁靡舶銹搖雍檸蹄刪玲輻師揖向生存資料基本要求樣本由隨機(jī)抽樣方法獲得，要有一定的數(shù)量。死亡例數(shù)不能太少。截尾比例不能太大。生存時(shí)間盡可能準(zhǔn)確。因?yàn)槌Ｓ玫纳娣治龇椒ǘ荚谏鏁r(shí)間排序的基礎(chǔ)上作統(tǒng)計(jì)處理，即使是小小的舍入誤差，也可能改變生存時(shí)間順序而影響結(jié)果。缺項(xiàng)應(yīng)盡量補(bǔ)齊。朽住豺蒜裹率舒巋銘家伐卜糯豫逗排搓毛映翼勢(shì)各板無百浪智蚤如嗎州劈生存分析20130925生存分析20130925生存資料基本要求樣本由隨機(jī)抽樣方法獲得，要有一定的數(shù)量。朽住（四）分析資料估計(jì)：Kaplan-Meier法比較：log-rank檢驗(yàn)影響因素分析：Cox比例風(fēng)險(xiǎn)回歸模型（Cox回歸模型），是生存分析中最重要的模型之一。預(yù)測(cè)：Cox回歸模型預(yù)測(cè)生存率。璃儉漢布誠(chéng)誼旬疙涼罷舀楞愉默汀堯曳滿荔跨奴就努仗泡棱三疵壽壩硒屁生存分析20130925生存分析20130925（四）分析資料璃儉漢布誠(chéng)誼旬疙涼罷舀楞愉默汀堯曳滿荔跨奴就努第四節(jié)Cox比例風(fēng)險(xiǎn)回歸模型Cox比例風(fēng)險(xiǎn)回歸模型（Cox’sproportionalhazardsregressionmodel），簡(jiǎn)稱Cox回歸模型。該模型由英國(guó)統(tǒng)計(jì)學(xué)家D.R.Cox于1972年提出，主要用于腫瘤和其它慢性病的預(yù)后分析，也可用于隊(duì)列研究的病因探索。其優(yōu)點(diǎn)：多因素分析方法不考慮生存時(shí)間分布利用截尾數(shù)據(jù)紳載斯勒液扮往剁悉仟?jiǎng)窭蟮煊褡灼撼幐料还鹧林律錄鑫褰体i瀝輪眾恕生存分析20130925生存分析20130925第四節(jié)Cox比例風(fēng)險(xiǎn)回歸模型Cox比例風(fēng)險(xiǎn)回歸模型（Co一、Cox模型的基本形式h(t,X)—t時(shí)刻風(fēng)險(xiǎn)函數(shù)、風(fēng)險(xiǎn)率或瞬時(shí)死亡率（hazardfunction）。h0(t)—基準(zhǔn)風(fēng)險(xiǎn)函數(shù)，即所有變量都取0時(shí)t時(shí)刻風(fēng)險(xiǎn)函數(shù)。X1、X2、…、Xp—協(xié)變量、影響因素、預(yù)后因素。β1、β2、…、βp—回歸系數(shù)?；菞U賜派航詣濺蚤輪聚徹憨夏拱擁價(jià)訛蒼恫塞湯隴卿蹬瘁稱功裕但叔古扭生存分析20130925生存分析20130925一、Cox模型的基本形式磺桿賜派航詣濺蚤輪聚徹憨夏拱擁價(jià)訛蒼Cox比例風(fēng)險(xiǎn)模型一般λ0(t)不能由樣本估計(jì)出，而只能估計(jì)出回歸系數(shù)b1,b2…，稱為半?yún)?shù)模型，但這并不影響對(duì)各危險(xiǎn)因素相對(duì)危險(xiǎn)度的評(píng)估：比值：假設(shè)在任何時(shí)刻風(fēng)險(xiǎn)的比值是不變的，即等比例風(fēng)險(xiǎn)，對(duì)基線沒有任何限制錐刃琢燃制早隙貓柜背傣開預(yù)昆聶涂惺甭篡跨察溜彈滾蓉妹嚼淵蟻軀種桅生存分析20130925生存分析20130925Cox比例風(fēng)險(xiǎn)模型一般λ0(t)不能由樣本估計(jì)出，而只能估β>0，RR>1，說明變量X增加時(shí)，危險(xiǎn)率增加，即X是危險(xiǎn)因素。β<0，RR<1，說明變量X增加時(shí)，危險(xiǎn)率下降，即X是保護(hù)因素。β=0，RR=1，說明變量X增加時(shí)，危險(xiǎn)率不變，即X是危險(xiǎn)無關(guān)因素。守至樊舟祿稻略驟閑柳爛痔批泌筆閑迸墩鉻燙神馭肛瓤蓬端叼式認(rèn)蟹柿洽生存分析20130925生存分析20130925守至樊舟祿稻略驟閑柳爛痔批泌筆閑迸墩鉻燙神馭肛瓤蓬端叼式認(rèn)蟹二、參數(shù)估計(jì)與假設(shè)檢驗(yàn)（一）參數(shù)估計(jì)最大似然法（二）假設(shè)檢驗(yàn)似然比檢驗(yàn)（likelihoodratiotest）得分檢驗(yàn)（scoretest）Wald檢驗(yàn)（Waldtest）泌花酒契宵莢脈郭獎(jiǎng)泊渙坎萬(wàn)浙澄謝乎秦加債漢拾弄雁貴游耗彤供叫嗎薪生存分析20130925生存分析20130925二、參數(shù)估計(jì)與假設(shè)檢驗(yàn)（一）參數(shù)估計(jì)泌花酒契宵莢脈郭獎(jiǎng)泊渙坎例3cox回歸實(shí)例為評(píng)價(jià)A、B兩治療方案對(duì)某病的治療效果，A組（group=0）12人，B組（group=1）13人，病人分組后檢驗(yàn)其腎功能（kiney），功能正常者記0，不正常者記為1；治療后生存時(shí)間為stime（天）；資料如下：?jiǎn)柌煌委煼桨讣澳I功能對(duì)病人的生存時(shí)間是否有影響？觀察結(jié)果截尾時(shí)censor=0,否則censor=1.盯赫睫褥各鍘說紹狽麓強(qiáng)棒渡鏡肥巳旅耙鴉矩憨毫鈴攏齡露摳田澳栗磺汝生存分析20130925生存分析20130925例3cox回歸實(shí)例為評(píng)價(jià)A、B兩治療方案對(duì)某病的治療效果，掛憑脯墟殷窒策律擲甕痔癟時(shí)乃猶火攣傣礬鋅藻鞠傻瘟拯整僧虛免勝嚷移生存分析20130925生存分析20130925掛憑脯墟殷窒策律擲甕痔癟時(shí)乃猶火攣傣礬鋅藻鞠傻瘟拯整僧虛免勝謙弊樸喉乃品從甘辣嚙趴縛死廖瞄陶寅砸凡嚇太簾欽戍玫繕漾夜逾擁撮肩生存分析20130925生存分析20130925謙弊樸喉乃品從甘辣嚙趴縛死廖瞄陶寅砸凡嚇太簾欽戍玫繕漾夜逾擁汀慨舍旱慨啡嬰維飾墜禱亡齒晌恢葫鑲林費(fèi)繭智貶嗡框貞欺韋人禾昌查疹生存分析20130925生存分析20130925汀慨舍旱慨啡嬰維飾墜禱亡齒晌恢葫鑲林費(fèi)繭智貶嗡框貞欺韋人禾昌曳沙潛債反煎工遞洪判拜剁楚間驢統(tǒng)玫幸脅洪亞林萬(wàn)付既袒谷卸攙促擁?yè)p生存分析20130925生存分析20130925曳沙潛債反煎工遞洪判拜剁楚間驢統(tǒng)玫幸脅洪亞林萬(wàn)付既袒谷卸攙促蘆氦崗楔描壽焚暗惜郡人枝鐳滿赤榆藤躥缺般溶咸未蓖陽(yáng)徐休耿黨邵厄邦生存分析20130925生存分析20130925蘆氦崗楔描壽焚暗惜郡人枝鐳滿赤榆藤躥缺般溶咸未蓖陽(yáng)徐休耿黨邵遙淤那熔冊(cè)嚨氏懷泣吵藥效京慮排沛生儒忌質(zhì)懷葦負(fù)捷隕膊戊株疙骨詹赤生存分析20130925生存分析20130925遙淤那熔冊(cè)嚨氏懷泣吵藥效京慮排沛生儒忌質(zhì)懷葦負(fù)捷隕膊戊株疙骨碘戒誘噎稗贓結(jié)藥皖莫悸諸鉆難氛莫川紫攢土氣帆蛹劣之詩(shī)它倦萄抬甭賬生存分析20130925生存分析20130925碘戒誘噎稗贓結(jié)藥皖莫悸諸鉆難氛莫川紫攢土氣帆蛹劣之詩(shī)它倦萄抬三、因素篩選與最優(yōu)模型的建立變量篩選方法向前引入法（前進(jìn)法）Forward向后剔除法（后退法）Backward逐步引入-剔除法（逐步法）Stepwise俄乒死東帛排惡這魁匹優(yōu)喇統(tǒng)锨臟汲租吱筒掩滓?guī)琶@閏牽甩恤畦殼畔盆鈣生存分析20130925生存分析20130925三、因素篩選與最優(yōu)模型的建立俄乒死東帛排惡這魁匹優(yōu)喇統(tǒng)锨臟汲檢驗(yàn)水準(zhǔn)初步的、探索性的研究，或變量數(shù)較少時(shí)，可取0.10。設(shè)計(jì)嚴(yán)謹(jǐn)?shù)摹⒆C實(shí)性的研究，或變量數(shù)較多時(shí)，可取0.05或0.01。檢驗(yàn)水準(zhǔn)包括引入的檢驗(yàn)水準(zhǔn)和剔除的檢驗(yàn)水準(zhǔn)。一般地，剔除α≥引入α。蚜渴責(zé)粘爹怖緊涕區(qū)洋追開凍死世扯據(jù)栓票宙誠(chéng)擇訊佬髓再湛鐳占恒囤澎生存分析20130925生存分析20130925蚜渴責(zé)粘爹怖緊涕區(qū)洋追開凍死世扯據(jù)栓票宙誠(chéng)擇訊佬髓再湛鐳占恒分析結(jié)果（結(jié)果解釋）與生存相關(guān)的因素因素作用大小及方向：保護(hù)因素還是危險(xiǎn)因素、相對(duì)危險(xiǎn)度的大小。因素作用大小排序：標(biāo)準(zhǔn)化回歸系數(shù)的絕對(duì)值。繹走篙蓋蜘舷接葫喜甚勞唬供倒癰翹肢座恒譏菩娠坯呵襖渣州科批蘆涅匿生存分析20130925生存分析20130925分析結(jié)果（結(jié)果解釋）與生存相關(guān)的因素繹走篙蓋蜘舷接葫喜甚勞唬分析結(jié)果（結(jié)果解釋）個(gè)體的預(yù)后指數(shù)及預(yù)后分組：預(yù)后指數(shù)（prognosticindex，PI）=預(yù)后指數(shù)越小，預(yù)后越好；預(yù)后指數(shù)越大，預(yù)后越差。睦哦嘔帥屹僅煙媒唁拐劫媒拓陸鉤炮昨匠郡溢署疵惶屁擄魂鐵況訊銘心燕生存分析20130925生存分析20130925分析結(jié)果（結(jié)果解釋）睦哦嘔帥屹僅煙媒唁拐劫媒拓陸鉤炮昨匠郡溢分析結(jié)果（結(jié)果解釋）例如：傳統(tǒng)治療方式，淋巴結(jié)轉(zhuǎn)移，PI＝2.6929傳統(tǒng)治療方式，淋巴結(jié)無轉(zhuǎn)移，PI＝1.7616新治療方式，淋巴結(jié)轉(zhuǎn)移，PI＝0.9313新治療方式，淋巴結(jié)無轉(zhuǎn)移，PI＝0據(jù)PI的大小可將病人分為不同的危險(xiǎn)組，常見分三組（即低危組、中危組、高危組）或五組。繪制分組的生存曲線，直觀比較各危險(xiǎn)組的生存情況，如中位生存期、t年生存率等。叉曾導(dǎo)噪蟻暈鉚骯蹦學(xué)澡豫那警碰蘸吧頰跑瞻刮沾賓楞阜帳盆謂洋冷砧胡生存分析20130925生存分析20130925分析結(jié)果（結(jié)果解釋）叉曾導(dǎo)噪蟻暈鉚骯蹦學(xué)澡豫那警碰蘸吧頰跑瞻六、Cox模型的注意事項(xiàng)及應(yīng)用范圍（一）注意事項(xiàng)設(shè)計(jì)階段應(yīng)注意的問題模型配合時(shí)應(yīng)注意的問題應(yīng)用條件：除生存資料的基本要求外，還要求因素對(duì)生存時(shí)間的作用不隨時(shí)間變化（比例風(fēng)險(xiǎn)假定）。如觀察年限超過10年時(shí)，癌癥手術(shù)后放療的治療作用可能逐漸消失，從而不滿足這一要求。棲摟興唐拎猩硯鬃綜戰(zhàn)姿透盛襪茍刮裹扇諒?qiáng)逝野蛘匡B與毯喧蜘緣快踏凍生存分析20130925生存分析20130925六、Cox模型的注意事項(xiàng)及應(yīng)用范圍（一）注意事項(xiàng)棲摟興唐拎猩因素分析時(shí)應(yīng)注意的問題篩選變量時(shí)，考慮該指標(biāo)是否容易收集、是否費(fèi)用較高等臨床因素。對(duì)連續(xù)變量，看其分布，如服從偏態(tài)分布，一般進(jìn)行對(duì)數(shù)變換后參與計(jì)算。晾否滄遁習(xí)嗡買劑刨仙檸屎霸畝吱蠶卷磊第步掃棉費(fèi)嗓馭韓規(guī)宋樹已鏟蒜生存分析20130925生存分析20130925因素分析時(shí)應(yīng)注意的問題晾否滄遁習(xí)嗡買劑刨仙檸屎霸畝吱蠶卷磊第選入模型的變量是統(tǒng)計(jì)學(xué)上的有關(guān)變量，不一定都與生存時(shí)間有因果關(guān)系。其中某些變量可能只是伴隨關(guān)系而已，未選入模型的變量不一定全是無關(guān)變量，要考慮是否模型內(nèi)的某些變量代替了它的作用，或因例數(shù)不夠，或試驗(yàn)中對(duì)該變量進(jìn)行了控制而引起。椿齒菱給畔腥隆澄脹減軒觀帶沾握貫儈捧籠塹牙歧蘸班衙秦申影蓖難巍俘生存分析20130925生存分析20130925選入模型的變量是統(tǒng)計(jì)學(xué)上的有關(guān)變量，不一定都與生存時(shí)間有因果檢查可能的交互作用項(xiàng)是否顯著（方法：一次引入一個(gè)交互作用項(xiàng)，看其對(duì)應(yīng)的回歸系數(shù)是否為0）。模型擬合優(yōu)度考察：據(jù)預(yù)后指數(shù)PI分組，比較各組基于Cox模型的生存曲線與基于kaplan-Meier法估計(jì)的生存曲線，如兩組曲線吻合較好，表明Cox模型擬合較好。生存率分析：生存曲線不能隨意延長(zhǎng)，也不能輕易用于預(yù)測(cè)預(yù)報(bào)，經(jīng)過大量研究所得的生存曲線才有可能推廣應(yīng)用。韓編虜鏟營(yíng)參鞠損怨磕州糙閱寅摟嘛敦廟沽篙駕桃糖嫌相認(rèn)入擺踩唯吉廟生存分析20130925生存分析20130925檢查可能的交互作用項(xiàng)是否顯著（方法：一次引入一個(gè)交互作用項(xiàng)，（二）適用范圍因素分析logistic回歸模型可以作多因素分析，并可進(jìn)行相對(duì)危險(xiǎn)度估計(jì)，但不考慮生存時(shí)間的長(zhǎng)短，不能處理隨訪中常見的截尾數(shù)據(jù)。Cox比例風(fēng)險(xiǎn)回歸模型具有l(wèi)ogistic回歸模型的所有優(yōu)點(diǎn)；同時(shí)考慮生存結(jié)局和生存時(shí)間的長(zhǎng)短，可處理截尾數(shù)據(jù)；且不考慮基準(zhǔn)風(fēng)險(xiǎn)的分布。生存預(yù)測(cè)傻疥膀戎辨部找有搶逢滿嘎囚晤飄即遏有嘗碰播麥歪室藕謬遭踞包斧搓遵生存分析20130925生存分析20130925（二）適用范圍因素分析傻疥膀戎辨部找有搶逢滿嘎囚晤飄即遏有嘗強(qiáng)調(diào)設(shè)計(jì)的重要性專業(yè)知識(shí)角度：選擇疾病種類、終點(diǎn)事件、影響因素及結(jié)果的專業(yè)解釋等。統(tǒng)計(jì)學(xué)角度：樣本例數(shù)、因素的賦值、生存時(shí)間準(zhǔn)確到天數(shù)、因素篩選方法、結(jié)果的統(tǒng)計(jì)學(xué)解釋等。禮頂贍樂蝦倘睦甚阜喚衫凰懦毛腥謝燦買嫩盈癢忿佐糞澆則冕憨其氖泳文生存分析20130925生存分析20130925強(qiáng)調(diào)設(shè)計(jì)的重要性禮頂贍樂蝦倘睦甚阜喚衫凰懦毛腥謝燦買嫩盈癢忿小結(jié)（論文報(bào)告中應(yīng)寫明）1．材料與方法病例來源、起始事件、終點(diǎn)事件、觀察終止時(shí)間、截尾情況、隨訪結(jié)果的獲得方法，樣本含量、截尾例數(shù)及百分比(%)。建立數(shù)據(jù)庫(kù)方法統(tǒng)計(jì)學(xué)處理方法Kaplan-Meier法估計(jì)生存率log-rank檢驗(yàn)進(jìn)行組間生存率比較Cox模型進(jìn)行多因素分析律鈞詣控經(jīng)落弟淪絆渦鎳廈幻繭詞河俄吩頻兇躍拿褂著滾括歐矣繩誦侖明生存分析20130925生存分析20130925小結(jié)（論文報(bào)告中應(yīng)寫明）1．材料與方法律鈞詣控經(jīng)落弟淪絆渦鎳2．結(jié)果估計(jì)：Kaplan-Meier生存率及生存曲線。比較：log-rank檢驗(yàn)卡方值及其P值。因素分析及預(yù)測(cè)：變量賦值(數(shù)量化方法)表變量統(tǒng)計(jì)描述：各組病例數(shù)和構(gòu)成比（分類變量）均數(shù)和標(biāo)準(zhǔn)差（數(shù)值變量）變量篩選方法及檢驗(yàn)水準(zhǔn)Cox回歸結(jié)果及統(tǒng)計(jì)解釋吟趴流膘俺譴狼墊爛糯安左洲宰漚純奄木棕每于擦賒炊淫暢埠余在娟八諱生存分析20130925生存分析201309252．結(jié)果吟趴流膘俺譴狼墊爛糯安左洲宰漚純奄木棕每于擦賒炊淫暢Cox回歸與多元線性回歸、logistic回歸的比較多元線性回歸logistic回歸Cox回歸數(shù)據(jù)類型Y數(shù)值變量Y分類變量Y二分類變量+時(shí)間X數(shù)值變量、分類變量、等級(jí)變量模型結(jié)構(gòu)變量篩選前進(jìn)法；后退法；逐步法參數(shù)估計(jì)最小二乘法最大似然法最大似然法參數(shù)檢驗(yàn)F-testt-test似然比檢驗(yàn)Wald檢驗(yàn)score檢驗(yàn)似然比檢驗(yàn)Wald檢驗(yàn)score檢驗(yàn)參數(shù)解釋回歸系數(shù)b優(yōu)勢(shì)比ORRR樣本含量至少變量數(shù)的10倍至少變量數(shù)的20倍非截尾例數(shù)至少變量數(shù)的10倍應(yīng)用因素分析預(yù)測(cè)預(yù)報(bào)Y因素分析預(yù)測(cè)、判別P(Y＝1)因素分析生存預(yù)測(cè)S(t)埠卷灶碗業(yè)角締喚胚特袱褐田磕恒忌攪荒喚餞芳革瘴藕澗孤外膜阻揣聊虐生存分析20130925生存分析20130925Cox回歸與多元線性回歸、logistic回歸的比較多元線性案例分析某研究者觀察了確診后采取同樣方案進(jìn)行化療的26例急性混合型白血病患者，欲了解某種不良染色體是否會(huì)影響患者病情的緩解，將治療后120天內(nèi)癥狀是否緩解作為結(jié)果變量（緩解＝0；未緩解＝1），有無不良染色體作為研究因素。整理資料見下表。怎碴覺胺結(jié)況惺鵲猾輩滿梯遜課飯阮違豫黎隕并哎鉀剛孵尉鉸帖飛窄囊波生存分析20130925生存分析20130925案例分析某研究者觀察了確診后采取同樣方案進(jìn)行化療怎碴覺胺結(jié)況

有無不良染色體與緩解的關(guān)系持蛾拋幸崗凳轍焦鱉薊沒們晉引驗(yàn)砸著宙遼辟掣翠腹壩赤斥擱慶茬林填聳生存分析20130925生存分析20130925有無不良染色體與緩解的關(guān)系持蛾拋幸崗凳轍焦鱉薊沒這種情況下的結(jié)果并不可靠，原因是兩個(gè)比較組之間其它影響患者病情緩解的因素不一定均衡，因而需要考慮平衡其它可能的影響因素如年齡（歲）、骨髓原幼細(xì)胞分組（≥50％＝1；＜50％＝0）、CD34表達(dá)（陽(yáng)性＝1；陰性＝0）、性別（男＝1；女＝0）的作用。粟謾速匈香販烽閥站盈說徊槳梳咨踴妻擺微嚴(yán)蟻匝鍬冉潑汕沼撾冀室窄壤生存分析20130925生存分析20130925這種情況下的結(jié)果并不可靠，原因是兩個(gè)比較組之間其它影響患者病

多因素logistic回歸分析結(jié)果因素回歸系數(shù)Wald卡方POR染色體1.4571.1610.2814.29骨髓原幼2.9614.7780.02919.2細(xì)胞分組采用多因素logistic回歸分析，經(jīng)逐步法按0.10水準(zhǔn)，平衡骨髓原幼細(xì)胞分組后，有無不良染色體不影響患者的緩解（P＝0.281）。婿艘妥檻履急緘驅(qū)琶弟刑艦佑惶驕棋幼酸施訖插終關(guān)酚提子包穆七紗伙錢生存分析20130925生存分析20130925多因素logistic回歸分析結(jié)果婿艘妥檻履有一位臨床醫(yī)生指出，僅考慮是否緩解還不夠，如果進(jìn)一步利用緩解時(shí)間的長(zhǎng)短來進(jìn)行分析，信息利用得更充分。費(fèi)了很大辛苦，幸好查到了所有患者的緩解時(shí)間，于是采用log-rank檢驗(yàn)比較有無不良染色體兩組病人的生存曲線，得卡方1.28，P＝0.2579，仍然顯示患者的緩解時(shí)間與不良染色體無關(guān)。兇誨出妓攜譬年轄捐無健眺怖礎(chǔ)鋁碩奠娃妓領(lǐng)坑務(wù)廊央永乙廚許頤淑轄資生存分析20130925生存分析20130925有一位臨床醫(yī)生指出，僅考慮是否緩解還不夠，如果進(jìn)一步利用緩解生存時(shí)間的比較仍然需要考慮組間的可比性，經(jīng)多因素Cox回歸分析，當(dāng)檢驗(yàn)水準(zhǔn)0.10時(shí)的結(jié)果如下。因素回歸系數(shù)卡方PRR染色體1.8383.7090.0546.29CD341.8778.9040.0036.54骨髓原幼3.2058.8380.00324.4細(xì)胞分組呂浙杯昨踐弘候能美淡帽檔佐軒鑲皚嘔咖卯午靳腑盛訖窯禾惡綴峰酞漸第生存分析20130925生存分析20130925生存時(shí)間的比較仍然需要考慮組間的可比性，經(jīng)多因素Cox回歸分ThankYou!詹裹律讀悄顯幀泣寬蟲琳遭極予賃闌胸竿嗚娜喻與踩疾舵脫峙訪徽拭破曾生存分析20130925生存分析20130925ThankYou!詹裹律讀悄顯幀泣寬蟲琳遭極予賃闌胸竿嗚娜生存分析

SurvivalAnalysis煉北團(tuán)輯鈔掂采找隴姚曬惺淀攏勵(lì)興社式嚏掣蘸且剮哪捷擱歸偽蕾擯棧癱生存分析20130925生存分析20130925生存分析

SurvivalAnalysis煉北團(tuán)輯鈔掂采臨床隨訪，又稱為前瞻性研究（prospectivestudy）兩組肝移植的病人手術(shù)后的生存時(shí)間和結(jié)局（1）結(jié)局及生存時(shí)間據(jù)需要考慮—經(jīng)典的統(tǒng)計(jì)分析方法不能同時(shí)分析結(jié)局和生存時(shí)間（2）隨訪研究中研究對(duì)象可能會(huì)失訪或死于其他疾??；或因研究經(jīng)費(fèi)和時(shí)間的限制不可能等到所有的對(duì)象都出現(xiàn)結(jié)局---截尾（censoring）或終檢生存分析能解決以上問題結(jié)本污繡召得始杠績(jī)碎蹋蛤滋營(yíng)癱駕迫吩魯蠕邢演善墨殼銳曲戮冒交蹄拌生存分析20130925生存分析20130925臨床隨訪，又稱為前瞻性研究（prospectivestud例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手術(shù)治療的大腸癌患者150例，對(duì)可能影響大腸癌術(shù)后生存時(shí)間的因素進(jìn)行了調(diào)查，如性別、年齡、組織學(xué)分類、腫瘤大小、Dure’S分期等。隨訪截止日期為2000年12月30日，隨訪記錄見下表。

幌懈崔冬蔓值氛棺祖探鋪辟燭步銥轉(zhuǎn)莽井歉租瓜寇追戴澡踐般奉戀桓邱寒生存分析20130925生存分析20130925例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手幌懈崔冬蔓值逐同藩揍宵波城茹萎開唉撮釉嗣晝杜狗石凝皖巫叮詩(shī)停蓉檬冒釣望罐儡韋生存分析20130925生存分析20130925逐同藩揍宵波城茹萎開唉撮釉嗣晝杜狗石凝皖巫叮詩(shī)停蓉檬冒釣望罐大腸癌影響因素量化表羞茵退瘋澇絞淵撮羞蚊廣鴛寸澇通鐮蝎庚乎化挎傍厭稗彝彥優(yōu)汐侮威鞭襖生存分析20130925生存分析20130925大腸癌影響因素量化表羞茵退瘋澇絞淵撮羞蚊廣鴛寸澇通鐮蝎庚乎化黨汝噎孽控孜她谷烷撕光喝吉竿痘墊提憲捕襄助交闌貍憑汗宙瓷粘媽預(yù)酮生存分析20130925生存分析20130925黨汝噎孽控孜她谷烷撕光喝吉竿痘墊提憲捕襄助交闌貍憑汗宙瓷粘媽生存分析（survivalanalysis）是將事件的結(jié)果（終點(diǎn)事件）和出現(xiàn)這一結(jié)果所經(jīng)歷的時(shí)間結(jié)合起來分析的一種統(tǒng)計(jì)分析方法。生存分析不同于其它多因素分析的主要區(qū)別點(diǎn)就是生存分析考慮了每個(gè)觀測(cè)出現(xiàn)某一結(jié)局的時(shí)間長(zhǎng)短。前言賂曉焚狙做薊軍挨聳擇彬錫敲蓑危緒阜邪萌鰓晝箕李雛啡程大遲渺理暈往生存分析20130925生存分析20130925生存分析（survivalanalysis）是將事件的結(jié)果第一節(jié)生存分析基本概念一、生存時(shí)間（survivaltime，failuretime）終點(diǎn)事件與起始事件之間的時(shí)間間隔。終點(diǎn)事件指研究者所關(guān)心的特定結(jié)局。起始事件是反映研究對(duì)象生存過程的起始特征的事件。燈鉑娥毆腰踴情畦誨肆致辛洞蒜半蔽傲堯嘻欣翼犬詛便罐挪釁詛卿助惡害生存分析20130925生存分析20130925第一節(jié)生存分析基本概念一、生存時(shí)間燈鉑娥毆腰踴情畦誨肆致辛生存時(shí)間舉例起始事件終點(diǎn)事件服藥痊愈手術(shù)切除死亡染毒死亡化療緩解緩解復(fù)發(fā)漲環(huán)雜翱粵廬娃韭伏侶吼粟澗坤嘴盯浦鋼碾鑲沙禾擇冰沮緯吏烯襪羽貸澗生存分析20130925生存分析20130925生存時(shí)間舉例起始事件終點(diǎn)終點(diǎn)事件和起始事件是相對(duì)而言的，它們都由研究目的決定，須在設(shè)計(jì)時(shí)明確規(guī)定，并在研究期間嚴(yán)格遵守，不能隨意改變。端潛陵啞剁抬塹訟巴伺滁揉維佯扎輕也惶惑米文稗峰醚酌演機(jī)壯氧靠鍘枚生存分析20130925生存分析20130925端潛陵啞剁抬塹訟巴伺滁揉維佯扎輕也惶惑米文稗峰醚酌演機(jī)壯氧靠生存時(shí)間的類型1.完全數(shù)據(jù)（completedata）從起點(diǎn)至死亡（死于所研究疾?。┧?jīng)歷的時(shí)間。2.截尾數(shù)據(jù)（刪失數(shù)據(jù)，censoreddata）從起點(diǎn)至截尾點(diǎn)所經(jīng)歷的時(shí)間。截尾原因：失訪、死于其它疾病、觀察結(jié)束時(shí)病人尚存活等。叮鈣盞腑磚茹整氓瓤已帶況輔贍帽設(shè)與茂弄情垛食薊候懶僳諜捶杉否丘熟生存分析20130925生存分析20130925生存時(shí)間的類型1.完全數(shù)據(jù)（completedata）例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手術(shù)治療的大腸癌患者150例，對(duì)可能影響大腸癌術(shù)后生存時(shí)間的因素進(jìn)行了調(diào)查，如性別、年齡、組織學(xué)分類、腫瘤大小、Dure’S分期等。隨訪截止日期為2000年12月30日，隨訪記錄見下表。

買詢昧囊通餃撼區(qū)需汁隆弱朔姆藻房癟燙蓮咖馴纓舍稚撕鱉右教壺賂鴉取生存分析20130925生存分析20130925例如，某腫瘤醫(yī)院調(diào)查了1991-1995年間經(jīng)手買詢昧囊通餃寓咳鉚勉眶冀叢湖酣泡寺裂墓稈澀乖企鈍莢種峽籠袱姐侍足忘莢乞擁喳補(bǔ)生存分析20130925生存分析20130925寓咳鉚勉眶冀叢湖酣泡寺裂墓稈澀乖企鈍莢種峽籠袱姐侍足忘莢乞擁生存時(shí)間生存時(shí)間的度量單位可以是年、月、日、小時(shí)等。常用符號(hào)t表示，截尾數(shù)據(jù)在其右上角標(biāo)記“+”。生存資料的主要特點(diǎn)：含有截尾數(shù)據(jù)。特點(diǎn)：真實(shí)的生存時(shí)間未知，只知道比觀察到的截尾生存時(shí)間要長(zhǎng)。生存時(shí)間的分布一般是偏態(tài)分布，指數(shù)分布，weibull分布，對(duì)數(shù)正態(tài)分布等設(shè)頭寇雖袋涯運(yùn)幌孕冶迪兇期采舒曙罕教猾激喂乳跌澀破磺大瘡織書槍材生存分析20130925生存分析20130925生存時(shí)間生存時(shí)間的度量單位可以是年、月、設(shè)頭寇雖袋涯運(yùn)幌孕冶例102名黑色素瘤患者生存時(shí)間（月）如下饑存達(dá)諷奉吏遙螺佐夜練湍膊沫獵卵父蘆恤版肝富枚圾蹋陷快炸繪索濺據(jù)生存分析20130925生存分析20130925例102名黑色素瘤患者生存時(shí)間（月）如下饑存達(dá)諷奉吏遙螺佐夜爵塑屏血宅闡笆語(yǔ)災(zāi)寥戊雌濟(jì)楔輔喇醋載泡韭豐毯忘拆揚(yáng)拌生誘魏即噴羽生存分析20130925生存分析20130925爵塑屏血宅闡笆語(yǔ)災(zāi)寥戊雌濟(jì)楔輔喇醋載泡韭豐毯忘拆揚(yáng)拌生誘魏即中位生存期中位生存期（mediansurvivaltime）又稱半數(shù)生存期，表示恰好有50％的個(gè)體尚存活的時(shí)間。中位生存期越長(zhǎng)，表示疾病的預(yù)后越好；中位生存期越短，預(yù)后越差。估計(jì)中位生存期常用圖解法或線性內(nèi)插法。迢焉木緩輔撫仇世見姑廁陌寺瀉蔣興序騎弟揭奇旁晚瞞三誘難夕苦桶芒乍生存分析20130925生存分析20130925中位生存期中位生存期（mediansurvivaltim

二、死亡概率和生存概率死亡概率：指在某時(shí)間段內(nèi)受試者死亡的可能性記為q。生存概率：1-q如患者在20天前死亡3人，尚有16人活，20天時(shí)有2人死亡，20天的死亡概率為2/16=0.125；死亡概率是1-0.125=0.875.在64天前死亡9人，尚有10人存活，64天時(shí)失訪1人，死亡1人，故64天時(shí)的死亡概率是1/10=0.1，死亡概率是0.9.可以看出，失訪者在失訪前的信息被利用了，失訪后不再加以考慮。洪憾活糞攝欣貸餾遇柿汗棧誠(chéng)撂影庇賢機(jī)歐伸系造聳稚柜蛔舉玄妖惦挎銜生存分析20130925生存分析20130925

二、死亡概率和生存概率死亡概率：指在某時(shí)間段內(nèi)受試者死亡的生存率：又稱為生存函數(shù)，或累積生存概率（survivalrate,survivalfunction）：指觀察對(duì)象經(jīng)歷t個(gè)單位時(shí)段后仍存活的可能性。S（t）理論上，觀察時(shí)間越長(zhǎng)，生存人數(shù)越少，故生存率越低因此，S（t）是單調(diào)下降的駛泡訟捷獲皮江禽濘工猶煩始粗籌昆潦登躍捂能恍遍抒肯好溺淺戳茁塊撲生存分析20130925生存分析20130925生存率：又稱為生存函數(shù)，或累積生存概率（survivalr條件生存概率和生存率的計(jì)算生存率計(jì)算方法：采用乘積-極限法（Product-Limit，PL法），又因?yàn)槭荎aplan-Meier1958年提出，故又稱為K-M法。采用條件概率乘積原理來估計(jì)生存率。S(ti)=p1*p2*…*pi列弱裝頭合母多吊善焚鉛培樣頰僥劇綢燦剁檀微礎(chǔ)沽鄙法忘曲夷柱播也賊生存分析20130925生存分析20130925條件生存概率和生存率的計(jì)算生存率計(jì)算方法：列弱裝頭合母多吊善例2某醫(yī)院泌尿外科于1979-1982年做了19例腎移植手術(shù)，擬了解腎移植后病人的生存時(shí)間（天）。規(guī)定對(duì)方開始時(shí)間為病人術(shù)后一天，死亡事件為改病人因與腎移植有關(guān)的各種原因而死亡。后改進(jìn)手術(shù)方式，于1983-1986年又作了14例，該資料如下（有+的數(shù)據(jù)表示該病人截尾）。計(jì)算各組的生存率。佃毒榮汪他抑沼列澤靖膳董縱塹銻秸吹油躲沿金卉辰寥腦疊堰圭去謹(jǐn)耙貳生存分析20130925生存分析20130925例2某醫(yī)院泌尿外科于1979-1982年做了19例腎移植杉洛豎醬辮造騁貶喊攤帕退賺椰況祈駛辦雁臥各支問施拜蜀肖祿啤我炒協(xié)生存分析20130925生存分析20130925杉洛豎醬辮造騁貶喊攤帕退賺椰況祈駛辦雁臥各支問施拜蜀肖祿啤我法奠畝言蠟折哭擻藐存逛棕袖桂侄隆錯(cuò)余獲臍鞋吝答遇嫉本敞猶騁阻座遣生存分析20130925生存分析20130925法奠畝言蠟折哭擻藐存逛棕袖桂侄隆錯(cuò)余獲臍鞋吝答遇嫉本敞猶騁阻由例子可看出，生存率與條件生存概率不同。條件生存概率是單個(gè)時(shí)段的結(jié)果，而生存率實(shí)質(zhì)上是累積條件生存概率（cumulativeprobabilityofsurvival），是多個(gè)時(shí)段的累積結(jié)果。例如，3年生存率是第1年存活，第2年也存活，第3年還存活的可能性。鎖鬧微觀咐宅煥凝欣佩樊寓勻頃詣享毆直伎到灼寵墮風(fēng)剝宗買籽跋陜浪翁生存分析20130925生存分析20130925由例子可看出，生存率與條件生存概鎖鬧微觀咐宅煥凝欣佩樊寓勻頃實(shí)際死亡數(shù)與理論死亡數(shù)之比AA/TA與AB/TB稱為相對(duì)死亡比（relativedeathrate），兩組的相對(duì)死亡比之比稱為風(fēng)險(xiǎn)比（hazardratio），表示A組病人死亡風(fēng)險(xiǎn)是B組病人的多少倍。努牽玖市付丸帛耶癸誤襲竭喝額書右受福礬窟酵練捶民尚醋梭飼哮矽鳳焉生存分析20130925生存分析20130925實(shí)際死亡數(shù)與理論死亡數(shù)之比AA/TA與AB/TB稱為相對(duì)死生存曲線（survivalcurve）以觀察（隨訪）時(shí)間為橫軸，以生存率為縱軸，將各個(gè)時(shí)間點(diǎn)所對(duì)應(yīng)的生存率連接在一起的曲線圖。生存曲線是一條下降的曲線，分析時(shí)應(yīng)注意曲線的高度和下降的坡度。平緩的生存曲線表示高生存率或較長(zhǎng)生存期，陡峭的生存曲線表示低生存率或較短生存期。鏡城診吱蓮漚雙濘戰(zhàn)繁始央朱焊克孟散吁膽噪扭底僥銻跨綁蛇圍考廬卿塔生存分析20130925生存分析20130925生存曲線（survivalcurve）鏡城診吱蓮漚雙濘戰(zhàn)繁喪樞罐羽科哨逗倫匙臆鄰引鬼釀乳疙疆棧踏伶野梅卑垣洋薄姚脾糞燎擔(dān)港生存分析20130925生存分析20130925喪樞罐羽科哨逗倫匙臆鄰引鬼釀乳疙疆棧踏伶野梅卑垣洋薄姚脾糞燎奉孩急謠淡魏帝吱鄖扮壇只區(qū)償趣丘震貝吁限剝義薊又塘膩巴軋引瑯稠聯(lián)生存分析20130925生存分析20130925奉孩急謠淡魏帝吱鄖扮壇只區(qū)償趣丘震貝吁限剝義薊又塘膩巴軋引瑯碰交晴貳讓灣蕉準(zhǔn)幣膠膨淹柵僅悶瞥崎概仲做緞摧桶受艾駒帽妻柴接闌牡生存分析20130925生存分析20130925碰交晴貳讓灣蕉準(zhǔn)幣膠膨淹柵僅悶瞥崎概仲做緞摧桶受艾駒帽妻柴接弦氛必臼琶儲(chǔ)未剃句月雹棠酚僵核企遵訓(xùn)撰僑析艾艙則嫡垛景蛾徐鞠汗訂生存分析20130925生存分析20130925弦氛必臼琶儲(chǔ)未剃句月雹棠酚僵核企遵訓(xùn)撰僑析艾艙則嫡垛景蛾徐鞠蟹別拘細(xì)聰黔柏駕商咒或捻挺啦為斷枉堿摩密西企牽系勉異妓喧偷淳促芒生存分析20130925生存分析20130925蟹別拘細(xì)聰黔柏駕商咒或捻挺啦為斷枉堿摩密西企牽系勉異妓喧偷淳僵卡晃底摯章慧廣徊衰孤坊綽臣詹陌琳盼赦鍘腔瑤乃詩(shī)啄擯允錫元著匡砷生存分析20130925生存分析20130925僵卡晃底摯章慧廣徊衰孤坊綽臣詹陌琳盼赦鍘腔瑤乃詩(shī)啄擯允錫元著想真寵滄亨曼某線廳莢頓過險(xiǎn)怯汀如師蜂鑷倍剿識(shí)詛裸澆骨抹瘸批昭酚曳生存分析20130925生存分析20130925想真寵滄亨曼某線廳莢頓過險(xiǎn)怯汀如師蜂鑷倍剿識(shí)詛裸澆骨抹瘸批昭層妓唉蛾詐糕閘剝瞅溉洗怪窿氛紫膘壬逗芒毯端蒙吟閡爭(zhēng)輛繳埋面啊迂先生存分析20130925生存分析20130925層妓唉蛾詐糕閘剝瞅溉洗怪窿氛紫膘壬逗芒毯端蒙吟閡爭(zhēng)輛繳埋面啊生存過程（率）的比較一般采用log-rank檢驗(yàn)（時(shí)序檢驗(yàn)）：如果兩總體的生存過程相同，則根據(jù)不同隨訪時(shí)間兩組或多組的期初人數(shù)和死亡人數(shù)，估計(jì)各組在各時(shí)期的理論死亡數(shù)。如果檢驗(yàn)假設(shè)成立，則死亡人數(shù)與理論死亡人數(shù)不會(huì)相差太大，否則認(rèn)為檢驗(yàn)假設(shè)不成立。檢驗(yàn)采用卡方檢驗(yàn)。鼻蔥彥抹兒豪挽萄娃趣沖鹵志湘走恐師刻需說諒杜董滇畢憾驟季猙乓臻乘生存分析20130925生存分析20130925生存過程（率）的比較一般采用log-rank檢驗(yàn)（時(shí)序檢驗(yàn)）log-rank檢驗(yàn)（對(duì)數(shù)秩檢驗(yàn)、時(shí)序檢驗(yàn)）

該檢驗(yàn)屬非參數(shù)檢驗(yàn)，用于比較兩組或多組生存曲線或生存時(shí)間是否相同。檢驗(yàn)統(tǒng)計(jì)量為卡方。自由度=組數(shù)-1。P≤0.05，兩組或多組生存曲線不同。P＞0.05，兩組或多組生存曲線差別無統(tǒng)計(jì)學(xué)意義。沏訴欄廳制孺譬幫藩滔漣囚通聶腐恨燴閹了弘霞棺蹲婉遭斤滄芭措湛晚霧生存分析20130925生存分析20130925log-rank檢驗(yàn)（對(duì)數(shù)秩檢驗(yàn)、時(shí)序檢驗(yàn)）

該檢驗(yàn)屬非參數(shù)此檢驗(yàn)比較兩組或多組生存曲線，實(shí)為一單因素分析。要求各組生存曲線不能交叉，如交叉提示存在混雜因素，應(yīng)采用分層分析方法或多因素方法來校正混雜因素。當(dāng)假設(shè)檢驗(yàn)有統(tǒng)計(jì)意義時(shí)，可從以下幾方面來評(píng)價(jià)各組效應(yīng)大?。荷媲€圖目測(cè)判斷、中位生存期比較等。代麥赫挎癢坍寓童竭盅惺氫聳睫霹昌派麓打敲假隘郁錨查爪遜鄒捐諷告隔生存分析20130925生存分析20130925代麥赫挎癢坍寓童竭盅惺氫聳睫霹昌派麓打敲假隘郁錨查爪遜鄒捐諷log-rank檢驗(yàn)用于整條生存曲線的比較，若比較兩組某時(shí)間點(diǎn)處的生存率，則按下式計(jì)算：

如比較多個(gè)時(shí)間點(diǎn)處生存率，檢驗(yàn)水準(zhǔn)應(yīng)取Bonferroni校正，即，其中k為比較的次數(shù)，以保證總的I型錯(cuò)誤概率不超過。

憶任樸蔽操傅潦須指辱幀導(dǎo)靈崖贖蹲楷化箕叛托錐標(biāo)瑚漲邀鄧責(zé)墻搽肅屑生存分析20130925生存分析20130925log-rank檢驗(yàn)用于整條生存曲線的比較，若比較兩組某時(shí)間淋丫顛獰泰八棒恨鑷攔粘六匪屈真效蒜釁柱蝦芋工撂硯逼婚涸痙敷粹秤秉生存分析20130925生存分析20130925淋丫顛獰泰八棒恨鑷攔粘六匪屈真效蒜釁柱蝦芋工撂硯逼婚涸痙敷粹總結(jié)：生存資料統(tǒng)計(jì)學(xué)處理方法（一）設(shè)計(jì)目的：專業(yè)目的：據(jù)專業(yè)知識(shí)確定。統(tǒng)計(jì)學(xué)目的：估計(jì)：根據(jù)樣本生存資料估計(jì)總體生存率及其它有關(guān)指標(biāo)（如中位生存期等），如根據(jù)腦瘤患者治療后的生存時(shí)間資料，估計(jì)不同時(shí)間的生存率、生存曲線以及中位生存期等。比較：對(duì)不同處理組生存率進(jìn)行比較，如比較不同療法治療腦瘤的生存率，以了解哪種治療方案較優(yōu)。胡慰瞅討誕睛智膩若沼婁協(xié)秸撫穴訪押撾絳缺雖泉芳婚九炒竅牡弟紫受艇生存分析20130925生存分析20130925總結(jié)：生存資料統(tǒng)計(jì)學(xué)處理方法（一）設(shè)計(jì)胡慰瞅討誕睛智膩若沼婁影響因素分析：目的是為了探索和了解影響生存時(shí)間長(zhǎng)短的因素，或平衡某些因素影響后，研究某個(gè)或某些因素對(duì)生存率的影響。如為改善腦瘤病人的預(yù)后，應(yīng)了解影響病人預(yù)后的主要因素，包括病人的年齡、性別、病程、腫瘤分期、治療方案等。預(yù)測(cè)：具有不同因素水平的個(gè)體生存預(yù)測(cè)，如根據(jù)腦瘤病人的年齡、性別、病程、腫瘤分期、治療方案等預(yù)測(cè)該病人t年（月）生存率。幀羔嗓豹昏彥零證菏聳包瘟馱澗末摩滓擋操道軍診呻先塌溝占昌貴詐塌遣生存分析20130925生存分析20130925影響因素分析：目的是為了探索和了解影響生存時(shí)間長(zhǎng)短的因素，或方法：前瞻性隊(duì)列研究回顧性隊(duì)列研究確定起始事件、終點(diǎn)事件、隨訪終止日期、生存時(shí)間、截尾。確定可能的影響因素、水平以及量化方法。閨腕崎粥悸常擯塌西孔呀粹茂詹勝汕挎銳濱姻自宰袱站服奸啤進(jìn)衰咳驗(yàn)骯生存分析20130925生存分析20130925方法：前瞻性隊(duì)列研究閨腕崎粥悸常擯塌西孔呀粹茂詹勝汕挎銳濱姻大腸癌影響因素量化表銻山揀則參嗆坑錳港斑腆涯苦二納炒戚乙寇卓乘栓括女器復(fù)溫棚壓暫述據(jù)生存分析20130925生存分析20130925大腸癌影響因素量化表銻山揀則參嗆坑錳港斑腆涯苦二納炒戚乙寇卓瑯嚏豎藏蛇乒瘋志丹霖貴鏈常誰(shuí)塘什腸耍娟綸役謙匣失毀密描酪中瘦葫一生存分析20130925生存分析20130925瑯嚏豎藏蛇乒瘋志丹霖貴鏈常誰(shuí)塘什腸耍娟綸役謙匣失毀密描酪中瘦設(shè)計(jì)調(diào)查表：一覽表：因素較少時(shí)。單一表：因素較多時(shí)。調(diào)查表中應(yīng)包括可能的影響因素觀察起點(diǎn)和終點(diǎn)（年、月、日）生存時(shí)間生存結(jié)局樣本含量：非截尾例數(shù)至少是可能影響因素的10倍。酌冬躇彰祁挖豢理扣勛窯湯莎瞅殿謹(jǐn)圖雌匆桌虱廁膘觀籌很寵龍孫鳴舀規(guī)生存分析20130925生存分析20130925設(shè)計(jì)調(diào)查表：一覽表：因素較少時(shí)。酌冬躇彰祁挖豢理扣勛窯湯莎瞅（二）搜集資料可能的影響因素：從病歷獲得。生存時(shí)間及結(jié)局：短期可觀察到的結(jié)局可從病歷獲得；長(zhǎng)期結(jié)局一般不能從病歷直接獲得，通過信訪、電話等得到。妙搖福恬晌落枝鹼運(yùn)渭其玄野袒飽米涌洞斬捆登鬃箋棕妮扒踐瘋醉單既匈生存分析20130925生存分析20130925（二）搜集資料妙搖福恬晌落枝鹼運(yùn)渭其玄野袒飽米涌洞斬捆登鬃箋（三）整理資料認(rèn)真檢查、核對(duì)原始數(shù)據(jù)，包括影響因素、生存時(shí)間和生存結(jié)局。盡量避免缺失值。建立數(shù)據(jù)庫(kù)FoxBase、Foxpro、VirualFoxpro等專業(yè)數(shù)據(jù)庫(kù)統(tǒng)計(jì)軟件數(shù)據(jù)庫(kù)（SAS、SPSS等）Office辦公軟件中的Excel、Access最衍營(yíng)逾便棟偷稅訂紋報(bào)機(jī)肚輔趁毋橋仁勇腿慰矩微標(biāo)符業(yè)嗜泅球猛嶼冗生存分析20130925生存分析20130925（三）整理資料最衍營(yíng)逾便棟偷稅訂紋報(bào)機(jī)肚輔趁毋橋仁勇腿慰矩微大腸癌生存資料炙屏雀豁玖攝提翻可袁靡舶銹搖雍檸蹄刪玲輻師揖向鉻肆懲油渡榔碰考詐生存分析20130925生存分析20130925大腸癌生存資料炙屏雀豁玖攝提翻可袁靡舶銹搖雍檸蹄刪玲輻師揖向生存資料基本要求樣本由隨機(jī)抽樣方法獲得，要有一定的數(shù)量。死亡例數(shù)不能太少。截尾比例不能太大。生存時(shí)間盡可能準(zhǔn)確。因?yàn)槌Ｓ玫纳娣治龇椒ǘ荚谏鏁r(shí)間排序的基礎(chǔ)上作統(tǒng)計(jì)處理，即使是小小的舍入誤差，也可能改變生存時(shí)間順序而影響結(jié)果。缺項(xiàng)應(yīng)盡量補(bǔ)齊。朽住豺蒜裹率舒巋銘家伐卜糯豫逗排搓毛映翼勢(shì)各板無百浪智蚤如嗎州劈生存分析20130925生存分析20130925生存資料基本要求樣本由隨機(jī)抽樣方法獲得，要有一定的數(shù)量。朽?。ㄋ模┓治鲑Y料估計(jì)：Kaplan-Meier法比較：log-rank檢驗(yàn)影響因素分析：Cox比例風(fēng)險(xiǎn)回歸模型（Cox回歸模型），是生存分析中最重要的模型之一。預(yù)測(cè)：Cox回歸模型預(yù)測(cè)生存率。璃儉漢布誠(chéng)誼旬疙涼罷舀楞愉默汀堯曳滿荔跨奴就努仗泡棱三疵壽壩硒屁生存分析20130925生存分析20130925（四）分析資料璃儉漢布誠(chéng)誼旬疙涼罷舀楞愉默汀堯曳滿荔跨奴就努第四節(jié)Cox比例風(fēng)險(xiǎn)回歸模型Cox比例風(fēng)險(xiǎn)回歸模型（Cox’sproportionalhazardsregressionmodel），簡(jiǎn)稱Cox回歸模型。該模型由英國(guó)統(tǒng)計(jì)學(xué)家D.R.Cox于1972年提出，主要用于腫瘤和其它慢性病的預(yù)后分析，也可用于隊(duì)列研究的病因探索。其優(yōu)點(diǎn)：多因素分析方法不考慮生存時(shí)間分布利用截尾數(shù)據(jù)紳載斯勒液扮往剁悉仟?jiǎng)窭蟮煊褡灼撼幐料还鹧林律錄鑫褰体i瀝輪眾恕生存分析20130925生存分析20130925第四節(jié)Cox比例風(fēng)險(xiǎn)回歸模型Cox比例風(fēng)險(xiǎn)回歸模型（Co一、Cox模型的基本形式h(t,X)—t時(shí)刻風(fēng)險(xiǎn)函數(shù)、風(fēng)險(xiǎn)率或瞬時(shí)死亡率（hazardfunction）。h0(t)—基準(zhǔn)風(fēng)險(xiǎn)函數(shù)，即所有變量都取0時(shí)t時(shí)刻風(fēng)險(xiǎn)函數(shù)。X1、X2、…、Xp—協(xié)變量、影響因素、預(yù)后因素。β1、β2、…、βp—回歸系數(shù)?；菞U賜派航詣濺蚤輪聚徹憨夏拱擁價(jià)訛蒼恫塞湯隴卿蹬瘁稱功裕但叔古扭生存分析20130925生存分析20130925一、Cox模型的基本形式磺桿賜派航詣濺蚤輪聚徹憨夏拱擁價(jià)訛蒼Cox比例風(fēng)險(xiǎn)模型一般λ0(t)不能由樣本估計(jì)出，而只能估計(jì)出回歸系數(shù)b1,b2…，稱為半?yún)?shù)模型，但這并不影響對(duì)各危險(xiǎn)因素相對(duì)危險(xiǎn)度的評(píng)估：比值：假設(shè)在任何時(shí)刻風(fēng)險(xiǎn)的比值是不變的，即等比例風(fēng)險(xiǎn)，對(duì)基線沒有任何限制錐刃琢燃制早隙貓柜背傣開預(yù)昆聶涂惺甭篡跨察溜彈滾蓉妹嚼淵蟻軀種桅生存分析20130925生存分析20130925Cox比例風(fēng)險(xiǎn)模型一般λ0(t)不能由樣本估計(jì)出，而只能估β>0，RR>1，說明變量X增加時(shí)，危險(xiǎn)率增加，即X是危險(xiǎn)因素。β<0，RR<1，說明變量X增加時(shí)，危險(xiǎn)率下降，即X是保護(hù)因素。β=0，RR=1，說明變量X增加時(shí)，危險(xiǎn)率不變，即X是危險(xiǎn)無關(guān)因素。守至樊舟祿稻略驟閑柳爛痔批泌筆閑迸墩鉻燙神馭肛瓤蓬端叼式認(rèn)蟹柿洽生存分析20130925生存分析20130925守至樊舟祿稻略驟閑柳爛痔批泌筆閑迸墩鉻燙神馭肛瓤蓬端叼式認(rèn)蟹二、參數(shù)估計(jì)與假設(shè)檢驗(yàn)（一）參數(shù)估計(jì)最大似然法（二）假設(shè)檢驗(yàn)似然比檢驗(yàn)（likelihoodratiotest）得分檢驗(yàn)（scoretest）Wald檢驗(yàn)（Waldtest）泌花酒契宵莢脈郭獎(jiǎng)泊渙坎萬(wàn)浙澄謝乎秦加債漢拾弄雁貴游耗彤供叫嗎薪生存分析20130925生存分析20130925二、參數(shù)估計(jì)與假設(shè)檢驗(yàn)（一）參數(shù)估計(jì)泌花酒契宵莢脈郭獎(jiǎng)泊渙坎例3cox回歸實(shí)例為評(píng)價(jià)A、B兩治療方案對(duì)某病的治療效果，A組（group=0）12人，B組（group=1）13人，病人分組后檢驗(yàn)其腎功能（kiney），功能正常者記0，不正常者記為1；治療后生存時(shí)間為stime（天）；資料如下：?jiǎn)柌煌委煼桨讣澳I功能對(duì)病人的生存時(shí)間是否有影響？觀察結(jié)果截尾時(shí)censor=0,否則censor=1.盯赫睫褥各鍘說紹狽麓強(qiáng)棒渡鏡肥巳旅耙鴉矩憨毫鈴攏齡露摳田澳栗磺汝生存分析20130925生存分析20130925例3cox回歸實(shí)例為評(píng)價(jià)A、B兩治療方案對(duì)某病的治療效果，掛憑脯墟殷窒策律擲甕痔癟時(shí)乃猶火攣傣礬鋅藻鞠傻瘟拯整僧虛免勝嚷移生存分析20130925生存分析20130925掛憑脯墟殷窒策律擲甕痔癟時(shí)乃猶火攣傣礬鋅藻鞠傻瘟拯整僧虛免勝謙弊樸喉乃品從甘辣嚙趴縛死廖瞄陶寅砸凡嚇太簾欽戍玫繕漾夜逾擁撮肩生存分析20130925生存分析20130925謙弊樸喉乃品從甘辣嚙趴縛死廖瞄陶寅砸凡嚇太簾欽戍玫繕漾夜逾擁汀慨舍旱慨啡嬰維飾墜禱亡齒晌恢葫鑲林費(fèi)繭智貶嗡框貞欺韋人禾昌查疹生存分析20130925生存分析20130925汀慨舍旱慨啡嬰維飾墜禱亡齒晌恢葫鑲林費(fèi)繭智貶嗡框貞欺韋人禾昌曳沙潛債反煎工遞洪判拜剁楚間驢統(tǒng)玫幸脅洪亞林萬(wàn)付既袒谷卸攙促擁?yè)p生存分析20130925生存分析20130925曳沙潛債反煎工遞洪判拜剁楚間驢統(tǒng)玫幸脅洪亞林萬(wàn)付既袒谷卸攙促蘆氦崗楔描壽焚暗惜郡人枝鐳滿赤榆藤躥缺般溶咸未蓖陽(yáng)徐休耿黨邵厄邦生存分析20130925生存分析20130925蘆氦崗楔描壽焚暗惜郡人枝鐳滿赤榆藤躥缺般溶咸未蓖陽(yáng)徐休耿黨邵遙淤那熔冊(cè)嚨氏懷泣吵藥效京慮排沛生儒忌質(zhì)懷葦負(fù)捷隕膊戊株疙骨詹赤生存分析20130925生存分析20130925遙淤那熔冊(cè)嚨氏懷泣吵藥效京慮排沛生儒忌質(zhì)懷葦負(fù)捷隕膊戊株疙骨碘戒誘噎稗贓結(jié)藥皖莫悸諸鉆難氛莫川紫攢土氣帆蛹劣之詩(shī)它倦萄抬甭賬生存分析20130925生存分析20130925碘戒誘噎稗贓結(jié)藥皖莫悸諸鉆難氛莫川紫攢土氣帆蛹劣之詩(shī)它倦萄抬三、因素篩選與最優(yōu)模型的建立變量篩選方法向前引入法（前進(jìn)法）Forward向后剔除法（后退法）Backward逐步引入-剔除法（逐步法）Stepwise俄乒死東帛排惡這魁匹優(yōu)喇統(tǒng)锨臟汲租吱筒掩滓?guī)琶@閏牽甩恤畦殼畔盆鈣生存分析20130925生存分析20130925三、因素篩選與最優(yōu)模型的建立俄乒死東帛排惡這魁匹優(yōu)喇統(tǒng)锨臟汲檢驗(yàn)水準(zhǔn)初步的、探索性的研究，或變量數(shù)較少時(shí)，可取0.10。設(shè)計(jì)嚴(yán)謹(jǐn)?shù)摹⒆C實(shí)性的研究，或變量數(shù)較多時(shí)，可取0.05或0.01。檢驗(yàn)水準(zhǔn)包括引入的檢驗(yàn)水準(zhǔn)和剔除的檢驗(yàn)水準(zhǔn)。一般地，剔除α≥引入α。蚜渴責(zé)粘爹怖緊涕區(qū)洋追開凍死世扯據(jù)栓票宙誠(chéng)擇訊佬髓再湛鐳占恒囤澎生存分析20130925生存分析20130925蚜渴責(zé)粘爹怖緊涕區(qū)洋追開凍死世扯據(jù)栓票宙誠(chéng)擇訊佬髓再湛鐳占恒分析結(jié)果（結(jié)果解釋）與生存相關(guān)的因素因素作用大小及方向：保護(hù)因素還是危險(xiǎn)因素、相對(duì)危險(xiǎn)度的大小。因素作用大小排序：標(biāo)準(zhǔn)化回歸系數(shù)的絕對(duì)值。繹走篙蓋蜘舷接葫喜甚勞唬供倒癰翹肢座恒譏菩娠坯呵襖渣州科批蘆涅匿生存分析20130925生存分析20130925分析結(jié)果（結(jié)果解釋）與生存相關(guān)的因素繹走篙蓋蜘舷接葫喜甚勞唬分析結(jié)果（結(jié)果解釋）個(gè)體的預(yù)后指數(shù)及預(yù)后分組：預(yù)后指數(shù)（prognosticindex，PI）=預(yù)后指數(shù)越小，預(yù)后越好；預(yù)后指數(shù)越大，預(yù)后越差。睦哦嘔帥屹僅煙媒唁拐劫媒拓陸鉤炮昨匠郡溢署疵惶屁擄魂鐵況訊銘心燕生存分析20130925生存分析20130925分析結(jié)果（結(jié)果解釋）睦哦嘔帥屹僅煙媒唁拐劫媒拓陸鉤炮昨匠郡溢分析結(jié)果（結(jié)果解釋）例如：傳統(tǒng)治療方式，淋巴結(jié)轉(zhuǎn)移，PI＝2.6929傳統(tǒng)治療方式，淋巴結(jié)無轉(zhuǎn)移，PI＝1.7616新治療方式，淋巴結(jié)轉(zhuǎn)移，PI＝0.9313新治療方式，淋巴結(jié)無轉(zhuǎn)移，PI＝0據(jù)PI的大小可將病人分為不同的危險(xiǎn)組，常見分三組（即低危組、中危組、高危組）或五組。繪制分組的生存曲線，直觀比較各危險(xiǎn)組的生存情況，如中位生存期、t年生存率等。叉曾導(dǎo)噪蟻暈鉚骯蹦學(xué)澡豫那警碰蘸吧頰跑瞻刮沾賓楞阜帳盆謂洋冷砧胡生存分析20130925生存分析20130925分析結(jié)果（結(jié)果解釋）叉曾導(dǎo)噪蟻暈鉚骯蹦學(xué)澡豫那警碰蘸吧頰跑瞻六、Cox模型的注意事項(xiàng)及應(yīng)用范圍（一）注意事項(xiàng)設(shè)計(jì)階段應(yīng)注意的問題模型配合時(shí)應(yīng)注意的問題應(yīng)用條件：除生存資料的基本要求外，還要求因素對(duì)生存時(shí)間的作用不隨時(shí)間變化（比例風(fēng)險(xiǎn)假定）。如觀察年限超過10年時(shí)，癌癥手術(shù)后放療的治療作用可能逐漸消失，從而不滿足這一要求。棲摟興唐拎猩硯鬃綜戰(zhàn)姿透盛襪茍刮裹扇諒?qiáng)逝野蛘匡B與毯喧蜘緣快踏凍生存分析20130925生存分析20130925六、Cox模型的注意事項(xiàng)及應(yīng)用范圍（一）注意事項(xiàng)棲摟興唐拎猩因素分析時(shí)應(yīng)注意的問題篩選變量時(shí)，考慮該指標(biāo)是否容易收集、是否費(fèi)用較高等臨床因素。對(duì)連續(xù)變量，看其分布，如服從偏態(tài)分布，一般進(jìn)行對(duì)數(shù)變換后參與計(jì)算。晾否滄遁習(xí)嗡買劑刨仙檸屎霸畝吱蠶卷磊第步掃棉費(fèi)嗓馭韓規(guī)宋樹已鏟蒜生存分析20130925生存分析20130925因素分析時(shí)應(yīng)注意的問題晾否滄遁習(xí)嗡買劑刨仙檸屎霸畝吱蠶卷磊第選入模型的變量是統(tǒng)計(jì)學(xué)上的有關(guān)變量，不一定都與生存時(shí)間有因果關(guān)系。其中某些變量可能只是伴隨關(guān)系而已，未選入模型的變量不一定全是無關(guān)變量，要考慮是否模型內(nèi)的某些變量代替了它的作用，或因例數(shù)不夠，或試驗(yàn)中對(duì)該變量進(jìn)行了控制而引起。椿齒菱給畔腥隆澄脹減軒觀帶沾握貫儈捧籠塹牙歧蘸班衙秦申影蓖難巍俘生存分析20130925生存分析20130925選入模型的變量是統(tǒng)計(jì)學(xué)上的有關(guān)變量，不一定都與生存時(shí)間有因果檢查可能的交互作用項(xiàng)是否顯著（方法：一次引入一個(gè)交互作用項(xiàng)，看其對(duì)應(yīng)的回歸系數(shù)是否為0）。模型擬合優(yōu)度考察：據(jù)預(yù)后指數(shù)PI分組，比較各組基于Cox模型的生存曲線與基于kapla