第十四章 生存分析的Stata实现

《第十四章 生存分析的Stata实现》由会员分享，可在线阅读，更多相关《第十四章 生存分析的Stata实现（8页珍藏版）》请在装配图网上搜索。

1、第十四章 生存分析的Stata实现本章使用的STATA命令:设置生存分析数据结构stset生存时间变量,f（结局变量）作图（先执行stset命令）单个生存曲线sts graph多组生存曲线sts graph ,by（分组变量）生存率（先执行stset命令）单组生存率sts list多组生存率sts list, by（分组变量）Logrank 检验（先执行stset命令）sts test分组变量寿命表单组寿命表ltable生存时间结局变量，graph多组寿命表ltable生存时间 结局变量,by（分组变量）graph结局变量为 1 表示失效事件发生例 14-2 McKelvey et al（19

2、76） 收集了 3 期的某型淋巴瘤患者的生存时间（天）。分别是 6，19，32，42，42，43+，94，126+，169+，207，211+，227+，253，255+，270+，310+，316+，335+，346+。现用 Kaplan-Meier 法计算生存率。Stata 数据为:timed61191321421421430941126016902071211022702531255027003100316033503460Stata 命令为：stset time,failure(d) sts liststs graph结果为：failure _d: d = 1 analysis tim

3、e _t: timeTimeBeg.TotalNetSurvivorFunctionSt d.FailLos tError95% Conf.Int.619100.94740.05120.68120.99241918100.89470.07040.64080.97263217100.84210.08370.58650.94624216200.73680.10100.47890.88104314010.73680.10100.47890.88109413100.68020.10800.42140.842112612010.68020.10800.42140.842116911010.68020.1

4、0800.42140.842120710100.61210.11670.34970.79512119010.61210.11670.34970.79512278010.61210.11670.34970.79512537100.52470.12870.25700.73632556010.52470.12870.25700.73632705010.52470.12870.25700.73633104010.52470.12870.25700.73633163010.52470.12870.25700.73633352010.52470.12870.25700.73633461010.52470.

5、12870.25700.736357 AuoooLnNO00.0Kaplan-Meier survival estimate例 14 3 下面是来自于 Berkson & Gage(1950) 的一个研究队列。为了叙述 方便，把原来的出院后的生存时间改称为某恶性肿瘤术后生存时间。共有 374 名患者进入研究队列。表 14 3寿命表法计算生存率的计算用表(1)(2)(3)(4)(5)(6)(7)(8)(9)序号术后生存期初观察期内期内截尾校正期初死亡生存生存率年数例数死亡人数人数概率概率tndcn =n-c/2q=d/ncpi-qS(t)103749003740.24060.75940.7594

6、212847602840.26760.73240.5561322085102080.24520.75480.41984315725121510.16560.83440.350354120205117.50.17020.82980.290765957990.50.07730.92270.268276794974.50.05370.94630.253887661364.50.01550.98450.249898623559.50.05040.94960.2372109542551.50.03880.96120.22801110472126340.61760.38240.0872Stata 数据格式

7、为：timewd190127613511425152016717418119311021112111002003004120550690790830950105011260Stata 命令为：stset time weight=w,failure(d=1) ltable time d weight=w,graph结果为：IntervalBeg.TotalDeathsSt d.Los tSurvivalError95% Conf.Int.123749000.75940.02210.71270.7995232847600.55610.02570.50430.6048342085100.41980.

8、02550.36950.46924515725120.35030.02480.30200.3989561202050.29070.02390.24470.33816795790.26820.02350.22320.31507879490.25380.02330.20940.30048966130.24980.02330.20550.296591062350.23720.02320.19320.2839101154250.22800.02320.18420.274811124721260.08720.02100.05180.1339例 14 4 在例 14 3 中介绍了 McKelvey et

9、al 收集了 19 名 3 期淋巴瘤病 人的生存资料，该研究者还同时收集了 4 期淋巴瘤病人的生存数据共 61 例：4， 6，10，11，11，11，13，17，20，20，21，22，24，24，29，30，30，31，33 34，35，39，40，41+，43+，45，46，50，56，61+，61+，63，68，82，85，88， 89，90，93，104，110，134，137，160+，169，171，173，175，184，201，222 235+，247+，260+，284+，290+，291+，302+，304+，341+，345+。要比较 3 期淋巴 瘤患者和4 期淋巴瘤患者

10、的生存率是否有差异。Stata 数据格式为：timedgroup61119113211421142114301941112601169012071121101881289129012931210412110121341213712160021691217112173121751218412201122221223502247022600228402290022910230202304023410234502Stata 命令为：stset time,failure(d=1)sts test groupsts graph,by(group)结果：Log-rank test for equality of survivor functions|EventsEventsgroup |1observedexpected+1|2|184616.6937.31+Total |5454.00chi2(l) =6.71Prchi2 =0.0096P=0.0096，3 期淋巴瘤患者和 4期淋巴瘤患者的生存率有差别。QU157 AuLnNO00.0