生存分析.SAS.ppt

:研究者从某个规定的时刻开始对研究对象作随访,观察某种规定事件的发生。结局:死于本病---失效病例失访---删失病例或截尾资料(终检值)死于它病---删失病例或截尾资料(终检值)目的:估计生存时间(失效时间),---=数值设置可信区间估计的置信水平METHOD=PL|KM|LT|LIFE指定估计生存函数的方法,默认为Kaplan-Meier法(KM选项),PL为乘积极限法(相当于Kaplan-Meier),LT和LIFE都是用寿命表法NINTERVAL=数值指出用寿命表法计算生存函数时的区间个数。WIDTH=数值指出用寿命表法计算生存函数时的区间宽度INTERVALS=指定计算区间的终点(寿命表法)=(S|LS|P)输出各种估计量与时间或对数时间的图;S为生存率与时间作图;LS为-log(生存率)与时间作图;当选用寿命表法估计生存函数时,可用P输出概率密度函数与时间作图。可同时选多个图,如PLOTS=(S,LS)FORMCHAR(1,2,7)=‘字符串’用于定义图形的横轴和纵轴。1对应的字符表示纵轴,2对应的字符表示横轴,两者也用于表示生存曲线的各阶;7对应的字符表示坐标点。-1:optionsps=40;/*设定输出结果页的大小为40行*/datalife1;inputdays@@;censor=(days<0);/*设定删失标志:括号中为逻辑表达式,days<0如为真,则变量censor值为1,表示删失;否则为0,表示失效*/group=(_n_>12)+1;/*通过逻辑判别式进行分组:_n_为内部变量,表示记录序号。*/days=abs(days);/*取绝对值*/(续)cards;852586363220365452496-528-560-67613182370761801952**********-490-540;proclifetestplots=(S)FORMCHAR(1,2,7)='|-*';/*要求对生存率和时间作图*/timedays*censor(1);/*时间变量为days,变量censor值为1的为删失记录*/stratagroup;/*按group分层分析,可比较层间一致性(齐性)*/run;,当有其他可能对预后有影响的因素(协变量)存在时,还要考虑协变量的效应例15-2采用15-1的数据,增加一个协变量(有无肾损害)-2datalife2;inputdaysrenal@@;censor=(days<0);group=(_n_>12)+1;days=abs(days);cards;815205816316312200365045204960-5280-5600-6760131181231700760180019502100232030003960-4900-5400;proclifetestplots=(s);timedays*censor(1);stratagroup;testrenal;/*要求检验renal与生存时间是否有关*/run;