欢迎来到天天文库
浏览记录
ID:21982566
大小:234.50 KB
页数:16页
时间:2018-10-26
《第十四章生存分析的sas实现》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、第十四章生存分析的SAS实现例14-2McKelveyetal(1976)收集了3期的某型淋巴瘤患者的生存时间(天)。分别是6,19,32,42,42,43+,94,126+,169+,207,211+,227+,253,255+,270+,310+,316+,335+,346+。SAS分析程序datach14_2;定义数据集名inputtc@@;输入生存时间变量t和截尾情况变量ccards;开始输入数据611913214214214309411260169020712110227025312550270031003160
2、33503460;proclifetestplots=(s);利用lifetest过程进行生存分析并作生存函数图timet*c(0);指定时间变量和截尾变量并指出数据截尾时截尾变量的取值run;SAS软件输出结果Product-LimitSurvivalEstimatesSurvivalStandardNumberNumbertSurvivalFailureErrorFailedLeft0.0001.0000000196.0000.94740.05260.051211819.0000.89470.10530.0704217
3、32.0000.84210.15790.083731642.000...41542.0000.73680.26320.101051443.000*...51394.0000.68020.31980.1080612126.000*...611169.000*...610207.0000.61210.38790.116779211.000*...78227.000*...77253.0000.52470.47530.128786255.000*...85270.000*...84310.000*...83316.000*...
4、82335.000*...81346.000*...80NOTE:Themarkedsurvivaltimesarecensoredobservations.SummaryStatisticsforTimeVariabletQuartileEstimatesPoint95%ConfidenceIntervalPercentEstimate[LowerUpper)75...50.94.000.2542.00032.000.MeanStandardError181.70124.497NOTE:Themeansurvivalti
5、meanditsstandarderrorwereunderestimatedbecausethelargestobservationwascensoredandtheestimationwasrestrictedtothelargesteventtimeSummaryoftheNumberofCensoredandUncensoredValuesPercentTotalFailedCensoredCensored1981157.89SAS软件输出结果解释该结果包含四个部分:第一部分用乘积极限法估计了生存率(Surviva
6、l),死亡率(Failure),生存率的标准误(SurvivalStandardError),死亡例数(NumberFailed)和该时间点前的生存例数(NumberLeft)。其中带有*号的表示截尾;第二部分给出了关于生存时间的描述性统计量,包括75%,50%和25%分位数以及相应的95%的可信区间(95%ConfidenceInterval),还有均数(Mean)和标准误(StandardError)从结果可以看出平均生存时间为181.701天;第三部分列出了完全数据(Failed),截尾数据(Censored)的例数
7、,以及截尾数据占全部数据的百分比(PercentCensored)。最后是生存曲线图。教材中的说明现用Kaplan-Meier法计算生存率,步骤如下:(1)将所有生存时间按从小到大顺序排列(见表14-2第(2)列)并标上序号(第(1)列)。(2)列出各t时点前的存活病例数n(第(3)列)、各个时间点的死亡人数d(第(4)列)和截尾人数c(第(5)列)。(3)计算各t时刻的死亡概率(第(6)列)。例如生存时间为32天时,死亡概率为。(4)计算各t时刻的生存概率(第(7)列)。例如生存时间为32天时,生存概率为。(5)计算各t
8、时刻的生存率(第(8)列)。例如生存时间为32天时,生存率为,由此验证了在没有截尾数据的情况下,式(14-4)与式(14-5)是相同的。(6)以时间t为横指标,生存率为纵指标,作生存曲线图(图14-1)。表14-2Kaplan-Meier法计算生存率的计算用表 (1)(2)(3)(4)(5)(6)(
此文档下载收益归作者所有