获得平滑的危害估计

Question

我如何获得数据集或 table 以及当我运行例如

时显示的平滑危害估计和随时间显示的 95% 置信区间

ods graphics on;
proc lifetest data=melanoma34 plots=(h(cl));
  time weeks*censor(1);
run;

编辑这就是我运行

data melanoma34;
  infile '/folders/myfolders/amsus/data/melanoma34.dat';
  input weeks status$;
  if status='alive' then censor=1;
     else censor=0;
run;

Proc lifetest data=Melanoma34 outsurv=hazards;
  time weeks*censor(1);
ods output hazardplot=smoothedhazard;
run;

并且 hazards outsurv 不包含平滑的危害估计

。数据可以在/amsus/data/

下找到here

Answer 1

outsurv 数据集包含危险函数估计值和置信区间。

Proc lifetest data=Melanoma outsurv=hazards;

编辑：以上生成了风险函数的估计值，而不是平滑估计值，并且仅适用于 method=LT.

要获得平滑估计，使用 method=KM（默认），然后使用以下行：

ods output hazardplot=smoothedhazard;

最终代码：

data melanoma34;
  infile '/folders/myfolders/amsus/data/melanoma34.dat';
  input weeks status$;
  if status='alive' then censor=1;
     else censor=0;
run;

Proc lifetest data=Melanoma34 plots=h(cl));
  time weeks*censor(1);
ods output hazardplot=smoothedhazard;
run;

获得平滑的危害估计

Obtain the smoothed hazard estimates

sas

survival-analysis