Proc Lifetest
PROC LIFETEST 的最简单用法是请求生存时间样本的幸存者函数的非参数估计值。在这种情况下,只需要 PROC LIFETEST 语句和 TIME 语句。
一、概述
二、语法
1. LIFETEST 过程:
PROC LIFETEST <options>
BY variables;/***by约等于strata? 是的!by比strata更高效,strata比by输出的多,机器够用,无脑strata***/
FREQ variable </ option>; /***该变量包含每个观测值的发生频率*** 默认variable为整数,如果missing,则不适用该观测值。且不为整数时会被截断,当在“/”后指定NOTRUNCATE或NOTRUNC时,频率值不截断。NOTRUNC不适用method=FH***/
ID variables;(ID 语句ID语句提供了一个变量列表,其值确定product-limit、Breslow 或 Fleming-Harrington 估计值中的观测值)
STRATA variable <(list)> <…variable <(list)>> </ options>;(分层,定义指定一个变量或一组变量,用于分析的分组;如果在 STRATA 语句中还指定了 GROUP= 选项,则会执行分层检验以检验由 GROUP= 变量定义的 k 个样本,同时控制 STRATA 变量的影响。举例:strata Sex / group=Treatment(评估治疗效果,同时调整性别差异);)/***strata Age(5,20 to 50 by 10) Sex;这是6组**/
/***option list***/
/***Homogeneity Tests***/
/***GROUP= 指定分层检验的组变量***/
/***NODETAIL 禁止打印检验统计和协方差矩阵***/
/***NOTEST 禁止任何检验***/
/***TEST= 指定对应于各种权重函数的检验***/
/***TREND 要求进行趋势测试***/
/***Multiple Comparisons***/
/***ADJ= & ADJUST= 请求多次比较调整 ***/
/***DIFF= 指定要考虑的差异类型***/
/***Missing Strata Value***/
/***MISSING 允许缺失值作为有效的stratum值***/
/***Display Option***/
/***NOLABEL 在显示中使用STRATA 变量的名称***/
TEST variables;(指定数字协变量与survival time的关联列表
TIME variable <*censor(list)> </ option>;(必选,定义survival time和censoring indicator的变量,censoring 0=终点事件,1=censoring。KM法censor(1),逆KM法(0)即把删失和终点事件对调)
WEIGHT variable;
2.output(SAShelp download,自行阅读)
图1 该图列出了每个观测时间的生存估计值、失败率、估计值的标准误差、累计失败数以及研究中剩余的受试者数
图 2 显示了该组生存时间的汇总统计信息。它包含第 25、50 和 75 百分位的估计值以及相应的 95% 置信限。平均生存时间及其标准误差被低估了,因为最大的观测值被删失,估计值被限制在最大的事件时间。
图5:截尾和事件观察的数量汇总。该图列出了每一层的事件和删失的观测数,以及删失的观测数的百分比。
图6 乘积极限生存函数估计值与生存时间的关系图。两种治疗方法的主要区别在于存活时间更长。使用ATRISK=option 显示特定时间点处于风险中的受试者数量。
3. ODS output=datasets
- Product-Limit Survival Estimates :The ODS name of this table is ProductLimitEstimates.
- Breslow Survival Estimates:The ODS name of this table is BreslowEstimates.
- Fleming-Harrington Survival Estimates:The ODS name of this table is FlemingEstimates.
- Quartile Estimates:The ODS name of this table is Quartiles.
- Mean Estimate:The ODS name of this table is Means.
- Life-Table Survival Estimates:The ODS name of this table is LifetableEstimates.
- Summary of the Number of Censored and Uncensored Values:The ODS name of this table is CensoredSummary.
- Rank Statistics:The ODS name of this table is HomStats.