l********l 发帖数: 10 | 1 我想fit一个Cox phreg.用的是如下的code:
proc tphreg data=fortest;
class trtgrp sex race;
model days*censor(0)=trtgrp sex race;
strata site;
output out=curve survival=est;
run;
这儿trtgrp=treatment group; days=survival time; censor=event/no event.
得到的survival probility 保存在curve,结果试图画图时发现curve dataset有以下问
题:
survival probility不是随survival time的增加而递减的阶梯函数。
不知道是什么地方出的问题,牛人们能给指点一下吗?
BTW,我用的是sas 9.1 | W**********E 发帖数: 242 | 2
我最近也碰到类似问题。我的想法:
1.
K-M CURVE 是UNADJUSTED. 但是你用TPHREG的话,你估计的SURVIVAL PROBABILITIES
ADJUSTED FOR COVARIATES(S(T)=S(T)0^EXP(BETA*X))。 SURVIVAL
PROBABILITIES 不仅仅随时间变化,也随这些COVARIATES数值的变化而变化.
打个比方,1月内生存概率S(T)应该比第2个月的高。但由于这些协变量X1,X2,
X3。那么不一定,有可能接受治疗+男性+白人第2月内比不治疗+女性+黑人第1月内
的生存概率反而高。
但你作图又无法考虑这些协变量的不同值的话,就出现上下波动,虽然趋势可能也是往下。
2.
用tphreg这个参数来输出SURVIVAL PROBABLITIES。你自己仔细看MANUAL
baseline out=;
【在 l********l 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我想fit一个Cox phreg.用的是如下的code: : proc tphreg data=fortest; : class trtgrp sex race; : model days*censor(0)=trtgrp sex race; : strata site; : output out=curve survival=est; : run; : 这儿trtgrp=treatment group; days=survival time; censor=event/no event. : 得到的survival probility 保存在curve,结果试图画图时发现curve dataset有以下问 : 题:
| l********l 发帖数: 10 | 3 谢谢解答!看了一些reference,感觉清楚了一些。 |
|