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

《第十四章生存分析的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;开始输入数据6119132142142143094112601690207121102270253125502

2、7003100316033503460;proclifetestplots=(s);利用lifetest过程进行生存分析并作生存函数图timet*c(0);指定时间变量和截尾变量并指出数据截尾时截尾变量的取值run;SAS软件输出结果Product-LimitSurvivalEstimatesSurvivalStandardNumberNumbertSurvivalFailureErrorFailedLeft0.0001.0000000196.0000.94740.05260.051211819.00

3、00.89470.10530.070421732.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.

4、000*...84310.000*...83316.000*...82335.000*...81346.000*...80NOTE:Themarkedsurvivaltimesarecensoredobservations.SummaryStatisticsforTimeVariabletQuartileEstimatesPoint95%ConfidenceIntervalPercentEstimate[LowerUpper)75...50.94.000.2542.00032.000.MeanStan

5、dardError181.70124.497NOTE:ThemeansurvivaltimeanditsstandarderrorwereunderestimatedbecausethelargestobservationwascensoredandtheestimationwasrestrictedtothelargesteventtimeSummaryoftheNumberofCensoredandUncensoredValuesPercentTotalFailedCensoredCensored

6、1981157.89SAS软件输出结果解释该结果包含四个部分：第一部分用乘积极限法估计了生存率（Survival），死亡率（Failure），生存率的标准误（SurvivalStandardError），死亡例数（NumberFailed）和该时间点前的生存例数（NumberLeft）。其中带有*号的表示截尾；第二部分给出了关于生存时间的描述性统计量，包括75%,50%和25%分位数以及相应的95%的可信区间（95%ConfidenceInterval），还有均数（Mean）和标准误（Standard

7、Error）从结果可以看出平均生存时间为181.701天；第三部分列出了完全数据（Failed），截尾数据（Censored）的例数，以及截尾数据占全部数据的百分比（PercentCensored）。最后是生存曲线图。教材中的说明现用Kaplan-Meier法计算生存率，步骤如下：(1)将所有生存时间按从小到大顺序排列（见表14-2第（2）列）并标上序号（第（1）列）。(2)列出各t时点前的存活病例数n（第（3）列）、各个时间点的死亡人数d（第（4）列）和截尾人数c（第（5）列）。(3)计算各t时刻的

8、死亡概率（第（6）列）。例如生存时间为32天时，死亡概率为。(4)计算各t时刻的生存概率(第（7）列)。例如生存时间为32天时，生存概率为。(5)计算各t时刻的生存率(第（8）列)。例如生存时间为32天时，生存率为，由此验证了在没有截尾数据的情况下，式(14-4)与式(14-5)是相同的。(6)以时间t为横指标，生存率为纵指标，作生存曲线图（图14－1）。表14-2Kaplan-Meier法计算生存率的计算用表　　　(1)(2)(3)(4)(5)(6)(