生存及cox分析

《生存及cox分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、生存及COX分析SAS第十课：随访资料的生存分析一一非参数法与半参数cox比例风险模型m晓燕分享生存分析方法人体上可分为三类：非参数法.参数法和半参数方法，与之相对应，SAS收藏座］團☆AJ提供了三个程序步用于生存分析，它们是：LIFETEST过程提供非参数分析方法，用乘积极限法(Productlimitmethod)和寿命表法(Lifetablemethod)估计牛:存率和中位生存时间等；用対数秩检验(Log-ranktest)>Wilcoxon检验和似然比检验等做分组比较。该过程主耍用于估计半存率及进行单因素分析。IJFEREG过稈提

2、供指数模型、WeibulI模型、Gompertz模型等参数分析方法。PHREG过程提供半参数Cox比例风险模型分析。木章只介绍常用的IJFETEST过程和PHREG过程。§10.1LIFETEST过程及其应用10.1.1语法格式LIFETEST过程的语法格式如下：PROCLIFETEST[选项];TIME〈生存时间变量*截尾指示变量(数值)>;[TEST〈分组变量名列〉；STRATA〈分组变量名列〉；FREQ〈变量名列〉；BY〈变量名列>;]10.1.2语法说明PROC和TIME语句为必需的，其他语句都町以省略。TIME语句，为必需语句，

3、定义生存时间和截尾指示变量【过程选项】1.METHOD二方法指定佔计生存率所用的方法：PL要求用乘积极限法(即Kaplam-Meier法)佔计生存率并计算中位生存时间等，为缺省方法。LT要求用寿命表法估计生存率等。INTERVALS=(初值TO终值BY步长)只能在指定分析方法为寿命表法时使用。川寿命表法分析吋，程序会白动给定生存时间的区间。如果人为规定生存吋间的分组区间，则需用该选项指定。步长的缺省值为1。2.PLOTS二绘图类型耍求输出生存分析图。可供输出的图形冇：S对生存函数S(t)做图，横、纵坐标分别为t、S(t)oLS对-LOGS

4、(t)做图，横、纵坐标分别为t、-L0GS(t)oLLS对LOG(-LOGS(t))做图，横、纵坐标分别为LOG(t)、LOG(-L0GS(t))oII对风险函数做图，横、纵坐标分别为t、H(t)o【TIME语句】TIME语句川于定义生存时间和截尾指示变量。对截尾指示变量可以指定发生失效事件的数值，默认失效事件用0来表示，截尾事件用1来表示。[STRATA语句和TEST语句】STRATA语句定义牛存率比较的分组变量，TEST语句处义牛存率比较的分组变量或

5、办变量。STRATA语句在这里的作用和BY语句类似，都是要求按分组变量名列进行分析，

6、在计算生存率时各组分开计算；TEST语句定义需检验的变量，即生存时间与该变量是否有关，如果它后而定义的变量为数值变量，则把该变量当作协变量检验与生存时间的关系。如果它定义的为分组变量，则分组比较生存时间有无差别。10.1.3应用实例例10.1观察两组卵巢腺癌患者的病程天数如下。请用対数秩检验比较两组的牛:存期差异有无统计学意义，并作生存率曲线。（医统P.321,7.2题）A组（低恶性高分化癌）：2829175195309377+393+421+447+452709+744+770+1106+1206++B组（高恶性低分化癌）：348813

7、7199280291299300309351358369370375382392429+4511119+解：程序如下：dataa.yt7_2;inputtift<0thencensor=l;elsecensor二0;if_N_<16thengroup二'A';elsegroup二'B';t=abs（t）;cards;2829175195309-377......382392-429451-1119proclifetestmethod=plplots=(s,Is);timet*censor(1);stratagroup;proclifete

8、stmethod=ltplots=(s,h);timet*censor(1);stratagroup;tun;本例为不分组资料，生存时间有截尾值，应用乘积极限法和寿命表法分别作生存分析，结局变量只有一个，即生存时间，负值表示截尾，用if语句完成转换，产生新变量CENSOR,其取值1表示截尾，0表示不截尾。用两个PROC语句，分别用不同的方法作生存分析，按GROUP变量分组统计。IIFETEST过程的默认输出结果很多，主要有：%1生存率及其标准误。%1生存时间四分位数。%1生存时间的均值及其标准误。%1Wilcoxcon2统计量。当生存时间

9、符合对数正态分布时该检验统计效率较高。⑤Log-rank2统计量。当牛:存时间符合Wcibul1分布或属于比例风险模型时该检验统计效率较高。本例程序运行的主要结果如2第一部分，PL法对牛存资料