获得平滑的危害估计值

Question

如何获得具有平滑的危险估计值的数据集或表格以及我运行时显示的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

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

Proc lifetest data=Melanoma outsurv=hazards;

编辑： 以上内容生成危险函数的估计值，而不是平滑估计值，仅适用于方法= 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;