Survival Plot with dashed line.

cypher85 — Wed, 20 Jun 2012 18:20:38 GMT

Hello,

I am trying to make a survival plot with a dashed line by modifying the Stat.Lifetest.Graphics.ProductLimitSurvival template. I cannot seem to figure out how to do this.

Any assistance would be greatly appreciated!

Here is some data and code as well as the template I am working with.

data ds;

input time event;

cards;

0.5 0

1 1

2 1

2 0

3 1

4 1

5 0

6 1

7 0

8 1

9 1

10 0

12 1

14 0

14 1

17 0

20 1

21 0

run;

proc lifetest data=ds plots=survival(cl);

time time*event(0);

run;

/******************************************************/

proc template;

define statgraph Stat.Lifetest.Graphics.ProductLimitSurvival / store =

SASUSER.TEMPLAT;

dynamic NStrata xName plotAtRisk plotCensored plotCL plotHW plotEP

labelCL labelHW labelEP maxTime StratumID classAtRisk plotBand

plotTest GroupName yMin Transparency SecondTitle TestName pValue;

BeginGraph;

if (NSTRATA=1)

if (EXISTS(STRATUMID))

entrytitle " " " for " STRATUMID;

else

entrytitle "Testing"; /*modifies Title*/

endif;

if (PLOTATRISK)

entrytitle "with Number of Subjects at Risk" / textattrs=

GRAPHVALUETEXT;

endif;

layout overlay / xaxisopts=(label = "x axis label" /*modifies x axis label*/ shortlabel=XNAME offsetmin=.05 linearopts

=(viewmax=MAXTIME /*modifies x max*/)) yaxisopts=(label="y axis label" /*modifies y axis label*/

shortlabel="Survival" linearopts=(viewmin=0 viewmax=1 /*modifies y axis min/max*/

tickvaluelist=(0 .2 .4 .6 .8 1.0 /*modifies y axis tick marks*/)));

if (PLOTHW=1 AND PLOTEP=0)

bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / modelname

="Survival" fillattrs=GRAPHCONFIDENCE name="HW" legendlabel=

LABELHW;

endif;

if (PLOTHW=0 AND PLOTEP=1)

bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / modelname

="Survival" fillattrs=GRAPHCONFIDENCE name="EP" legendlabel=

LABELEP;

endif;

if (PLOTHW=1 AND PLOTEP=1)

bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / modelname

="Survival" fillattrs=GRAPHDATA1 datatransparency=.55 name="HW"

legendlabel=LABELHW;

bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / modelname=

"Survival" fillattrs=GRAPHDATA2 datatransparency=.55 name="EP"

legendlabel=LABELEP;

endif;

if (PLOTCL=1)

if (PLOTHW=1 OR PLOTEP=1)

bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME /

modelname="Survival" display=(outline) outlineattrs=

GRAPHPREDICTIONLIMITS name="CL" legendlabel=LABELCL;

else

bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME /

modelname="Survival" fillattrs=GRAPHCONFIDENCE name="CL"

legendlabel=LABELCL;

endif;

stepplot y=SURVIVAL x=TIME / name="Survival" rolename=(_tip1=

ATRISK _tip2=EVENT) tip=(y x Time _tip1 _tip2) legendlabel=

"Survival";

if (PLOTCENSORED=1)

scatterplot y=CENSORED x=TIME / markerattrs=(symbol=plus /*modifies marker symbol*/) name=

"Censored" legendlabel="Censored";

endif;

if (PLOTCL=1 OR PLOTHW=1 OR PLOTEP=1)

discretelegend "Censored" "CL" "HW" "EP" / location=outside

halign=center;

else

if (PLOTCENSORED=1)

discretelegend "Censored" / location=inside autoalign=(

bottomleft);

endif;

if (PLOTATRISK=1)

innermargin / align=bottom;

blockplot x=TATRISK block=ATRISK / repeatedvalues=true display=

(values) valuehalign=start valuefitpolicy=truncate

labelposition=left labelattrs=GRAPHVALUETEXT valueattrs=

GRAPHDATATEXT (size=7pt) includemissingclass=false;

endinnermargin;

endif;

endlayout;

else

entrytitle " ";

if (EXISTS(SECONDTITLE))

entrytitle SECONDTITLE / textattrs=GRAPHVALUETEXT;

endif;

layout overlay / xaxisopts=(shortlabel=XNAME offsetmin=.05 linearopts

=(viewmax=MAXTIME)) yaxisopts=(label="Recurrence-free survival probability"

shortlabel="Survival" linearopts=(viewmin=0 viewmax=1

tickvaluelist=(0 .2 .4 .6 .8 1.0)));

if (PLOTHW)

bandplot LimitUpper=HW_UCL LimitLower=HW_LCL x=TIME / group=

STRATUM index=STRATUMNUM modelname="Survival" datatransparency=

Transparency;

endif;

if (PLOTEP)

bandplot LimitUpper=EP_UCL LimitLower=EP_LCL x=TIME / group=

STRATUM index=STRATUMNUM modelname="Survival" datatransparency=

Transparency;

endif;

if (PLOTCL)

if (PLOTBAND)

bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / group=

STRATUM index=STRATUMNUM modelname="Survival" display=(outline)

;

else

bandplot LimitUpper=SDF_UCL LimitLower=SDF_LCL x=TIME / group=

STRATUM index=STRATUMNUM modelname="Survival" datatransparency=

Transparency;

endif;

stepplot y=SURVIVAL x=TIME / group=STRATUM index=STRATUMNUM name=

"Survival" rolename=(_tip1=ATRISK _tip2=EVENT) tip=(y x Time

_tip1 _tip2);

if (PLOTCENSORED)

scatterplot y=CENSORED x=TIME / group=STRATUM index=STRATUMNUM

markerattrs=(symbol=plus);

endif;

if (PLOTATRISK)

innermargin / align=bottom;

blockplot x=TATRISK block=ATRISK / class="CLASSATRISK"

repeatedvalues=true display=(label values) valuehalign=start

valuefitpolicy=truncate labelposition=left labelattrs=

GRAPHVALUETEXT valueattrs=GRAPHDATATEXT (size=7pt)

includemissingclass=false;

endinnermargin;

endif;

DiscreteLegend "Survival" / title=GROUPNAME location=outside;

if (PLOTCENSORED)

if (PLOTTEST)

layout gridded / rows=2 autoalign=(TOPRIGHT BOTTOMLEFT BOTTOM)

border=true BackgroundColor=GraphWalls:Color Opaque=true;

entry "+ Censored";

if (PVALUE < .0001)

entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4));

else

entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4));

endif;

endlayout;

else

layout gridded / rows=1 autoalign=(TOPRIGHT BOTTOMLEFT BOTTOM)

border=true BackgroundColor=GraphWalls:Color Opaque=true;

entry "+ Censored";

endlayout;

endif;

else

if (PLOTTEST)

layout gridded / rows=1 autoalign=(BOTTOMLEFT TOP BOTTOM)

border=true BackgroundColor=GraphWalls:Color Opaque=true;

if (PVALUE < .0001)

entry TESTNAME " p " eval (PUT(PVALUE, PVALUE6.4));

else

entry TESTNAME " p=" eval (PUT(PVALUE, PVALUE6.4));

endif;

endlayout;

endif;

endlayout;

endif;

EndGraph;

end;

run;

Re: Survival Plot with dashed line.

ballardw — Wed, 20 Jun 2012 20:22:09 GMT

Line style by default is controlled by the value of linestyle in the GRAPHDATA1, GRAPHDATA2 etc. of the style. It may be easiest to modify that parameter in a custom style inherited from your preferred style.

Re: Survival Plot with dashed line.

DanH_sas — Thu, 21 Jun 2012 00:23:44 GMT

Add LINEATTRS=(pattern=dashed) to this line:

stepplot y=SURVIVAL x=TIME / name="Survival" rolename=(_tip1=

ATRISK _tip2=EVENT) tip=(y x Time _tip1 _tip2) legendlabel=

"Survival" lineattrs=(pattern=dashed);

There are two stepplots in this template, do NOT change the one with the GROUP= option. Let me know if this works for you.

Thanks!
Dan

Re: Survival Plot with dashed line.

cypher85 — Thu, 21 Jun 2012 12:11:23 GMT

Worked Great!! Thanks!

topic Re: Survival Plot with dashed line. in Graphics Programming

Survival Plot with dashed line.

Re: Survival Plot with dashed line.

Re: Survival Plot with dashed line.

Re: Survival Plot with dashed line.