-
LIFETEST Procedure
(
1<
/p>
)
(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=
确定输入的数据集;
p>
OUTSURV=
OUTTEST=
确定输
出数据集,其内包括生存函数估计和置信区间;
确定输出数据集,其内包括秩检验统计;
Nonparametric Estimation
用寿命表估计时用来确定区间端点;
INTERVALS=
例如:
intervals=5,10
to 30 by 10
;
则表示区
间
{[0,5),[5,10),[10,20),[20,30),[30,∞)}<
/p>
NELSON
表示加上
Nelson-Aalen
估计,正常默认只有乘积项估计;
确定计算生存估计的方法,总共两种,一种是
PL
(乘积项估计方法),一种是
METHOD=
LT<
/p>
(寿命表方法);
NINTERVAL
=
WIDTH=
确定采用寿命表方法来估计的区间个数;
确定区间宽度;
Confidence Limits for Survivorship
ALPHA=
区间估计的置信水平;
BANDMAXTIME=
BANDMINTIME=
CONFBAND=
置信区间最大时间点;
置信区间最小时间点;
确定在
OUTSURV=
的数据集里
置信区间的类型,分为两种:
EP
和
H
W,
confband=all
(表示同时包含两种)
/EP/HW
;
表示利用生存函数来求置信范围的方法,有下面几种:
ASINSQRT
反正弦平方根转化
LOGLOG
对数的对数转换
CONFTYPE=
LINEAR
线性转换
LOG
对数转换
LOGIT
多元逻辑转换
Line
Printer Plots
(表示打印输出
plot
)
FORMCHAR(1,2,7,9)=<
/p>
确定打印的
plot
坐标轴设置;
;
LINEPR
INTER
MAXTIME=
明确打印图形;
< br>
确定
plot
的最大时间点;
NOCENSPLOT
PLOTS=
禁止
plot
上的缺失值的出现;
p>
确定需要呈现的
plot
:
C
、
S
、
LS
、
LLS
< br>、
H
、
P
;
ODS Graphics
(利
用
ODS
来输出
plot
)
MAXTIME=
确定
plot
的最大时间点;
确定需要呈现的
plot
:
C
、
S
、
L
S
、
LLS
、
H
、
P
;在这里可以确定核平滑估计的
方法:
U
(均匀核)、
B
(二重权核)、
E
(
Epanechnikov
核)
;
同时可以得到
确定的点的函数估计值。
PLOTS=
例如:
plots=(survi
val(atrisk=100 to 350 by 50 cb=all/ep/hw
strata=unpack/overlay/panel test) logsu
rv)/
这个
atrisk
就是需要确
定的
点,
cb=
表示在图上呈现的生存
函数的置信区间,默认是
hw strata=
表示图形
程是否呈现在同一个图上,默认为
overlay
,
test
为同质性检验
/
p>
plots=hazard(kernel=B)/
表示危险函数核
平滑估计方法为二重权核
/
Traditional
Graphics
ANNOTATE=
确定图形模板数据集;
CENSOREDSYMBOL=
确
定缺失数据在图中的符号;
DESCRIPTION=
对图形进行描述;
EVENTSYMBOL
=
确定事件观察值的符号;