获取SAS中的生存函数值

Question

我从指数分布生成了一个随机样本并对它们进行了排序，使它们从最低值到最高值，给出了我的订单统计信息。现在我需要在这些数字处获得生存函数的值，并根据订单统计数据绘制它们。我似乎无法弄清楚如何在SAS中获取这些生存值的列表，因此我可以绘制它们。

Answer 1

观察到的事件时间的幸存者函数估计可以从SAS的LIFETEST程序中获得。

proc lifetest data = yourdata outsurv = survest;
    time time2event*eventindicator(0);
run;

OUTSURV =选项创建的输出数据集包含幸存者函数估计值。请注意，默认情况下，这会给出幸存者函数的Kaplan-Meier估计值。

LIFETEST自动生成估计的幸存者功能随时间变化的图。如果要根据订单统计信息绘图，请在计算订单统计信息后使用SGPLOT过程中的STEP语句。

proc sgplot data = survandorder;
    step x = order y = survival;
run;