Overlaying life test plots with phreg

[sjmosiman]

Hi all,

I'm trying to stratify my survival curves however my output gives me three separate plots rather than overlaying them.  I've seen the use of baseline/rowid but I can't seem to make that work. Also I cant use lifetest procedure because I'm using a continuous variable anyone ahs suggestions they would be much appreciated, my current code is below:

ods graphics on;
proc phreg data=merged11 plot(overlay=stratum)=(survival);
where rs>34;
class chemoYN (param=ref ref='N');
model dx_lastcontact_death_months * Puf_Vital_status(1) = chemoYN |absBenefitTAM1;
hazardratio chemoyn / at(absBenefitTAM1= (5 10 15 20 25 30)) units=5;
strata T1b1c;
run;

[OsoGris · Accepted solution]

Since you have no BASELINE statement or COVARIATES= data set you have only one set of covariates being plotted by default. This is the reference level of the CLASS variables and the means of any continuous covariates. As such, you might want to try the OVERLAY=BYROW option which "displays, for each covariate set, a separate plot containing the curves for all the strata." This is according to the PHREG documentation on the OVERLAY option.