生存分析（三十）

 

前边小胖主要介绍了Cox PH model以及Stratified Cox model的一些基本内容，下边我们简单来看一下其SAS程序。Cox mode是通过SAS程序中的PROC PHREG来实现的。

 

（1）Cox PH model

 

PROC  PHREG;

  MODEL  TIME*STATUS(0)=TREAT X1 X2 / TIES=EXACT  RL;

RUN;

 

解释：

 

（1）TIME 代表从随机化到发生event或censored的时间

（2）STATUS中1代表发生event，0代表censored

（3）TREAT代表治疗组，X1,X2代表其他预后变量或者说协变量

（4）TIES=EXACT选项在很多event同时发生的情况下比较适用，默认的为TEST=BRESLOW

（5）RL选项可以使SAS输出hazard ratio的95％可信区间

 

 

（2）Stratified Cox model

 

PROC  PHREG;

  MODEL  TIME*STATUS(0)=TREAT X1 X2 / TIES=EXACT  RL;

  STRATA CENTER;

RUN;

 

解释：

 

假设CENTER不满足PH assumption，TREAT， X1， X2满足PH assumption，我们则可以采用上述的Stratified Cox model。与Cox PH model不同的是：

 

（1）MODEL中只包含满足PH assumption的变量

（2）多了一个STRATA语句，用来对不满足PH assumption的变量进行分层