turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- PHREG

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

03-10-2014 01:18 PM

in PROC PHREG, is it possible to obtain a "failure" plot instead of a survival plot?

for example, below will produce an adjusted survival curve.

proc phreg data=x plots(overlay)=s;

class trt;

model;

baseline covariates=x out=y SURVIVAL=Survival /diradj group=trt;

run;

what I'm interested in is a failure curve to show cumulative probability of failure

In PROC LIFETEST when can request to change the default survival to failure curve. Is this available in PHREG?

proc lifetest data=x plot=(survival ( failure));

..

run;

Thanks

Accepted Solutions

Solution

03-10-2014
01:23 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

03-10-2014 01:23 PM

Hi,

No, the probability plot of failure is not possible to obtain directly from phreg - not as far as I know.

But, it is not more difficult to obtain than to calculate 1-survivalprobability.

data y;

set y;

prob_of_failure=1-survival;

run;

Jacob

All Replies

Solution

03-10-2014
01:23 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

03-10-2014 01:23 PM

Hi,

No, the probability plot of failure is not possible to obtain directly from phreg - not as far as I know.

But, it is not more difficult to obtain than to calculate 1-survivalprobability.

data y;

set y;

prob_of_failure=1-survival;

run;

Jacob

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

03-11-2014 08:18 AM

Thanks for response. I couldn't get it to work using any of the proc options. However, I was able to manipulate the template for the survival plot (Stat.Phreg.Graphics.Survival), and specifying (y=eval(1-SURVIVAL) ) in the seriesplot statement.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-25-2016 09:33 AM

That is a good solution for those comfortable with modifying the underlying ODS Graphics template for the SurvivalPlot.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-25-2016 09:32 AM

JacobSimonsen is correct. You can post-process the BASELINE data set to compute failure = 1 - survival and then plot failure versus time with PROC SGPLOT (for an ODS Graphics type plot) or PROC GPLOT (for a traditional graphics type plot).