给出具有数天事件的数据,以及结果,如下:
data pretend ;
do subject = 1 to 1000 ;
fup_time = round(uniform(83386)*500, 1) ;
select(round(uniform(778523)*5, 1)) ;
when(1) outcome = 'cholys' ;
when(2) outcome = 'death' ;
when(3) outcome = 'tx end' ;
when(4) outcome = 'vascul' ;
otherwise outcome = 'reop' ;
end ;
output ;
end ;
label
fup_time = "The day on which ::outcome:: occurred"
outcome = "What the subject's last observed event was"
;
run ;
为每个结果生成曲线的最简单方法是什么,这些结果显示在第X天仍然可以观察到样本的比例,结果是什么?
我试过了:
proc lifetest data = pretend plots = all notable ;
strata outcome ;
time fup_time*censored(1) ;
run ;
在哪里审查'只要结果是'tx end'就设置为1或者“死亡”。我非常喜欢产生的产品极限生存曲线,除了死亡和生产线之外。 tx末端完全平坦,y = 1.0。
我实际上根本不想在这里做任何推断 - 只想要漂亮的照片。有一个简单的方法吗?