加载中…LIFETEST Procedure（1）
(2014-06-09 23:12:45)
标签：
saslifetest |
分类： 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 |
|
|---|---|
|
确定输入的数据集； |
|
|
确定输出数据集，其内包括生存函数估计和置信区间； |
|
|
确定输出数据集，其内包括秩检验统计； |
|
|
Nonparametric Estimation |
|
|
用寿命表估计时用来确定区间端点； 例如：intervals=5,10 to 30 by 10； 则表示区间{[0,5),[5,10),[10,20),[20,30),[30,∞)} |
|
|
表示加上 |
|
|
确定计算生存估计的方法，总共两种，一种是PL（乘积项估计方法），一种是LT（寿命表方法）； |
|
|
确定采用寿命表方法来估计的区间个数； |
|
|
确定区间宽度； |
|
|
Confidence Limits for Survivorship |
|
|
区间估计的置信水平； |
|
|
置信区间最大时间点； |
|
|
置信区间最小时间点； |
|
|
确定在OUTSURV= 的数据集里置信区间的类型，分为两种：EP和HW, confband=all（表示同时包含两种）/EP/HW； |
|
|
表示利用生存函数来求置信范围的方法，有下面几种：
ASINSQRT
LINEAR
LOG
LOGIT |
|
|
Line Printer Plots（表示打印输出plot） |
|
|
确定打印的plot坐标轴设置； |
|
|
明确打印图形； |
|
|
确定plot的最大时间点； |
|
|
禁止plot上的缺失值的出现； |
|
|
确定需要呈现的plot：C、S、LS、LLS、H、P； |
|
|
ODS Graphics（利用ODS来输出plot） |
|
|
确定plot的最大时间点； |
|
|
确定需要呈现的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 |
|
|
确定图形模板数据集； |
|
|
确定缺失数据在图中的符号； |
|
|
对图形进行描述； |
|
|
确定事件观察值的符号； |
|
|
确定保存输出图形的目录位置； |
|
|
确定图形模板数据集，与group配合使用，每一个不同group，其模板都不一样，如果没有group，则和annotate=数据集有着相同的效果，不能和LINEPRINTER语句以及ODS
GRAPHIC语句共存； |
|
|
确定图形最大时间点； |
|
|
确定需要呈现的plot：C、S、LS、LLS、H、P； |
|
|
Control Output（控制生成数据的输出） |
|
|
禁止在output窗口呈现； |
|
|
禁止输出生存函数估计； |
|
|
表示只呈现每一个确定区间的最小时间点的乘积项估计； |
|
|
确定乘积项估计呈现的时间点； |
|
|
确定 |
|
|
Miscellaneous（其它参数） |
|
|
确定生存时间四分位点的置信水平 |
|
|
表示允许缺失值作为一个stratum水平 |
|
|
sets the tolerance for testing singularity of covariance matrix of rank statistics确定秩检验统计的协方差矩阵的奇异点的可包容值，也就是； |
|
|
在OUTSURV=数据集中输出生存函数估计的标准误差； |
|
|
利用此时间点来估计平均生存时间和其标准误差，有三种：event、observed、number（具体的数字，要求大于等于event时间点）； |
|
喜欢
0
赠金笔