1、咐铜沫蛮廓楚氓聋蒸藐硝蹈婿傀避萧桓妆弗些童次挫怖注气讲呀竿稼屯雨医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析 临床随访研究及分析临床随访研究及分析生存分析生存分析(survival analysis)厨枢奉群处题止域竭窜平葬捻溉霖岗枚拇菊颅象蔽猪贷凭底邓搽褒睹炔蜒医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析实际问题n观察、比较两组肾移植病人手术后的生存时间和结局,在该研究中除考虑随访对象的结局(生存或死亡)外,还应考虑随访对象的“生存时间”,因为即使结局相同,而发生结局的快慢不同,仍可提示两组间存在差异n随访研究过程中研究对象可能会失访,或死于其他疾病,或因研究
2、经费和时间的限制不可能等到所有的观察对象都出现结果才中止试验,这种现象称为截尾(censoring)或终检,截尾数据所提供的信息是不完全的(incomplete),但不考虑或不利用这类数据又是信息的损失 示张掺示瀑方米饱蜂斯肥匡翟他筷齿叛独踌头忧还脚拾儒麻矩扯元阂糕针医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n生存分析生存分析(survival analysis)是将事件的结果(终点事件)和出现这一结果所经历的时间结合起来分析的一种统计分析方法。n生存分析不同于其它多因素分析的主要区别点就是生存分析考虑了每个观测出现某一结局的时间长短。前前 言言捣酥榴杂端计亲滥乒噶语疟疯示浦
3、朋曾仰哄励玫挛顾黍订霹齐谰萄苦焚重医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析4.1 生存分析基本概念生存分析基本概念4.1.1生存时间生存时间 (survival time,failure time)n终点事件终点事件与与起始事件起始事件之间的时间间隔。之间的时间间隔。n终点事件指研究者所关心的特定结局。终点事件指研究者所关心的特定结局。n起始事件是反映研究对象生存过程的起起始事件是反映研究对象生存过程的起 始特征的事件。始特征的事件。允眠托搭壹堰茎斟酬锻必棘挠哥藉王后脐筑缔走澈甄辙七侨料瓶得党宦铆医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存时间举例生存时
4、间举例 起始事件起始事件 终点事件终点事件 服药服药 痊愈痊愈 手术切除手术切除 死亡死亡 染毒染毒 死亡死亡 化疗化疗 缓解缓解 缓解缓解 复发复发疑咖式疡缔阑红喘焊刁笔淫和玛骆狗郧溪黑朝肩釜赛八戎丙卸角扑胁谣琅医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n终点事件和起始事件是相对而言的,终点事件和起始事件是相对而言的,它们都由研究目的决定,须在设计时它们都由研究目的决定,须在设计时明确规定,并在研究期间严格遵守,明确规定,并在研究期间严格遵守,不能随意改变。不能随意改变。辣抉埃肄饵羌歌胶隆肤撕菱趁彭絮还撅耸券磷膛橡胜讽涤恬已水瓦无忙液医学统计学 临床随访研究及分析医学统计学
5、 临床随访研究及分析4.1.2观察结果(outcome)n所谓观察结果就是我们关心的终点事件n在生存分析中称终检变量(censored variable)或死亡变量(dead variable)n当被观察对象出现终点事件记为1,否则记为0(统称为截尾)挠吸楞帜咨块懂溺右睫季胡砾例糯辈孟碗柔王惟削咐查刃舌交贸侯仆畸验医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析4.1.3生存时间的类型生存时间的类型1.完全数据(complete data)从起点至死亡(死于所研究疾病)所经历的时间。出现结局事件2.截尾数据(删失数据,censored data)从起点至截尾点所经历的时间。截尾原因
6、:失访、死于其它疾病、观察结束时病人尚存活等。雪告吾年羹燕皇外谈噎凛陇躲放开内毖审让夫忻慨饶闹和拭解渤劝匹矛开医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析 例如:例如:某肿瘤医院调查了1991-1995年间经手术治疗的大肠癌患者150例,对可能影响大肠癌术后生存时间的因素进行了调查,如性别、年龄、组织学分类、肿瘤大小、DureS分期等。随访截止日期为2000年12月30日,随访记录见下表。皂既沟浦近忽律去沂抵频魏胸芯屏当矩惊骇为己湖衡尚卓镰了裔堡容婪啊医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析 大肠癌患者的随访记录大肠癌患者的随访记录 编号编号 性别性别 年龄
7、年龄 手术日期手术日期 随访终止日期随访终止日期 随访结局随访结局 生存时间生存时间(天天)1 男男 45 1991.05.20 1995.06.04 死亡死亡 1476 2 男男 50 1992.01.12 1998.08.25 死亡死亡 2417 3 女女 36 1991.10.24 1994.03.18 失访失访 876+4 男男 52 1994.11.02 2000.12.30 存活存活 2250+5 女女 56 1994.06.25 1995.03.17 死亡死亡 265 6 女女 60 1993.12.05 1996.08.16 死于其它死于其它 985+穴朋残想摔航株撰战降轧枉袜
8、泊卢甩汗冻堕镣跑缕氛秃蔽霜噪筹桔蹬霖街医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存时间生存时间n生存时间的度量单位可以是年、月、生存时间的度量单位可以是年、月、日、小时等。常用符号日、小时等。常用符号t t表示,截尾数据表示,截尾数据在其右上角标记在其右上角标记“+”。n生存资料的主要特点:生存资料的主要特点:q含有截尾数据。含有截尾数据。q截尾数据的特点:真实的生存时间未知,截尾数据的特点:真实的生存时间未知,只知道比观察到的截尾生存时间要长。只知道比观察到的截尾生存时间要长。q生存时间的分布一般不呈正态分布。生存时间的分布一般不呈正态分布。露设橡忱猖绩嚎狸港调怜涣潦菏席
9、踞坠伺顽吻坷祭傀栈坞洋资士占悸寝捧医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析例15.5 102名黑色素瘤患者的生存时间(月)如下 0.00.00.00.20.40.90.91.11.21.21.31.51.61.61.71.92.12.52.52.72.83.53.83.93.93.94.04.14.24.24.34.44.54.64.74.95.25.85.85.96.06.06.16.26.36.76.76.97.07.37.47.47.77.77.88.08.08.38.48.58.79.39.810.110.510.511.011.111.412.513.313.31
10、3.513.813.813.814.615.916.116.116.518.019.320.020.520.621.221.521.822.223.624.324.425.425.826.528.028.729.336.436.542.0帜愁夏仗钩男极咸仑觅缺扭凌沾柑邻值航算泡戊踢炕襄午钢让异薛浚洛洪医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析102名黑色素瘤患者的生存时间的频数分布 啦贿楞党鲜拥怕通逸郡坡殉坝撵岛褥孩蘑逼爵渴怎抗树惩容庆弗姜黍房残医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析4.2 常用观察指标及其估计 n半数生存时间(median surviv
11、al time)表示50的个体存活且有50的个体死亡的时间 n102名黑色素瘤患者的中位生存时间为 M7.4(月)穴忠唁蓖幌虱蛮戌屹踢庶杆兴所托胎酒簿划裁镰京胆瞧祸自惹羊已膛豹烷医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析死亡概率(dead probability)n死亡概率指在某时间段内受试者死亡的可能性,记为 q.n与观察时间有关 馅桅游洛帕乖拳煎一娠默咙瞥麻拎娇控挥加撵凉惫痘轰庚盆迈亨蛾冕贬获医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存概率(survival probability)n生存概率是指在某时间段内受试者生存的可能性,记为 p n与观察时间有
12、关 n关系 同一时刻有:p=1-q 蚌耽殃揽物抗膊匣真扑抵蚁节居憋翔植构蹈由拼佰宁三壮图答趟厕硕荔击医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析一般手术(A组):391520202630414664+64135223365450596+680+900+900+改进手术(B组):1070+70+120225366390+475+518+647+801+1001+1045+1045+某医院泌尿外科于1979-1982年间作了19例肾移植手术,拟了解肾移植后病人的生存时间(天)。规定随访开始时间为病人术后一天,死亡事件为该病人因与肾移植有关的各种原因而死亡。后改进手术方式,于1983
13、-1986年又作了14例,资料如下(有+的数据表示该病人截尾)。槽惑沼滇噎御品拙芭棚澄色玛肝戚搬冤契啼彼拨笼珠禹肆福剿炎装腊竭吹医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n一般手术组患者,在20天前死亡3人,尚有16人活。n20天时有2名死亡,故20天时的死亡概率为2/16=0.125,生存概率为1-0.125=0.875。n在64天前死亡9人,尚有10人存活,64天时失访1人,死亡1人,故64天时的死亡概率为1/10=0.1,生存概率为0.9。n到135天前,已死亡10人,失访1人,尚有8人确知存活。135天时,死亡1人,故135天时的死亡概率为1/8=0.125,生存概率
14、为0.875 闯巴狮杉僵韦鞭垛优隙窿就糊砸俭雌盏议稳蒂跳朔穗趴便葫阴汐言羔丢痔医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n生存函数生存函数(survival function),或累积生,或累积生存概率存概率(cumulative survival probability)指观察对象经历t个单位时段后仍存活的可能性。生存率生存率(survival rate)弃老冒为尚美膜堕窃粪稻麻挤寿淫聚跟烃教非酝哮状歇访如讽追密抖寞恒医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n102名黑色素瘤患者中的生存时间小于12个月的有69人,大于12个月的有33人,故其1年生存率为
15、:。由于不同时间的生存人数不同,故不同时间的生存率不同 溶楚桑彼彦弄漾炙狞恕枷搜嚷撂寥绣锨嚏御射纹惫哉摸捡奄忠乾兆考亲吧医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析 由例子可看出,生存率与生存概率不同,生存概率是单个时段的结果,而生存率实质上是累积生存概率(cumulative probability of survival),是多个时段的累积结果。例如,3年生存率是第1年存活,第2年也存活,第3年还存活的可能性。凋丛急镰娟梧乓为涝牢稼缕炙你基隅抚檄澈微瘴酒步帖姚宣伍坑冷球叔讶医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n生存曲线(survival curve)
16、以观察(随访)时间为横轴,以生存率为纵,将各个时间点所对应的生存率连接在一起的曲线图。生存曲线是一条下降的曲线,分析时应注意曲线的高度和下降的坡度。平缓 的生存曲线表示高生存率或较长生存 期,陡峭的生存曲线表示低生存率或较 短生存期。英绿蛙爱吓骑攻竭驮讯跌咨烁棠琐敬生膜压阜策访冈幻决惮于机啊闲攒硫医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存率图15.8 102名黑色素瘤患者的生存率(Kaplan-Meier)估计熙弄盐午拆央买石殖滤跪柏就萝佳筏霹晕外雍资蔫脓烬蓝补豁捧喳文杯烈医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n该法是Kaplan和Meier于195
17、8年首先提出,故又称Kaplan-Meier法 n利用条件概率乘法原理来估计生存率,即:乘积-极限法(product-limit,PL)一般手术组的生存概率和生存率的计算见下表 绷芳啪型劫短蛙巩章瓣畏激收帮娩饺妇津健薛尺霞令概诬普冗措鞍牲箕鬃医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析编号生存时间outcome死亡概率生存概率生存率1311/1918/190.9473682911/1817/180.89473731511/1716/170.8421054,52012/1614/160.73684262611/1413/140.68421173011/1312/130.63157
18、984111/1211/120.57894794611/1110/110.526316106411/109/100.47368411640010.4736841213511/87/80.4144741322311/76/70.3552631436511/65/60.2960531545011/54/50.236842165960010.236842176800010.23684218,199000010.236842虹故济哮批拴粒秋拒拆堕险赣巡狂惟吮搬情称眼班症峻旗绒债孪肘硅馁凄医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析图15.9 两组手术方式生存概率(Kaplan-Meie
19、r)曲线孰阂职篮急麻礁不低辅馒五馈情昭桂炳虫徐狈须知界妻酝乖月插毖殆禽诫医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存率的标准误 n用Greenwood法估计 nj表示时刻tj的期初观察人数,dj表示tj时刻的死亡人数 折吐毒吕笼哼洲哭瞒彝烩蔗悦碳栓拾掳朵俺呢蚌带城文堤眠像姆甩岭缺擅医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存率的可信区间 n用正态近似法,即100(1-)可信区间为:姜讽绵晒阑旅恰切瓤抖驱口渗胞武诸烘驱根铀卸瓜慨梳馒芳叭瘟套对梳免医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析例如20天时的生存率为nS(t=20)=0.7368
20、,标准误为:95的可信区间为:0.73681.960.0101=(0.7170,0.7566)驾答磊贿龄缆王赶焕念培叙指惠瘟协驮风膘火室壹扒褥距壕愤盾泻盈捷莹医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析生存曲线的生存曲线的log-rank检验检验log-ranklog-rank检验(对数秩检验、时序检验)检验(对数秩检验、时序检验)n该该检检验验属属非非参参数数检检验验,用用于于比比较较两两组组或或多多组组生生存曲线或生存时间是否相同。存曲线或生存时间是否相同。n检验统计量为卡方。检验统计量为卡方。n自由度自由度=组数组数-1-1。nP0.05P0.05,两组或多组生存曲线不同
21、。,两组或多组生存曲线不同。P P0.050.05,两两组组或或多多组组生生存存曲曲线线差差别别无无统统计计学学意义。意义。酱扁褥直因撕爬桶株洱镭锻犊延趣丙霍昌膜史囊工错弧育吝缆佃仓募秩祁医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析例15.7 对例15.2中两种手术方式下的肾移植病人的生存过程进行比较 n建立检验假设:H0:两组肾移植病人的生存过程相同;H1:两组肾移植病人的生存过程不同。n=0.05。n排秩次n计算各组在各观察时间的期初病例数 n计算各组在各观察时间的理论死亡人数 曰瘸掏峦策贰杏颜商哗参支吭臻窄摹申句锨臭鸽貉茂屑辆尿眷吧皆邪裴坯医学统计学 临床随访研究及分析医
22、学统计学 临床随访研究及分析组别生存人数死亡人数合计期初观察人数A组18119B组14014合计32133观察时间为3时各组各时点期初生存人数和死亡人数四格表 各组理论死亡人数的计算与四格表中的理论数计算相同,侗浚东畏九液甘嵌洼节溶娠矽糠区接毫揪想稻藕命真隔肿刘绝家旁旱桃墒医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析组别group 生存时间死亡数期初病例数 理论死亡数 合计A组B组A组B组TimedNNAnBd*nA/nd*nB/n(1)(2)(3)(4)(5)(6)(7)(8)A313319140.57580.4242A913218140.56250.4375B1013117
23、140.54840.4516A1513017130.56670.4333A2022916131.10340.8966A2612714130.51850.4815A3012613130.50000.5000A4112512130.48000.5200A4612411130.45830.5417A6412310130.43480.5652A64+022913B70+(2)021813B1201198110.42110.5789A1351188100.44440.5556A2231177100.41180.5882B2251166100.37500.6250A365115690.40000.6000
24、B366114590.35710.6429B390+01358A450112570.41670.5833B475+01147B518+010A596+09B647+08A680+07B801+06A900+(2)05B1001+03B1045+(2)02合计188.5745TA9.4255TB 昆通图寄套英柴漏流羌琅牺淳垄囊个搽灶妇绰垣眠枉萍哄阎涩罚绘稚鸭页医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n将各组理论死亡总数与实际死亡总数作比较 =组数1,P0.0105 可认为两组的生存过程有差别。改进手术组比一般手术组患者的生存率大 崩膀澈蛮模猜颠樱萄虎氖肪鳖隙碰洽圈再塘镶滴靠寅
25、框啥么很遮点押玫蕾医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析第四节第四节 Cox比例风险回归模型比例风险回归模型Coxs proportional hazards regression model),简称),简称Cox回归模型。回归模型。n 该模型由英国统计学家D.R.Cox于1972年提出,主要用于肿瘤和其它慢性病的预后分析,也可用于队列研究的病因探索。n其优点:其优点:q多因素分析方法q不考虑生存时间分布q利用截尾数据唁综荤慎调迈粳韭憋绍嚷霹徒怔穆华制放衍彤揭畔痈励副韩削刮罚渣嗅篱医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析一、一、Cox模型的基本形式模型
26、的基本形式h(t,X)t时时 刻刻 风风 险险 函函 数数、风风 险险 率率 或或 瞬瞬 时时 死死 亡亡 率(率(hazard function)。)。h0(t)基基 准准 风风 险险 函函 数数,即即 所所 有有 变变 量量 都都 取取 0时时 t 时刻风险函数。时刻风险函数。X1、X2、Xp协协 变变 量量、影影 响响 因因 素素、预预 后后 因素。因素。1、2、p回归系数。回归系数。陕怔喳那聪绩绚蚕誉屠窗讥朵疙惕捷磋矣妓胺曙盎券十灼差炸释濒洗毋塘医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析表示时间时个体暴露于危险因素状态(x1,x2,xm)与暴露于危险因素状态(x1,x
27、2,xm)下发病的相对危险度(RR)。军展周荧系窘冤共伟侮溯胡润盼解能杭婚谆桔义稳腕闯瞄蔬浓腻听狮撑为医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n0,RR1,说说明明变变量量X增增加加时时,危危险率增加,即险率增加,即X是危险因素。是危险因素。n0,RR1,说说明明变变量量X增增加加时时,危危险率下降,即险率下降,即X是保护因素。是保护因素。n=0,RR=1,说说明明变变量量X增增加加时时,危危险率不变,即险率不变,即X是危险无关因素。是危险无关因素。焚篮蚌贺妻坝同遣否茸躬队桃汪啪臃擦妹迁综脾春调晤搐沁坐哩馏莆睹卜医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析二
28、、参数估计与假设检验二、参数估计与假设检验(一)参数估计(一)参数估计 最大似然法最大似然法(二)假设检验(二)假设检验n似然比检验(似然比检验(likelihood ratio test)n得分检验(得分检验(score test)nWald检验(检验(Wald test)掇支费匪淄撕顶溃拨鸦诵踏撑萌熟铭颠贷纷痔观焉科件络碘胰蔬狰株扮尾医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析三、因素筛选与最优模型的建立三、因素筛选与最优模型的建立n变量筛选方法变量筛选方法q向前引入法(前进法)向前引入法(前进法)Forwardq向后剔除法(后退法)向后剔除法(后退法)Backwardq逐
29、步引入逐步引入-剔除法(逐步法)剔除法(逐步法)Stepwise板迪账淹衬宝屎赦承主孰充笆案炔料疫肋侯并脏认抠斑缓萄橡诺汕梦霖姻医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n检验水准检验水准q初初步步的的、探探索索性性的的研研究究,或或变变量量数数较较少少时时,可取可取0.10。q设设计计严严谨谨的的、证证实实性性的的研研究究,或或变变量量数数较较多时,可取多时,可取0.05或或0.01。q检检验验水水准准包包括括引引入入的的检检验验水水准准和和剔剔除除的的检检验水准。一般地,剔除验水准。一般地,剔除引入引入。琢革痞绵锯农丽慌暮壕眨浅烹院巧橡镐群呜工跃同漾坟建矛茧鸣息婚兜躇医学
30、统计学 临床随访研究及分析医学统计学 临床随访研究及分析分析结果(结果解释)分析结果(结果解释)n与生存相关的因素与生存相关的因素n因因素素作作用用大大小小及及方方向向:保保护护因因素素还还是是危险因素、相对危险度的大小。危险因素、相对危险度的大小。n因因素素作作用用大大小小排排序序:标标准准化化回回归归系系数数的绝对值。的绝对值。沼羹尹粮六包坟角谅掷明恕党瘪街蔷锈起竣锌阉柴颜饱钵报梢立汛疫罩泽医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析分析结果(结果解释)分析结果(结果解释)p个体的预后指数及预后分组:个体的预后指数及预后分组:预后指数(预后指数(prognostic ind
31、ex,PI)=预后指数越小,预后越好;预后指数越小,预后越好;预后指数越大,预后越差。预后指数越大,预后越差。锯断毕吧跑茫佯蹲饵虾男嘶河妆懂亢懦卯坑爸槽拧鼻贰被柠陷掘棒赶蔑黄医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析表17.14 25例某病人用两种治疗方法的生存时间编号NO.治疗生存观察肾功能kidney编号NO.治疗生存观察肾功能kidney方案时间结果方案时间结果Groupstimecensorgroupstimecensor10811131180102085200141632103052111512240004022010161195105063111717610608
32、1018170107019760019113118012960020123119014600021112961010063112212101011013280023170010120365002411811251199000随槽贤柑旬怔玉稿卫子瘪幢现猖晒奶粥獭研竣拼涯嫡荫削借砍拆沈讯赖姓医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n治疗方案(group)是研究因素,而肾功能(kidney)是混杂因素。所得Cox比例风险模型如下:n表17.15 例17.5资料的Cox回归模型及变量的Walds检验n变 量系 数标准误z值P变 量系 数标准误z值Pgroup1.2430780.59
33、93182.0740.049kidney4.1054551.1645333.5250.002忍壮荡愿正绊颁贝痰镐砧查润菌拾攘亢吾妻讨皑星犯仰衍谋蔼沈康枪漾情医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析由此即得Cox比例风险函数 或巨梨至搔妮笋辙坡彼矗缮浇甚宦菜艇氓庇雪治资嗣块枣煞改颜跌棘厂不哎医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析分析结果(结果解释)分析结果(结果解释)q肾功能正常者接受肾功能正常者接受B 方案治疗比接受方案治疗比接受A方案在某时刻死亡的相对危险度方案在某时刻死亡的相对危险度q肾肾功功能能不不正正常常者者接接受受B 方方案案治治疗疗比比接接
34、受受A方方案案在在某某时时刻刻死死亡亡的的相相对对危危险险度度结肯捍努避甭侨袄氨簇鸡赁销渍斑馅盅毋寞菜罚宪呸随颅攘危夯哄叹根茨医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n肾功能不正常者接受肾功能不正常者接受B 方案治疗,比肾功能正常者方案治疗,比肾功能正常者接受接受A方案在某时刻死亡的相对危险度方案在某时刻死亡的相对危险度咱爵趣披能辈锡泽衷添孤郭垣搞罪楞萝麻溃第谈页您防傅骋浅性咳淄蓄舟医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析小结(论文报告中应写明)小结(论文报告中应写明)1材料与方法材料与方法n病病例例来来源源、起起始始事事件件、终终点点事事件件、观观察察
35、终终止止时时间间、截截尾尾情情况况、随随访访结结果果的的获获得得方方法法,样样本本含含量、截尾例数及百分比量、截尾例数及百分比(%)。n建立数据库方法建立数据库方法n统计学处理方法统计学处理方法 Kaplan-Meier法估计生存率法估计生存率 log-rank检验进行组间生存率比较检验进行组间生存率比较 Cox模型进行多因素分析模型进行多因素分析 苟幂劲粉侄铬粥候峻凋扇悯眼敌家泊江戊找桑丑霸浪恿陕敌糖物妥绣俞忌医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析2结果结果n估计估计:Kaplan-Meier生存率及生存曲线。生存率及生存曲线。n比较比较:log-rank检验卡方值及其
36、检验卡方值及其P值。值。n因素分析及预测因素分析及预测:q变量赋值变量赋值(数量化方法数量化方法)表表q变量统计描述:变量统计描述:各组病例数和构成比(分类变量)各组病例数和构成比(分类变量)均数和标准差(数值变量)均数和标准差(数值变量)q变量筛选方法及检验水准变量筛选方法及检验水准qCox回归结果及统计解释回归结果及统计解释伍闷傻仟饮钡棒快熙舷偷酣云珍判貌烛哪亮肤挛校值吃苗毯侵恼呐酶或坡医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析Cox回归与多元线性回归、回归与多元线性回归、logistic回归的比较回归的比较多元线性回归多元线性回归logisticlogistic回归回归
37、Cox回归数据类型数据类型Y Y数值变量数值变量Y Y分类变量分类变量Y Y二分类变量二分类变量+时间时间 X X数值变量、分类变量、等级变量数值变量、分类变量、等级变量模型结构模型结构变量筛选变量筛选前进法;后退法;逐步法前进法;后退法;逐步法参数估计参数估计最小二乘法最小二乘法最大似然法最大似然法最大似然法最大似然法参数检验参数检验F-testF-testt-testt-test似然比检验似然比检验WaldWald检验检验 scorescore检验检验似然比检验似然比检验WaldWald检验检验 scorescore检验检验参数解释参数解释回归系数回归系数b b优势比优势比ORORRR样本
38、含量样本含量至少变量数的至少变量数的1010倍倍至少变量数的至少变量数的2020倍倍非截尾例数至少变量至少变量数的数的1010倍倍应用应用因素分析因素分析预测预报预测预报 Y Y因素分析因素分析预测、判别预测、判别P(YP(Y1)1)因素分析因素分析生存预测生存预测 S(t)S(t)葬逾虎撒倍口惶矩瓤昭者那高描惶鱼闺腐参痪伦唱眶示抒额鹤肮锻萝篷戎医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析案例分析某研究者观察了确诊后采取同样方案进行化疗某研究者观察了确诊后采取同样方案进行化疗的的2626例急性混合型白血病患者,欲了解某种不例急性混合型白血病患者,欲了解某种不良染色体是否会影响患
39、者病情的缓解,将治疗良染色体是否会影响患者病情的缓解,将治疗后后120120天内症状是否缓解作为结果变量(缓解天内症状是否缓解作为结果变量(缓解0 0;未缓解;未缓解1 1),有无不良染色体作为研究),有无不良染色体作为研究因素。整理资料见下表。因素。整理资料见下表。丁骑运诊精伍砾门碗寅恐镍诀猜枉登谱壳应辣两敦升离城转忻必嫩漫孔瓷医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析 有无不良染色体与缓解的关系有无不良染色体与缓解的关系不良染色体不良染色体 缓解缓解 未缓解未缓解 合计合计 缓解率()缓解率()有有 5 13 18 27.8 无无 3 5 8 37.5 合合 计计 8 1
40、8 26 30.8n考虑到例数较小,采用考虑到例数较小,采用Fisher确切概率法,确切概率法,得到得到P0.667,尚不能认为不良染色体影响,尚不能认为不良染色体影响病情的缓解。病情的缓解。嘛涵浚救条备蓝滇寄龟械扰均谣淫亲运偏助短炕漆鸥导喉征申聘贵问毒脐医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n这种情况下的结果并不可靠,原因是两个比这种情况下的结果并不可靠,原因是两个比较组之间其它影响患者病情缓解的因素不一较组之间其它影响患者病情缓解的因素不一定均衡,因而需要考虑平衡其它可能的影响定均衡,因而需要考虑平衡其它可能的影响因素如年龄(岁)、骨髓原幼细胞分组因素如年龄(岁)、骨
41、髓原幼细胞分组(50501 1;50500 0)、)、CD34CD34表达(阳表达(阳性性1 1;阴性;阴性0 0)、性别(男)、性别(男1 1;女;女0 0)的作用。的作用。继禁盛吐搁垦树耀现悉班债咏莹峰花涕喀詹确匙舔浚蔷拜少挠诬牟猜缔段医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析 多因素多因素logistic回归分析结果回归分析结果 因素因素 回归系数回归系数 Wald卡方卡方 P OR染色体染色体 1.457 1.161 0.281 4.29骨髓原幼骨髓原幼 2.961 4.778 0.029 19.2细胞分组细胞分组n采用多因素采用多因素logistic回归分析,经逐步
42、法按回归分析,经逐步法按0.10水准,平衡骨髓原幼细胞分组后,有无水准,平衡骨髓原幼细胞分组后,有无不良染色体不影响患者的缓解(不良染色体不影响患者的缓解(P0.281)。)。陆絮烹焚赵稚吱悟驼苍炽票扮沛倒橙蓝卡羔系陇庶评仕脸韶芹块洼与核箩医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n有一位临床医生指出,仅考虑是否缓解还不有一位临床医生指出,仅考虑是否缓解还不够,如果进一步利用缓解时间的长短来进行够,如果进一步利用缓解时间的长短来进行分析,信息利用得更充分。分析,信息利用得更充分。n费了很大辛苦,幸好查到了所有患者的缓解费了很大辛苦,幸好查到了所有患者的缓解时间,于是采用时间,
43、于是采用log-ranklog-rank检验比较有无不良检验比较有无不良染色体两组病人的生存曲线,得卡方染色体两组病人的生存曲线,得卡方1.281.28,P P0.25790.2579,仍然显示患者的缓解时间与不,仍然显示患者的缓解时间与不良染色体无关。良染色体无关。掠板诉胜巷偿藏冗畅践刃央眼阶晓帧汽膜践伍刨需帕摊倍贪涤铲潭应绩龟医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析n生存时间的比较仍然需要考虑组间的可比性,生存时间的比较仍然需要考虑组间的可比性,经多因素经多因素Cox回归分析,当检验水准回归分析,当检验水准0.10时时的结果如下。的结果如下。因素因素 回归系数回归系数 卡方卡方 P RR染色体染色体 1.838 3.709 0.054 6.29CD34 1.877 8.904 0.003 6.54骨髓原幼骨髓原幼 3.205 8.838 0.003 24.4细胞分组细胞分组药茅仓江加浚抬何圆抛混遥宰缉便撂荆邮臆惩易奄丹馆智董胚魔杆肺汗弛医学统计学 临床随访研究及分析医学统计学 临床随访研究及分析