topic Re: How do I zoom in the survival graph with many group for analysis in Graphics Programming

How do I zoom in the survival graph with many group for analysis

RAVI2000 — Fri, 29 Oct 2021 03:58:47 GMT

Hello Everyone,

I have created a survival graph with multiple groups as below.

As you look here it is a bit clustered at the top of the graph.

Is there any procedure or way to zoom the graph? so can analyze properly.

ODS PDF file='survival_Ethnicity white race by SES.pdf';
ods graphics on / imagefmt=png imagename= "Surival for Ethnicity by White Race and SES";
proc lifetest data = adsl3 method=KM plots=s(test) 
		timelist = 0 1 2 5 10 15 20 25 30 CONFTYPE= linear;
	time Years*death(0);
	strata race_ses_white;
title1 "Overall Survival Ethnicity by White Race by SES";
run ;

Re: How do I zoom in the survival graph with many group for analysis

whymath — Fri, 29 Oct 2021 05:13:39 GMT

height= option at ods graphics can modify the height of graph.
And locate the legend at right bottom inside the axis area may helps, too.

Re: How do I zoom in the survival graph with many group for analysis

ballardw — Fri, 29 Oct 2021 09:49:59 GMT

May also want to set the graph to make a wider display using Ods graphics /width = <something>.

Or create a data set and use that data for SGPLOT where you can apply more controls.

Re: How do I zoom in the survival graph with many group for analysis

Ksharp — Fri, 29 Oct 2021 12:28:28 GMT

Check this macro and option "height=7in" :

https://communities.sas.com/t5/SAS-Communities-Library/Kaplan-Meier-Survival-Plotting-Macro-NEWSURV/ta-p/479747

Re: How do I zoom in the survival graph with many group for analysis

RAVI2000 — Tue, 02 Nov 2021 19:08:38 GMT


options nodate nonumber; 
ods escapechar='^';
ods listing close;
ods graphics on;
ods output survivalplot = sj.survplot;
proc lifetest data = sj.adsl3 method = km plots=s(test) 
		timelist = 0 1 2 5 10 15 20 25 30 conftype = linear;
	time years * death(0);
	strata ses;
run;


ods tagsets.rtf file = "/home/Projects/programs/exhibit_01.rtf";
ods graphics on / height = 8in width = 9in;
title1 j=c h=12pt font = arial "Overall Survival by SES";

proc sgplot data = sj.survplot;
		step x=time y=survival /group=stratum name='s' lineattrs=(pattern=solid);
	xaxis label="Years" values=(0 to 5 by 1);
	yaxis label="Survival Probability" values=(0 to 1 by 0.2);
	scatter x=time y=censored / markerattrs=(symbol=plus) name='c';
	scatter x=time y=censored / markerattrs=(symbol=plus) group=stratum;
	keylegend 'c'/location=inside position=topright;
	keylegend 's'/location=inside position=bottomleft;
run; 

ods tagsets.rtf close;
ods listing;

I have given height and width option in ODS but still not getting desired output zoomed.

I still feel the graph is clustered. Any other options @ballardw @Ksharp @whymath ? Thankyou!

Re: How do I zoom in the survival graph with many group for analysis

djrisks — Tue, 02 Nov 2021 19:57:55 GMT

Hi @RAVI2000,

Perhaps you can use the MIN option in the YAXIS statement to make the graph go from 0.4 to 1.

Please see the example below:

options nodate nonumber; 
ods escapechar='^';
ods listing close;
ods graphics on;
ods output survivalplot = sj.survplot;
proc lifetest data = sj.adsl3 method = km plots=s(test) 
		timelist = 0 1 2 5 10 15 20 25 30 conftype = linear;
	time years * death(0);
	strata ses;
run;


ods tagsets.rtf file = "/home/Projects/programs/exhibit_01.rtf";
ods graphics on / height = 8in width = 9in;
title1 j=c h=12pt font = arial "Overall Survival by SES";

proc sgplot data = sj.survplot;
		step x=time y=survival /group=stratum name='s' lineattrs=(pattern=solid);
	xaxis label="Years" values=(0 to 5 by 1);
	yaxis label="Survival Probability" values=(0.4 to 1 by 0.2);
	scatter x=time y=censored / markerattrs=(symbol=plus) name='c';
	scatter x=time y=censored / markerattrs=(symbol=plus) group=stratum;
	keylegend 'c'/location=inside position=topright;
	keylegend 's'/location=inside position=bottomleft;
run; 

ods tagsets.rtf close;
ods listing;

Re: How do I zoom in the survival graph with many group for analysis

JeffMeyers — Tue, 02 Nov 2021 20:14:13 GMT

To add to this you might want to turn off the censor markers (not plot the SCATTER statements above). Your graphs have a lot of censor indicators making the lines into what Dr. Therneau would call "Wooly Bear" lines. If you remove them it will be easier to see the lines for comparison and they aren't adding much information to the plot currently.

Re: How do I zoom in the survival graph with many group for analysis

JeffMeyers — Tue, 02 Nov 2021 20:15:27 GMT

Nice to see a shout out to my macro. @RAVI2000 if you want to try using it feel free to ask me any questions you may have.

Re: How do I zoom in the survival graph with many group for analysis

RAVI2000 — Tue, 02 Nov 2021 20:36:39 GMT

I wonder I cannot remove the censor as this is a Spatial data analysis. I have a lot of data. Please correct me.

Re: How do I zoom in the survival graph with many group for analysis

JeffMeyers — Tue, 02 Nov 2021 20:49:27 GMT

I don't think they add much information to the plot and make the lines hard to read. Another option you can do is to add a Number of cumulative censors to the bottom of the graph (like a patient's at risk type of table) if you want to keep track of where and how many censors have occurred.

Re: How do I zoom in the survival graph with many group for analysis

Ksharp — Wed, 03 Nov 2021 11:49:34 GMT

Maybe could try this :

yaxis label="Survival Probability" values=(0.5 to 1 by 0.1);

Re: How do I zoom in the survival graph with many group for analysis

Ksharp — Wed, 03 Nov 2021 11:53:04 GMT

Your macro is awesome . The problem I meet is can't display Chinese . and I have to use option FONTFAMILY= 宋体 to display Chinese characters .

And I need change "+ Censored" in legend into Chinese "+ 删失" .and I have to modify your macro original code to get that .

Re: How do I zoom in the survival graph with many group for analysis

RAVI2000 — Wed, 03 Nov 2021 15:14:55 GMT

Hello @djrisks,

Can you convert this into a macro to generate multiple outputs? meaning how can I give multiple destination files?

Thanks!