文档介绍：生存分析的主要内容
生存过程的描述
生存时间的分布特点
估计生存率
估计平均生存时间
绘制生存率曲线
生存过程的统计推断
生存时间的影响因素研究
生存分析和cox比例风险模型
2021/1/25
1
生存分析需具备的资料
生存时间或称时间变量
观察结果或称截尾变量
暴露变量及混杂因素
生存分析和cox比例风险模型
2021/1/25
2
肾移植例
一般手术(A组):
3 9 15 20 20 26 30 41 46 64+ 64 135
223 365 450 596+ 680+ 900+ 900+
改进手术(B组):
10 70+ 70+ 120 225 366 390+ 475+
518+ 647+ 801+ 1001+ 1045+ 1045+
生存分析和cox比例风险模型
2021/1/25
3
手术方式
生存时间
观察结果
treat
time
outcome
1
3
1
1
9
1
1
41
1
1
46
1
1
64
1
1
64
0
……
……
……
1
900
0
2
10
1
……
……
……
2
801
0
2
1001
0
2
1045
0
2
1045
0
生存分析和cox比例风险模型
2021/1/25
4
生存资料的定义(1)
将数据定义为生存资料格式
stset 时间变量 [，failure(截尾变量)]
failure (截尾变量) stata以所有不等于0的非缺失值为预期的出现结果；
若failure选择项缺省，则stata默认所有纪录均出现结果
生存分析和cox比例风险模型
2021/1/25
5
生存资料的定义(2)
四个变量---定义数据库后，stata系统自动产生的
_st /* 数据库中该条记录是否被定义为生存资料
_d /* 该条记录是否出现预期结果
_t /* 观察对象被随访的时间
_t0 /* 观察对象第一次被观察到的时间(开始过程的时间为0)
生存分析和cox比例风险模型
2021/1/25
6
stset time,failure(outcome)
failure event: outcome != 0 & outcome < .
obs. time interval: (0, time]
exit on or before: failure
------------------------------------------------------------------
33 total obs.
0 exclusions
------------------------------------------------------------------
33 obs. remaining, representing
18 failures in single record/single failure data
11370 total analysis time at risk, at risk from t = 0
earliest observed entry t = 0
last observed exit t = 1045
生存分析和cox比例风险模型
2021/1/25
7
生存资料的描述
stci ,[by(分组变量)选择项]
p(#) /* 生存时间的百分位数
level(#) /* 可信区间的可信度
media /* 计算中位生存时间，为默认选项
rmean /* 计算平均生存时间
生存分析和cox比例风险模型
2021/1/25
8
. stci,by(treat)
failure _d: outcome
analysis time _t: time
| no. of
treat | subjects 50% Std. Err. [95% Conf. Interval]
-------+----------------------------------------