topic Re: how to convert it to a macro to generate multiple graphs in png format? in Graphics Programming

how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Wed, 03 Nov 2021 14:33:59 GMT

Hello Everyone,

I am trying to convert my below code into a macro where it generates multiple png images. I have a question on how to get the png formats?

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 / reset height = 6in width = 8in noborder;
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 valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Years" values=(0 to 5 by 1);
	yaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Survival Probability" values=(0.5 to 1 by 0.1);
	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;

As you see it is currently generating just one rtf file now. I have to generate multiple png files.

Thanks!

Re: how to convert it to a macro to generate multiple graphs in png format?

Kurt_Bremser — Wed, 03 Nov 2021 14:52:43 GMT

Use the OUTPUTFMT= option of the ODS GRAPHICS statement.

Re: how to convert it to a macro to generate multiple graphs in png format?

Kurt_Bremser — Wed, 03 Nov 2021 14:56:02 GMT

RTF embeds the images; if you want separate image files, use ODS HTML or similar.

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Wed, 03 Nov 2021 14:59:50 GMT

If I convert above code into a macro call, will that generate multiple png files?

Re: how to convert it to a macro to generate multiple graphs in png format?

ballardw — Wed, 03 Nov 2021 15:08:03 GMT

What is the rule for "multiple graphs"? How do you identify which graph goes where?

The only thing I can guess might be desirable for individual graphs would be the Stratum variable but you are using that in a group role.

If you use a BY STRATUM; (assumes the data set is sorted by Stratum) you will get one graph for each level of Stratum which would create separate PNG files if the proper ODS Graphics options are set. You would likely want to remove Stratum from the group role where used.

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Wed, 03 Nov 2021 15:10:49 GMT

I have different group data to be analyzed. I'm afraid to take it out.

Re: how to convert it to a macro to generate multiple graphs in png format?

Reeza — Wed, 03 Nov 2021 16:07:14 GMT

1. Does the code posted above do what you want? Before you automate anything it should work as you need. If not, solve that issue first.
2. Once you have 1, show how it would work for two or three iterations if you copied/pasted it. Then post back and we can help you convert it to a macro.

Re: how to convert it to a macro to generate multiple graphs in png format?

ballardw — Wed, 03 Nov 2021 16:25:23 GMT

@RAVI2000 wrote:
I have different group data to be analyzed. I'm afraid to take it out.

So it appears, without any rule of how to create multiple graphs, that creating them actually isn't needed.

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Wed, 03 Nov 2021 17:15:41 GMT

Hello @Reeza ,

Thanks for reaching out. Here are the iteration of graphs and their corresponding outputs. I have highlighted the code where ever I have changed it while producing the graphs.

options nodate nonumber; 
ods escapechar='^';
ods listing close;
ods graphics on;
ods output survivalplot = sj.survplot1;
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 / reset height = 6in width = 8in noborder;
title1 j=c h=12pt font = arial "Overall Survival by SES";

proc sgplot data = sj.survplot1;
		step x=time y=survival /group=stratum name='s' lineattrs=(pattern=solid);
	xaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Years" values=(0 to 5 by 1);
	yaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Survival Probability" values=(0.5 to 1 by 0.1);
	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;


ods listing close;
ods graphics on;
ods output survivalplot = sj.survplot2;
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 race_ses;
run;


ods tagsets.rtf file = "/home/Projects/programs/exhibit_02.rtf";
ods graphics / reset height = 7.5in width = 9in noborder;
title1 j=c h=12pt font = arial "Survival by Race and SES";

proc sgplot data = sj.survplot2;
		step x=time y=survival /group=stratum name='s' lineattrs=(pattern=solid);
	xaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Years" values=(0 to 5 by 1);
	yaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Survival Probability" values=(0.3 to 1 by 0.1);
	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=outside;
run; 
ods tagsets.rtf close;
ods listing;



ods listing close;
ods graphics on;
ods output survivalplot = sj.survplot3;
proc lifetest data = sj.adsl4 method = km plots=s(test) 
		timelist = 0 1 2 5 10 15 20 25 30 conftype = linear;
	time years * death(0);
	strata race_ses_min;
run;

ods tagsets.rtf file = "/home/Projects/programs/exhibit_03.rtf";
ods graphics / reset height = 6in width = 8in noborder;
title1 j=c h=12pt font = arial "Survival by Minority Race and SES";

proc sgplot data = sj.survplot3;
		step x=time y=survival /group=stratum name='s' lineattrs=(pattern=solid);
	xaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Years" values=(0 to 5 by 1);
	yaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Survival Probability" values=(0.3 to 1 by 0.1);
	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;

Also some of the other output formats are png images.

Thanks!

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Wed, 03 Nov 2021 17:14:37 GMT

Hello @ballardw,

I have posted my codes for different graphs and their output. I still do have more graphs with different strata variables.

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Wed, 03 Nov 2021 17:19:08 GMT

I have highlighted other elements that are changing across the code. Not sure why they are not appearing here. Some of them include in and out datasets like survplot1,2,3, height, title, strata variables, keylegend location, values for Y-AXIS etc.,

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Thu, 04 Nov 2021 01:58:28 GMT

Hello @Reeza,

Can you help by taking a look into those codes and output and how can I convert them into macros?

Re: how to convert it to a macro to generate multiple graphs in png format?

Reeza — Thu, 04 Nov 2021 15:03:08 GMT

Nothing is highlighted for me.

Re: how to convert it to a macro to generate multiple graphs in png format?

Reeza — Thu, 04 Nov 2021 15:09:58 GMT

Tutorial on converting a working program to a macro

This method is pretty robust and helps prevent errors and makes it much easier to debug your code. Obviously biased, because I wrote it 🙂 https://github.com/statgeek/SAS-Tutorials/blob/master/Turning%20a%20program%20into%20a%20macro.md


%macro create_graphs(input_data = ,
                                      output_data= ,
                                      output_file=,
                                       variable = );
options nodate nonumber; 
ods escapechar='^';
ods listing close;
ods graphics on;
ods output survivalplot = sj.&outputfile;
proc lifetest data = sj.&input_data method = km plots=s(test) 
		timelist = 0 1 2 5 10 15 20 25 30 conftype = linear;
	time years * death(0);
	strata ses &variable;
run;

ods tagsets.rtf file = "/home/Projects/programs/&output_file..rtf";
ods graphics / reset height = 6in width = 8in noborder;
title1 j=c h=12pt font = arial "Overall Survival by SES";

proc sgplot data = sj.&output_data.;
		step x=time y=survival /group=stratum name='s' lineattrs=(pattern=solid);
	xaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Years" values=(0 to 5 by 1);
	yaxis valueattrs=(size=11.2pt) labelattrs=(size=12.5pt) label="Survival Probability" values=(0.5 to 1 by 0.1);
	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;

%mend;

%create_graphs(input_data = adsl3 , output_data= survplot1, output_file=exhibit_01, variable=ses);
%create_graphs(input_data = adsl4 , output_data= survplot2, output_file=exhibit_02, variable=race_ses);

Probably not everything but enough to get started and give you the idea. Please follow the tutorial linked above to see explanations and details.

This may not work because I have no data to test it on, if it doesn't and you cannot debug it yourself please post your code, log with macro debugging options on (MPRINT, SYMBOLGEN) and explain what is wrong.

Re: how to convert it to a macro to generate multiple graphs in png format?

RAVI2000 — Thu, 04 Nov 2021 16:38:31 GMT

This was a very helpful tutorial who are struggling with macros. I have learnt some new things from you @Reeza. Thank you so much!

Re: how to convert it to a macro to generate multiple graphs in png format?

Reeza — Thu, 04 Nov 2021 16:41:44 GMT

You should also check out this macro. It's old so haven't tested it lately but there's example testing code at the bottom.
https://gist.github.com/statgeek/d3bce2a9e2ef0523db9d