I have read SAS threads on exporting multiple graphs which use the same predictor and response variables, but I need to export multiple graphs into a single JPG file (or EPS or TIF) when a different predictor variable is used each time. In addition, there is a group variable for some but not all graphs. Here is the code to output to a Word file (not allowed by journal to which I am submitting): /*muscle-liver figure without interaction*/ proc glm data=som.pittonly ; class EEFEMALE (ref="1: Yes") prediabetes; options fmtsearch=(SOM.formats); model log_MRMLIVFAT = CRPERCSC D1AGE EEFEMALE ACTMMVPAM prediabetes/ solution clparm tolerance; where ACTMMVPAM >1 and CRPERCSC>0.00001; estimate 'CRPERCSC' CRPERCSC 0.0680996; store contcont; output out=outpred p=p lclm=lower uclm=upper; run; quit; DATA INPLOTDATA; ACTMMVPAM=179.4505; d1age=75; EEFEMALE=0; prediabetes=0; do CRPERCSC = .10 to .50 by .1; OUTPUT; END; RUN; PROC PLM RESTORE=CONTCONT; SCORE DATA=INPLOTDATA OUT=PLOTDATAnoint PRED=FIT UCLM=U LCLM=L ; RUN; proc sort data = PLOTDATAnoint; by CRPERCSC ; run; /*muscle-liver with WC interaction*/ proc glm data=som.pittonly ; class EEFEMALE (ref="1: Yes") high_WC prediabetes; options fmtsearch=(SOM.formats); model log_MRMLIVFAT = CRPERCSC|high_WC D1AGE EEFEMALE ACTMMVPAM prediabetes/ solution clparm tolerance; where high_wc ne . and ACTMMVPAM >1 and CRPERCSC>0.00001; store contcont; output out=outpred p=p lclm=lower uclm=upper; run; quit; DATA INPLOTDATA; ACTMMVPAM=179.4505; d1age=75; EEFEMALE=0; prediabetes=0; do CRPERCSC = 0.1 to 0.5 by 0.1; DO high_WC=0, 1; OUTPUT; END; end; RUN; PROC PLM RESTORE=CONTCONT; SCORE DATA=INPLOTDATA OUT=PLOTDATA PRED=FIT UCLM=U LCLM=L ; RUN; proc sort data = PLOTDATA; by CRPERCSC high_WC; run; proc format; value high_wc 0='Normal' 1='High'; run; /*MFI adjusted figure NO INTERACTION*/ proc glm data= som.pittonly; class EEFEMALE (ref="1: Yes") prediabetes; options fmtsearch=(SOM.formats); model log_MRMLIVFAT = MRMATHMFI D1AGE EEFEMALE ACTMMVPAM prediabetes/ solution clparm tolerance; store contcont; run; DATA INPLOTDATA; ACTMMVPAM=179.4505; d1age=75; EEFEMALE=0; prediabetes=0; do MRMATHMFI = 0.0 to 0.16 by 0.02; OUTPUT; END; RUN; PROC PLM RESTORE=CONTCONT; SCORE DATA=INPLOTDATA OUT=PLOTDATA_MFI PRED=FIT UCLM=U LCLM=L ; RUN; proc sort data = PLOTDATA_MFI; by MRMATHMFI ; run; ODS _ALL_ CLOSE; ods graphics on /noborder /*width= 410 px height=300 px*/ width=4in height=3in imagefmt=tiff; options nocenter nonumber nodate orientation=landscape; ods word columns=2 file = "C:\Users\Figure 1.docx" STARTPAGE=no; /*MUSCLE-LIVER NO INTERACTION*/ proc sgplot data=PLOTDATAnoint noborder noautolegend ; format CRPERCSC percent9.4; band x=CRPERCSC lower=l upper=u / transparency=0.3 ; series x=CRPERCSC y=fit / markers markerattrs=(symbol="circlefilled" size=5) lineattrs=(pattern=1 color=black); /*lineattrs=(thickness=2)*/ series x=CRPERCSC y=u / lineattrs=(pattern=1 color=black); series x=CRPERCSC y=l / lineattrs=(pattern=1 color=black); xaxis label = "Muscle Mass (%)" values = (.10 to .50 by .10) LABELATTRS=(size=12pt family="Arial") valueattrs=(size=12pt family="Arial"); yaxis label = "Log Liver Fat (%)" values = (-5.0 to 0 by 0.50) LABELATTRS=(size=12pt family="Arial") valueattrs=(size=12pt family="Arial"); INSET "P=0.82" "Per-SD Adjusted Beta 95%CI:" "-0.87 (-8.1, 6.9)" / TEXTATTRS = (SIZE=12 family="Arial"); run; /*MFI-LIVER NO INTERACTION*/ proc sgplot data=PLOTDATA_MFI noborder noautolegend ; format MRMATHMFI percent9.4; band x=MRMATHMFI lower=l upper=u / transparency=0.3 ; series x=MRMATHMFI y=fit / markers markerattrs=(symbol="circlefilled" size=5) lineattrs=(pattern=1 color=black); /*lineattrs=(thickness=2)*/ series x=MRMATHMFI y=u / lineattrs=(pattern=1 color=black); series x=MRMATHMFI y=l / lineattrs=(pattern=1 color=black); xaxis label = "Mean Anterior Thigh Muscle Fat Infiltration (%)" values = (0.0 to 0.16 by 0.02) LABELATTRS=(size=12pt family="Arial") valueattrs=(size=11pt family="Arial"); yaxis label = "Log Liver Fat (%)" values = (-5.0 to 0 by 0.50) LABELATTRS=(size=12pt family="Arial") valueattrs=(size=12pt family="Arial"); INSET "P=0.0003" "Per-SD Adjusted Beta 95%CI:" "15.2 (6.8, 24.3)" / TEXTATTRS = (SIZE=12 family="Arial"); run; proc sort data = outpred; by high_WC CRPERCSC ; run; /*MUSCLE-LIVER WITH INTERACTION*/ proc sgplot data=PLOTDATA noborder ; format high_WC high_WC. CRPERCSC percent9.4; label high_WC = "Waist circumference"; band x=CRPERCSC lower=l upper=u / group=high_WC transparency=0.3; styleattrs datacolors=(blue orange ); series x=CRPERCSC y=fit / group=high_WC markers markerattrs=(symbol="circlefilled" size=5 ) lineattrs=(pattern=1 color=black); /*lineattrs=(thickness=2)*/ series x=CRPERCSC y=u / group=high_WC lineattrs=(pattern=1 color=black); series x=CRPERCSC y=l / group=high_WC lineattrs=(pattern=1 color=black); xaxis label = "Muscle Mass (%)" values = (.1 to .5 by .1) LABELATTRS=(size=12pt family="Arial") valueattrs=(size=12pt family="Arial"); yaxis label = "Log Liver Fat (%)" values = (-5.0 to 0 by 0.50) LABELATTRS=(size=12pt family="Arial") valueattrs=(size=12pt family="Arial"); INSET ("(*ESC*){UNICODE P}" = "Interaction P=0.001" ) / TEXTATTRS = (SIZE=12 family="Arial"); run; ods word close;
... View more
