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）０＾EXP（BETA＊X））。 SURVIVAL
PROBABILITIES 不仅仅随时间变化，也随这些COVARIATES数值的变化而变化.
打个比方，１月内生存概率S（T）应该比第２个月的高。但由于这些协变量X１，X２，
X３。那么不一定，有可能接受治疗+男性+白人第２月内比不治疗+女性+黑人第１月内
的生存概率反而高。
但你作图又无法考虑这些协变量的不同值的话，就出现上下波动，虽然趋势可能也是往下。
2.
用tphreg这个参数来输出SURVIVAL PROBABLITIES。你自己仔细看MANUAL
baseline out=;
【在 l********l 的大作中提到】
: 我想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,感觉清楚了一些。 |
|
