I extracted this code from one of my 2018 PharmaSUG examples. My actual example is more complicated. Come to Seattle to see the entire presentation. My abstract is below the code. The code might have a few superfluous lines left, but it illustrates the principles involved. As you note, there is no way to have an axis table and a 0,0 origin, but there is a way to create the at risk table without using an axis table. It relies on the graph modification method that I call "highly customized graphs". I have written about it at SAS Global Forum, PharmaSUG, in the Graphically Speaking blog, https://blogs.sas.com/content/graphicallyspeaking/ , and elsewhere including here: https://support.sas.com/documentation/prod-p/grstat/9.4/en/PDF/odsadvg.pdf
data simdata(drop=n);
Treatment = 'Drug ';
do until(n = 500);
Month = -100 * log(uniform(368));
status = not (uniform(368) > 0.10 or (Month > 15 and uniform(368) > 0.2));
if Month < 18 then do; output; n + 1; end;
end;
n = 0;
Treatment = 'Placebo';
do until(n = 500);
Month = -150 * log(uniform(368));
status = not (uniform(368) > 0.16 or (Month > 15 and uniform(368) > 0.1));
if Month < 18 then do; output; n + 1; end;
end;
run;
proc template;
delete Stat.Lifetest.Graphics.ProductLimitFailure / store=sasuser.templat;
source Stat.Lifetest.Graphics.ProductLimitFailure / file='tpl.tpl';
quit;
* Modify template for full-size graph;
data _null_;
infile 'tpl.tpl' end=eof;
input;
* Add PROC TEMPLATE statement;
if _n_ = 1 then call execute('options nosource; proc template;');
* Remove offset, make the Y axis start precisely at zero;
_infile_ = tranwrd(_infile_, 'yaxisopts=(', 'yaxisopts=(offsetmin=0 offsetmax=0 ');
* Remove offset, make the X axis start precisely at zero;
_infile_ = tranwrd(_infile_, 'offsetmin=.05', 'offsetmin=0 offsetmax=0');
* Remove the at-risk table, change Boolean expression to unconditional false;
_infile_ = tranwrd(_infile_, 'PLOTATRISK=1', '0');
* Find beginning of LAYOUT OVERLAY statement;
if find(_infile_, 'layout overlay / xaxisopts') then lo + 1;
* Find end of LAYOUT OVERLAY statement;
if lo and index(_infile_, ';') then do;
* Reset layout flag;
lo = 0;
* Add padding;
_infile_ = tranwrd(_infile_, ';', 'pad=(bottom=10% left=5%);');
* Write original or modified template line;
call execute(_infile_);
%* Add ANNOTATE statement to LAYOUT OVERLAY and give it an ID;
_infile_ = 'annotate / id="id";';
end;
* Write original, new, or modified template line;
call execute(_infile_);
* End PROC TEMPLATE, display source again;
if eof then call execute('quit; options source;');
run;
* Create data object for graph (failureplot=fp);
ods document name=MyDoc (write);
ods graphics on;
proc lifetest data=simdata notable
plot=survival(failure nocensor atrisk(outside atrisktickonly)=0 to 18 by 3);
time Month * status(0);
strata treatment;
ods output failureplot=fp homtests=p;
run;
ods document close;
* List document contents, need to copy the path from the output;
proc document name=MyDoc;
list / levels=all;
quit;
* Using failure plot path from the ODS document, output dynamic variables;
proc document name=MyDoc;
ods exclude dynamics;
ods output dynamics=dynamics;
obdynam \Lifetest#1\FailurePlot#1;
quit;
* SG annotation data set;
data anno(drop=atrisk tatrisk stratum stratumnum left probchisq);
retain Function 'Text' /* All annotations add text */
X1Space 'DataValue ' /* Most X coordinates are data values */
Y1Space 'GraphPercent' /* Most Y coordinates are graph percentages */
ID 'id' /* Constant ID. Matches ANNOTATE stmt ID. */
Width 100 /* Width=100% ensures no splitting */
left -3; /* Extrapolated data value for label column */
/* of axis table */
* Extract axis table from graph data object;
set fp(drop=_: censored event survival time);
by stratumnum;
if _n_ eq 1 then do;
Label = 'No. at Risk'; /* Header for axis table row labels */
x1 = left; /* Ad hoc X coordinate */
y1 = 9; /* Ad hoc Y coordinate */
Anchor = 'Left '; /* Left-justify header */
TextWeight = 'Bold'; /* Bold text */
TextStyleElement = 'GraphDefault'; /* Text style element */
output; /* Write header */
textweight = ' '; /* Back to normal font */
TextSize = .; /* Restore default text size */
end;
/* Set colors for axis table rows */
textstyleelement = cats('GraphData', stratumnum);
textweight = ' '; /* Normal font (not bold) */
if first.stratumnum then do; /* Set up axis table row labels */
label = stratum; /* Row label */
x1 = left; /* Ad hoc X coordinate */
y1 = (3 - stratumnum) * 3; /* Ad hoc Y coordinate */
anchor = 'Left'; /* Left-justify header */
output; /* Write row label */
anchor = 'Right'; /* Right-justify */
y1space = 'DataValue'; /* Use data value Y coordinate */
x1 = 97.5; /* Ad hoc Y coordinate */
y1 = ifn(stratumnum eq 1, 75, 90); /* Ad hoc Y coordinates */
TextSize = 7; /* Text size for ad hoc curve labels */
x1space = 'GraphPercent'; /* Graph percentage X coordinates */
output; /* Write out ad hoc curve labels */
textsize = .; /* Restore default text size */
x1space = 'DataValue'; /* Restore X coordinate space */
end;
if n(tatrisk); /* If this is part of the at-risk table */
label = put(atrisk, 6. -L); /* Format numeric value to string */
x1 = tatrisk; /* X coordinate from data set */
y1 = (3 - stratumnum) * 3; /* Ad hoc Y coordinate */
y1space = 'GraphPercent'; /* Y coordinate is graph percentage */
anchor = 'Center'; /* Center values relative to ticks */
output; /* Write out body of axis table */
run;
* Display SG annotation data set;
proc print noobs; run;
ods html body='b.html';
* Write SGRENDER code, populate dynamics;
data _null_;
* Read dynamics;
set dynamics(where=(label1 ne '___NOBS___')) end=eof;
* Write PROC statement, start DYNAMIC statement;
if _n_ = 1 then do;
call execute('proc sgrender data=fp sganno=anno ' ||
'template=Stat.Lifetest.Graphics.ProductLimitFailure;');
call execute('dynamic');
end;
* Write name/value pairs. Numeric: name=formatted-value.
Character: name=quoted formatted-value;
if cvalue1 ne ' ' then
call execute(catx(' ', label1, '=',
ifc(n(nvalue1), cvalue1,
quote(trim(cvalue1)))));
* End the DYNAMIC statement, end the PROC SGRENDER call;
if eof then call execute('; run;');
run;
ods html close;
proc template;
delete Stat.Lifetest.Graphics.ProductLimitFailure / store=sasuser.templat;
quit;
Creating a Publication-Quality Graph Embedded in Another Graph
Warren F. Kuhfeld, SAS Institute Inc., Cary, NC
John King, Ouachita Clinical Data Services, Inc.
Several recent papers in the New England Journal of Medicine (NEJM) display graphs inside of graphs (Marso et al.
2016b, a; Zinman et al. 2015; Holman et al. 2017). Typically, both graphs are Kaplan-Meier plots. The larger graph
has a Y axis that ranges from 0 to 100. The failure curves occupy a relatively small part of the graph and extend
perhaps no more than 10%–20% of the way up the Y axis. Inside this graph, in the open white space, you can place
another graph in which the Y axis extends from 0 to slightly above the maximum value. The outer graph shows the
failure probability in the full 0%–100% context, whereas the inner graph zooms in to the results and provides additional
text and annotation, such as curve labels, the hazard ratio, and p-values. The outer graph also displays the number of
subjects at risk. The format for the NEJM is quite specific. This paper illustrates how to use ODS Graphics to create a
graph inside a graph that is suitable for publication in the NEJM.
