继续浏览精彩内容
慕课网APP
程序员的梦工厂
打开
继续
感谢您的支持,我会继续努力的
赞赏金额会直接到老师账户
将二维码发送给自己后长按识别
微信支付
支付宝支付

Proc Phreg

乌然娅措
关注TA
已关注
手记 64
粉丝 21
获赞 12

SAS day 17: Proc Phreg

One day, my boss took a glance at a table with Hazard Ratio and Median Survival Time then he told me the program set the reference group in Proc Phreg wrong.

It turns out he was correct after validating the program. However, I was very curious about how did he figure it out by an Augenblick. Then he shared with us some knowledge about Hazard Ratio and Proc Phreg.

Background knowledge

Hazard Ratio: HRs represent instantaneous risk over the study time period.

Median Survival Time: The time after which 50 percent of people with a particular condition are still living, and 50 percent have died.

Median survival time and Hazard Ratio

if Median(A) > Median(B) then HR(A) < HR(B), vice versa.

Without Loss of Generality, if the median survival time of group A is significantly higher than group B, then Hazard Ratio of group A is smaller than Hazard Ratio of group B, vice versa.

[caption id=“attachment_1056” align=“alignnone” width=“750”]image

qimono / Pixabay[/caption]

**Example  **

Now we will demonstrate Proc Phreg with Hazard Ratio with Sashelp.BMT dataset.

Basic Proc Phreg Syntax

Proc phreg data=a;
class x ; /*ref(x) only character value, the lowest character value is reference*/
model x*y= a/rl ; /*can input more than one variable on the RHS of the equation. rl:risk limit (eßt)*/
by z;/*optional*/ 
id m;/*optional*/
run;
  1. Read in the sashelp.BMT dataset and only consider AML Low Group and AML High Group
    Note: Hazard Ratio is usually designed for two-level comparison,
    we will introduce 3-level comparison next time.
data aml;
set sashelp.bmt;
if group="AML-Low Risk" then ref="0";
if group="AML-High Risk" then ref="1";
if ref in("0","1");
run;

image

2.  Apply Proc Phreg to generate the hazard ratio

proc phreg data=aml outest=est61 covout;
class ref;
model T*status(0)=ref/rl;
run;

image

Results:

From the output, we have Hazard Ratio=0.39, which means the probability a patient in AML high Group survives until today and will survive until tomorrow, compared with the AML Low Group Reduced by 0.39.

Validation:

As we mentioned earlier, if Hazard Ratio <1 then we expected a longer survival median time in AML Low Group. The following Kaplan Meier plot demonstrated the AML Low Group indeed has a larger survival median.

image

proc lifetest data=aml plots=survival;
time T*status(0);
strata Group;
run;

Thanks very much to Joe wang to share the Proc Phreg knowledge with me, we will show the 3-level Hazard Ratio in Proc Phreg next time!

Happy studying! 🐰

打开App,阅读手记
0人推荐
发表评论
随时随地看视频慕课网APP