生存資料有哪些常用的分析方法?

SP的工作中經(jīng)常會(huì)遇到生存分析,臨床試驗(yàn)常用的生存分析方法包括:對(duì)生存過(guò)程進(jìn)行描述的Kaplan-Meier法,對(duì)生存率進(jìn)行組間比較的log-rank檢驗(yàn), 以及對(duì)生存過(guò)程的影響因素進(jìn)行分析的Cox比例風(fēng)險(xiǎn)模型。
1、Kopion-Meier法
由Kaplan 和Meier于1958年提出的,故稱(chēng)Kopion-Meier法,又稱(chēng)乘積限估計(jì)法(productlimit method, 簡(jiǎn)稱(chēng)PL法)。主要用于生存資料的生存率估計(jì), 是一種非參數(shù)估計(jì)方法。其基本原理是基于概率的乘法法則,通過(guò)各時(shí)期的人群的生存概率的累計(jì)相乘得到活過(guò)一定時(shí)期的生存率,從而計(jì)算累計(jì)死亡率。以生存率(或死亡率)為縱坐標(biāo),時(shí)間為橫坐標(biāo).并以階梯形的折線(xiàn)連接各點(diǎn),得到的圖為生存率(或死亡率)的Kaplan-Meier 圖。相應(yīng)的SAS 程序如下:
ODS GRAPHICS ON;

PROC LIFETEST DATA=MyData;
  /*繪制KM生存曲線(xiàn)。顯示指定觀察時(shí)間時(shí)的樣本量和假設(shè)檢驗(yàn)結(jié)果*/
  PLOTS=(SURVIVAL(ATRISK=0 TO 20 BY 2)TEST)MAXTIME=20;
  TIME Stime*Outcome(0);
  STRATA Group;
RUN;

ODS GRAPHICS OFF;
2、log-rank檢驗(yàn)
log-rank檢驗(yàn)又稱(chēng)時(shí)序檢驗(yàn)。其基本思想是:如果兩個(gè)或多個(gè)總體的生存過(guò)程相同,則根據(jù)不同隨訪(fǎng)時(shí)間兩組或多組的期初人數(shù)和死亡人數(shù),估計(jì)各組在各時(shí)期的理論死亡數(shù)。如果檢驗(yàn)假設(shè)成立,則實(shí)際死亡數(shù)與理論死亡數(shù)不會(huì)相差太大,否則認(rèn)為檢驗(yàn)假設(shè)不成立。實(shí)際死亡數(shù)與理論死亡數(shù)的比較用卡方檢驗(yàn)。相應(yīng)的SAS 程序如下:
PROC LIFETEST DATA=Mydata;
  TIME Stime*Outcome(0);*stime為生存時(shí)間,outcome為指標(biāo)變量,此處0表示截尾;
  STRATA Group;         *分組變量;
RUN;
3、 Cox比例風(fēng)險(xiǎn)模型
英國(guó)統(tǒng)計(jì)學(xué)家Cox用偏似然原理巧妙地回避了求基線(xiàn)風(fēng)險(xiǎn)而解決估計(jì)回歸系數(shù)的問(wèn)題,無(wú)須對(duì)基線(xiàn)風(fēng)險(xiǎn)作任何限制,相應(yīng)的風(fēng)險(xiǎn)模型被稱(chēng)為Cox 比例風(fēng)險(xiǎn)模型(Coxproportional hazards model),又稱(chēng)Cox 回歸。Cox回歸的相應(yīng)SAS程序如下:
PROC PHREG DATA=MYDATA;
  MODEL Stime*Outcome(0)=Group/RL; *RL計(jì)算可信區(qū)間;
RUN;

請(qǐng)前往:http://lygongshang.com/TeacherV2.html?id=166