LIFETEST Procedure（1）
(2014-06-09 23:12:45)
转载▼
标签： sas lifetest | 分类： SAS笔记 |
PROC LIFETEST ;
BY variables ;
FREQ variable ;
ID variables ;
STRATA variable <(list)> <...variable <(list)>> ;
TEST variables ;
TIME variable <*censor(list)> ;
proc lifetest选项:
Input and Output Data Sets | |
DATA= | 确定输入的数据集； |
OUTSURV= | 确定输出数据集，其内包括生存函数估计和置信区间； |
OUTTEST= | 确定输出数据集，其内包括秩检验统计； |
Nonparametric Estimation | |
INTERVALS= | 用寿命表估计时用来确定区间端点； 例如：intervals=5,10 to 30 by 10； 则表示区间{[0,5),[5,10),[10,20),[20,30),[30,∞)} |
NELSON | 表示加上 Nelson-Aalen 估计，正常默认只有乘积项估计； |
METHOD= | 确定计算生存估计的方法，总共两种，一种是PL（乘积项估计方法），一种是LT（寿命表方法）； |
NINTERVAL= | 确定采用寿命表方法来估计的区间个数； |
WIDTH= | 确定区间宽度； |
Confidence Limits for Survivorship | |
ALPHA= | 区间估计的置信水平； |
BANDMAXTIME= | 置信区间最大时间点； |
BANDMINTIME= | 置信区间最小时间点； |
CONFBAND= | 确定在OUTSURV= 的数据集里置信区间的类型，分为两种：EP和HW, confband=all（表示同时包含两种）/EP/HW； |
CONFTYPE= | 表示利用生存函数来求置信范围的方法，有下面几种： ASINSQRT 反正弦平方根转化 LOGLOG 对数的对数转换 LINEAR 线性转换 LOG 对数转换 LOGIT 多元逻辑转换 ； |
Line Printer Plots（表示打印输出plot） | |
FORMCHAR(1,2,7,9)= | 确定打印的plot坐标轴设置； |
LINEPRINTER | 明确打印图形； |
MAXTIME= | 确定plot的最大时间点； |
NOCENSPLOT | 禁止plot上的缺失值的出现； |
PLOTS= | 确定需要呈现的plot：C、S、LS、LLS、H、P； |
ODS Graphics（利用ODS来输出plot） | |
MAXTIME= | 确定plot的最大时间点； |
PLOTS= | 确定需要呈现的plot：C、S、LS、LLS、H、P；在这里可以确定核平滑估计的方法：U（均匀核）、B（二重权核）、E（Epanechnikov核）;同时可以得到确定的点的函数估计值。 例如：plots=(survival(atrisk=100 to 350 by 50 cb=all/ep/hw strata=unpack/overlay/panel test) logsurv)/这个atrisk就是需要确定的点，cb=表示在图上呈现的生存函数的置信区间，默认是hw strata=表示图形程是否呈现在同一个图上，默认为overlay，test为同质性检验/ plots=hazard(kernel=B)/表示危险函数核平滑估计方法为二重权核/ |
Traditional Graphics | |
ANNOTATE= | 确定图形模板数据集； |
CENSOREDSYMBOL= | 确定缺失数据在图中的符号； |
DESCRIPTION= | 对图形进行描述； |
EVENTSYMBOL= | 确定事件观察值的符号； |
GOUT= | 确定保存输出图形的目录位置； |
LANNOTATE= | 确定图形模板数据集，与group配合使用，每一个不同group，其模板都不一样，如果没有group，则和annotate=数据集有着相同的效果，不能和LINEPRINTER语句以及ODS GRAPHIC语句共存； |
MAXTIME= | 确定图形最大时间点； |
PLOTS= | 确定需要呈现的plot：C、S、LS、LLS、H、P； |
Control Output（控制生成数据的输出） | |
NOPRINT | 禁止在output窗口呈现； |
NOTABLE | 禁止输出生存函数估计； |
INTERVALS= | 表示只呈现每一个确定区间的最小时间点的乘积项估计； |
TIMELIST= | 确定乘积项估计呈现的时间点； |
| variable什么意思中文 REDUCEOUT | 确定 OUTSURV=只只包含INTERVAL= or TIMELIST=的观测值； |
Miscellaneous（其它参数） | |
ALPHAQT= | 确定生存时间四分位点的置信水平 |
MISSING | 表示允许缺失值作为一个stratum水平 |
SINGULAR= | sets the tolerance for testing singularity of covariance matrix of rank statistics确定秩检验统计的协方差矩阵的奇异点的可包容值，也就是； |
STDERR | 在OUTSURV=数据集中输出生存函数估计的标准误差； |
TIMELIM= | 利用此时间点来估计平均生存时间和其标准误差，有三种：event、observed、number（具体的数字，要求大于等于event时间点）； |
LIFETEST Procedure（2）
(2014-06-11 19:09:54)
转载▼
标签： sas lifetest | 分类： SAS笔记 |
TIME STATEMENT：
time variable*censor(list);
variable是失败时间变量名；
例如：time T*flag(1,2);
T是观测的事件发生时间，flag代表删失变量，（1,2）分别表示删失，对应的T表示删失值；
STRATA Statement |
STRATA variable <(list)> <...variable < (list)> > ;
例如：strata age(5,20 to 50 by 10) sex;
表示age里面有6个levels：
{(-∞,5),[5,20),[20,30),[30,40),[40,50),[50,∞)};
options:
group= 表示对分层检验进行分组；
nodetail= 禁止输出统计检验和协方差矩阵
notest 禁止任何检验
test= 确定检验权函数，总共有这么几种：Fleming（p1， p2）/logrank/none/lr/modpeto/peto/wilcoxon/tarone，默认为test=（logrank wilcoxon lr）
missing 表示缺失值为一个水平
多重比较：
adjust= 确定一个多重比较的检验方法：bon/dunnett/scheffe/sidak/get/tukey/simulate
diff= 确定多重比较的一个基组
例如：DIFF=CONTROL('large' 'F');表示（large，F）这一组与其他组进行多重统计检验
发表评论