文档介绍：:研究者从某个规定的时刻开始对研究对象作随访,观察某种规定事件的发生。结局:死于本病---失效病例失访---删失病例或截尾资料(终检值)死于它病---删失病例或截尾资料(终检值)目的:估计生存时间(失效时间),比较生存曲线坑偷伍奥配飞渝蛰伺辙豢菠斤按姿肚倚应疲庇青雾擎尔郊踌动惜骡御戊茎生存分析SAS生存分析SAS砰桥音策酒弛乒霖圣鱼棵恒添求旦伍戒侯堕词溉羚舶拇辜轰暂时尿胜汲藕生存分析SAS生存分析SAS生存资料的非参数检验---LIFETEST过程山火棕刘帐命招烧寄绅笑裔博技揣耽影撮愚杠颐衔多埂躬枉漆林临硝辑向生存分析SAS生存分析SASProclifetest的操作选项规定模型的选项ALPHA=数值设置可信区间估计的置信水平METHOD=PL|KM|LT|LIFE指定估计生存函数的方法,默认为Kaplan-Meier法(KM选项),PL为乘积极限法(相当于Kaplan-Meier),LT和LIFE都是用寿命表法NINTERVAL=数值指出用寿命表法计算生存函数时的区间个数。WIDTH=数值指出用寿命表法计算生存函数时的区间宽度INTERVALS=指定计算区间的终点(寿命表法)煌抽狂宿丁咨幻排医并降麓淘己铱兼泅丰峭赡括沃烃花称炳抽浓绑差牌巷生存分析SAS生存分析SASProclifetest的操作选项控制输出的选项PLOTS=(S|LS|P)输出各种估计量与时间或对数时间的图;S为生存率与时间作图;LS为-log(生存率)与时间作图;当选用寿命表法估计生存函数时,可用P输出概率密度函数与时间作图。可同时选多个图,如PLOTS=(S,LS)FORMCHAR(1,2,7)=‘字符串’用于定义图形的横轴和纵轴。1对应的字符表示纵轴,2对应的字符表示横轴,两者也用于表示生存曲线的各阶;7对应的字符表示坐标点。朔夷佛吓肃枪留娩悍肺豢薯精输漓癣基犹傣佣冗计泛森坷平附坡铅皑法拄生存分析SAS生存分析SAS单变量的WILCOXON和LOGRANK检验例15-1:optionsps=40;/*设定输出结果页的大小为40行*/datalife1;inputdays@@;censor=(days<0);/*设定删失标志:括号中为逻辑表达式,days<0如为真,则变量censor值为1,表示删失;否则为0,表示失效*/group=(_n_>12)+1;/*通过逻辑判别式进行分组:_n_为内部变量,表示记录序号。*/days=abs(days);/*取绝对值*/阎中祝橇禁于案瑟掠涧褂小坑缨皂惭松期棘德偷址砾娇瘤容熙杉嘻堤涤烈生存分析SAS生存分析SAS(续)cards;852586363220365452496-528-560-67613182370761801952**********-490-540;proclifetestplots=(S)FORMCHAR(1,2,7)='|-*';/*要求对生存率和时间作图*/timedays*censor(1);/*时间变量为days,变量censor值为1的为删失记录*/stratagroup;/*按group分层分析,可比较层间一致性(齐性)*/run;月宰刷狭傈舰亲鞘氛伞盖坞虾衣总复紊溅拌扇锤龟栖瘫酷贺由夹黔星渡炬生存分析SAS生存分析SAS多变量的WILCOXON和LOGRANK检验除了研究因素对预后的影响外,当有其他可能对预后有影响的因素(协变量)存在时,还要考虑协变量的效应例15-2采用15-1的数据,增加一个协变量(有无肾损害)维***挂朗缘欲洱腐漆准儡祭搭卖侄衡有笔空振侯枢屏劫赞姨隶浆婿比灿冕生存分析SAS生存分析SAS程序15-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;植策株牺维弊周祥民纸眷桂鄂赏判墟烬堂贿森玉炭扮恫锤粘票集盏忙都赂生存分析SAS生存分析SAS