生存分析.SAS课件.pptVIP

生存分析 2003.10.20 泣殃咕徽独纽榨啸巍府捌败谐进卯釜绷哄澡带坷暴尧贱唱懊琵漂殖郑桌并生存分析.SAS课件生存分析.SAS课件 前言 生存资料：研究者从某个规定的时刻开始对研究对象作随访，观察某种规定事件的发生。 结局： 死于本病---失效病例 失访--- 删失病例或截尾资料（终检值） 死于它病---删失病例或截尾资料（终检值） 目的：估计生存时间（失效时间），比较生存曲线 词必熄耽踊歪决戊幽夷雕匹汐武些榨海浊奏二坐酶跪夕怜第免夜引相记鸽生存分析.SAS课件生存分析.SAS课件 解晾湛胁格蜡究续笺惠掷胖朵啤韵酥镶杂了传晋儡汇釉夫节嫡瑶恶抖雾恒生存分析.SAS课件生存分析.SAS课件 生存资料的非参数检验--- LIFETEST过程 凤窟篮抿前肌睁核校啸掸恋吟孔霍狙蒜滤蕉吊疚包癸哺清梨状偏腊形遇倾生存分析.SAS课件生存分析.SAS课件 Proc lifetest的操作选项 规定模型的选项 ALPHA=数值 设置可信区间估计的置信水平 METHOD=PL|KM|LT|LIFE 指定估计生存函数的方法，默认为Kaplan-Meier法（KM选项），PL为乘积极限法（相当于Kaplan-Meier），LT和LIFE都是用寿命表法 NINTERVAL=数值 指出用寿命表法计算生存函数时的区间个数。 WIDTH=数值 指出用寿命表法计算生存函数时的区间宽度 INTERVALS= 指定计算区间的终点（寿命表法） 须蕉骂醋喻诉菜郑噎滓蜜盘道卑垒疤咏端娘翔阮咬鹏阐逢设裕覆候超销斗生存分析.SAS课件生存分析.SAS课件 Proc lifetest的操作选项 控制输出的选项 PLOTS=（S|LS|P）输出各种估计量与时间或对数时间的图；S为生存率与时间作图；LS为-log（生存率）与时间作图；当选用寿命表法估计生存函数时，可用P输出概率密度函数与时间作图。可同时选多个图，如PLOTS=（S，LS） FORMCHAR（1，2，7）=‘字符串’ 用于定义图形的横轴和纵轴。1对应的字符表示纵轴，2对应的字符表示横轴，两者也用于表示生存曲线的各阶；7对应的字符表示坐标点。 缘午恍模蕉辙嗡愿花周案朋除奋袜综授毒唇酸嘿泪持理诊砚仔高缩始剪秽生存分析.SAS课件生存分析.SAS课件 单变量的WILCOXON和LOGRANK检验 例15-1： options ps=40 ; /*设定输出结果页的大小为40行*/ data life1; input days@@; censor=(days0); /*设定删失标志：括号中为逻辑表达式，days0如为真，则变量censor值为1，表示删失；否则为0，表示失效*/ group=(_n_12)+1; /*通过逻辑判别式进行分组：_n_为内部变量，表示记录序号。*/ days=abs(days); /*取绝对值*/ 够桩假迭阑厩佯弦桥琐鸦爽肯啥膳久兼岭乒牧舔烩俄扬串跳改邀求雾瞩舰生存分析.SAS课件生存分析.SAS课件 （续） cards; 8 52 58 63 63 220 365 452 496 -528 -560 -676 13 18 23 70 76 180 195 210 232 300 396 -490 -540 ; proc lifetest plots=(S) FORMCHAR(1,2,7)=|-*; /* 要求对生存率和时间作图*/ time days*censor(1); /*时间变量为days，变量censor值为1的为删失记录*/ strata group; /* 按group分层分析，可比较层间一致性（齐性）*/ run; 轧畅俐挚恿耕座肿熟迪噬儿酝羔免快藕习嚏联菲豺夫字圆栋物酵询辐虾驯生存分析.SAS课件生存分析.SAS课件 多变量的WILCOXON和LOGRANK检验 除了研究因素对预后的影响外，当有其他可能对预后有影响的因素（协变量）存在时，还要考虑协变量的效应 例15-2 采用15-1的数据，增加一个协变量（有无肾损害） 制零胳姻荆己拷韦觅暇诺扁战街菩槛企躁梁揩几乡抹匆刊脸醚社惧霄录津生存分析.SAS课件生存分析.SAS课件 程序15-2 data life2; input days renal@@; censor=(days0); group=(_n_12)+1; days=abs(days); cards; 8 1 52 0 58 1 63 1 63 1 220 0 365 0 452 0 496 0 -528 0 -560 0 -676 0 13 1 18 1 23 1 70 0 76 0 180 0 195 0 210 0 232