topic Re: proc lifetest add number of events at specific time points in Statistical Procedures

proc lifetest add number of events at specific time points

Li17 — Tue, 21 Feb 2017 14:36:27 GMT

Hello everyone,

I am new here and I am working on a project where I need to produce KM survival plots, which not only provide the number of patients at risk at specific time-points in a table below the graph but, in addition, I need to provide the number of events, which happend up to these timepoints.

I use the following code to produce the number of people at risk at the time points 0, 200, ... 1000 but I have no idea how I can add the number of events that happened at the same time-points. I know how I could add an additional table with the total number of events during follow-up, but I need to add the more detailed information for the specific time-points. Does anybody know how I can add this?

proc lifetest data=one plots=s(atrisk(outside(0.15))=0 to 1000 by 200); time totexp1 * outexp(0); run;

Thank you very much in advance.

Best regards,

Re: proc lifetest add number of events at specific time points

RW9 — Tue, 21 Feb 2017 14:48:49 GMT

You may be better off taking the output from lifetest and then running that (or doing it directly in) GTL. Then you can plot whatever data you want, for example these posts may help:

http://blogs.sas.com/content/graphicallyspeaking/2014/02/09/survival-plot/

http://blogs.sas.com/content/graphicallyspeaking/2011/11/28/survival-plot-with-external-at-risk-table/

Illustrates plotting lifetest outputs with other information below.

Re: proc lifetest add number of events at specific time points

Li17 — Mon, 27 Feb 2017 10:44:58 GMT

Thank you for the advice. I'll try that.

Re: proc lifetest add number of events at specific time points

XG — Thu, 23 Feb 2023 22:31:08 GMT

Hi Li, how are you. I ran into the same questions and wondering if you solve this problem and could you share your solution/code?

Thank you very much!

Re: proc lifetest add number of events at specific time points

Reeza — Thu, 23 Feb 2023 22:47:13 GMT

Something like the data in this chart?

https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/statug/statug_lifetest_examples02.htm

FYI - responding to an almost 6 year old thread isn't super useful. It's better to start your own thread.