Graphics Programming topics

Axis table overlapping the first value

new510 — Wed, 21 Jul 2021 23:26:37 GMT

Hi All,

I am trying to add axis table, and the label(stratum) overlaps the first value. I am using the following code for the axis table.

layout overlay / walldisplay=none xaxisopts=(display=none OFFSETMIN=0);
axistable x=TATRISK value=ATRISK /title="Number At Risk" display=(label) valueattrs=(size=8pt) class=stratum colorgroup=stratum VALUEHALIGN=right;
endlayout;

Example fig:

The label was supposed to be Treat11, and the first value is truncated to 41 instead of 2841.

I really appreciate any help you can provide.

a percentage-based bar graph

dustychair — Sun, 18 Jul 2021 05:52:44 GMT

Hi all,

I have a categorical variable in my data. I am trying to draw a bar graph for this variable based on the percents. The below code is working but the sum of percentages is not 100%. I don't know how it is calculating this. Is there an easy way to draw a percentage-based bar graph?

Thank you,

proc sgplot data=alt_listening;
title "Proficiency Level: List 1-2";
vbar listening_pl / stat=pct transparency=0.25 fillattrs=(color=black) ;
yaxis label= 'Percent' ;
xaxis label= 'Proficiency Level';
by cluster_listening;
run;

Fun With SAS ODS Graphics: Happy Ba(r)stille Day!

tc — Wed, 14 Jul 2021 03:30:42 GMT

Flag of France - Take 3!

A quick SAS ODS Graphics take on the Flag of France for Bastille Day - third time's the charm!

* Fun With SAS ODS Graphics: Happy Ba(r)stille Day! (vertical bar chart Flag of France)
  Based on colors/dimensions at en.wikipedia.org/wiki/Flag_of_France;

data flag; do c=1 to 3; output; end;                                     * Generate flag data - tricolour of blue/white/red;

ods graphics / reset height=4in width=6in antialias;                     * 2:3 proportion (height:width);

proc sgplot data=flag;                                                   * Plot #1: Generate bar chart;
vbar c;

proc sgplot data=flag;                                                   * Plot #2: Generate bar chart with colors of Flag of France;
vbar c / colorresponse=c colormodel=(cx0055A4 cxFFFFFF cxEF4135);

proc sgplot data=flag aspect=.66667 noautolegend pad=0 noborder nowall;  * Plot #3: Generate bar chart as tricolour Flag of France;
vbar c / barwidth=1 colorresponse=c colormodel=(cx0055A4 cxFFFFFF cxEF4135) nooutline baselineattrs=(thickness=0);
xaxis display=none min=.5 max=3.5 offsetmin=0 offsetmax=0 type=linear;   * Suppress axis/padding;
yaxis display=none min=0 max=1 offsetmin=0 offsetmax=0;
run;

Flag of France - Take 1!Flag of France - Take 2!

Create and IMAGE file of output from PROC REPORT

data_null__ — Fri, 09 Jul 2021 20:54:16 GMT

Can I create output from PROC REPORT as an image (PNG). I want to annotate a figure (SGPLOT) with a small table of statistics, using IMAGE function seems like an easy way, but I also don't want to create the image file by hand. I don't think INSET or AXISTABLE are suitable that's why I think IMAGE would be nice.

How to make tickstyle=inside work for BOTH major and minor tick marks?

JonathanNitzan — Fri, 09 Jul 2021 17:00:26 GMT

1. I'm using proc template, trying to have my tick marks plotted inside rather than outside. I'm using option tickstyle=inside, but this option seems to affect only the major tick marks, not the minor ones.

2. I'm using SAS 9.4 1M7.

3. Any suggested solution for how to have BOTH major and minor tick marks inside?

Thank you.

How to reduce spacing in Shewhart chart?

annaabimbola — Thu, 08 Jul 2021 23:05:25 GMT

Hi All,

I'm want to create an XSchart.

The chart was split into 4 instead of 1. Please how do I make this come out as just 1 chart?

The data points for the subgroup variable are 09/01/2020, 10/01/2020,...01/01/2021

Please see the code.

title "Length of Stay";
symbol v=dot;

proc shewhart data=work.nowsbaseline_sorted;
	xschart neonatald*measureperiodid/haxis= '01SEP2020'd to '01MAR2021'd by 31;
run;

Any help is greatly appreciated.

Thanks in advance!

combine Barchart with Circles

pmksasbi — Tue, 06 Jul 2021 21:22:40 GMT

HI ,

can anyone please give some example code for the vertical bar graph with circles plotting on it.

Thanks,

Mahesh

Fun With SAS ODS Graphics: Patriotic Pom Pom July 4th Wreath (Bubble+Scatter Plots)

tc — Sun, 04 Jul 2021 03:52:02 GMT

Patriotic Pom Pom Wreath

A SAS ODS Graphics bubble plot take on a neat Etsy Patriotic Pom Pom Wreath. Have a Happy 4th, all!

* Fun With SAS ODS Graphics: Patriotic Pom Pom July 4th Wreath (Bubble+Scatter Plots)
  Inspired by Patriotic Pom Pom wreath at etsy.com/listing/964348987/patriotic-pom-pom-wreath-red-white-and;

proc format;                              * Assign red/white/blue colors based on angles of wedges; 
value rwb -250--140='Blue' -270--260, 80-90, 20-40, -40--20, -100--80 ='Red'  
          50-70, -10-10, -70--50, -130--110='White'; 

data wreath;                              * Generate points for July 4th 'wreath';
retain size 1 txt 'HAPPY*4TH OF JULY!';
do a=90 to -270 by -10;                   * Generates points for balloons;
  do r=.7 to 1 by .15;                    * Three points for each angle;  
    x=r*cos(a*constant('pi')/180);
    y=r*sin(a*constant('pi')/180);        
    c=put(a,rwb.);                        * Assign color;
    output;
  end;
end;
x=.; y=.;                                 * Generate points for stars;
do a=-150 to -240 by -15;                 * In the blue wedge of wreath;                
  if r=.9 then r=.8; else r=.9;           * Alternate location of stars (radius=.8 or .9);          
  xs=r*cos(a*constant('pi')/180);
  ys=r*sin(a*constant('pi')/180);         
  output;
end;
xs=.; ys=.;                               * Generate point for text;
xt=0; yt=0;
output;
                                          * Patriotic Pom Pom Wreath (Bubble+Scatter Plots);
proc sgplot data=wreath aspect=1 noautolegend noborder;
styleattrs datacolors=(red white navy);
symbolchar name=uniStar char='2605'x;     * Unicode 5-pointed star;
bubble x=x y=y size=size /                /* Red, white and blue bubbles */
       dataskin=pressed group=c nooutline bradiusmin=.23in bradiusmax=.25in;   
scatter x=xs y=ys /                       /* White stars */                          
        markerattrs=(symbol=unistar size=36pt color=white);        
xaxis display=none; yaxis display=none;   * Look Ma, no axes!;
text x=xt y=yt text=txt / strip textattrs=(size=18pt color=navy weight=bold) 
     splitchar='*' splitpolicy=splitalways contributeoffsets=none; * Happy 4th of July!; 
run;

xaxis marks correlated with x-values

Jason2020 — Thu, 01 Jul 2021 14:06:08 GMT

Greetings,

I am using proc sgplot to model the relationship between bone level and follow up years in two treatment groups (A/B). The issue I have is that the tick marks of the x-axis are not spaced right. In other words, the 1 year and 6 years are spaced equally. I hope someone can show me how to space the marks so that they are correlated with the number of years. I included below a sample of my data, the program I used, and the resulting figure.

Thank you.

subject group followup level

02 A 1.0 4.7

02 A 4.5 4.6

02 A 6.0 4.9

02 A 0.0 1.6

03 A 2.0 2.3

03 A 6.0 1.4

03 A 0.0 1.5

07 B 1.0 1.1

07 B 2.0 1.4

07 B 4.5 1.3

07 B 6.0 1.3

07 B 0.0 1.1

08 A 2.0 1.0

08 A 4.5 1.0

proc sgplot data=p1;

vline followup / response=level group=group

stat=mean limitstat=stderr lineattrs=(thickness=6px);

yaxis label="Bone Level (mm)"

valueattrs=(size=18pt) labelattrs=(size=20pt weight=bold);

xaxis discreteorder=unformatted label="Follow up (years)"

valueattrs=(size=18pt) labelattrs=(size=20pt weight=bold);

keylegend /location=inside valueattrs=(size=18pt)

title="" titleattrs=(size=16pt);

run;

How to change font in a plot

marta25 — Thu, 01 Jul 2021 07:04:38 GMT

Hello,

It seems a very easy question but somehow I am struggling with this.

I just want to change the FONT (not font size, this I can do easily) of all text in a plot (in SGPLOT): axis labels and values, legend items...

I guess I have to define a new style, something like:

proc template;
define style my_style;
parent=styles.htmlblue;
class GraphValueText / font = ("<sans-serif>, <MTsans-serif>",20pt);
end;

How can I set the font to Times New Roman, or Arial, etc, instead of sans-serif?

Where can I find a list of available fonts?

And how do I apply the new style to the plot?

Could someone give an example (as easy as possible)?

proc lifetest not all tatrisk values are not being displayed.

S_RAVI — Thu, 01 Jul 2021 21:10:34 GMT

I have used proc lifetest and output the survival plot for atrisk and tatrisk to plot the km plot. My survival output dataset has tatrisk values as 0, 2, 4, 6. But when I output it using proc template and sgrender only 0, 2, 6 values are being displayed.

using annotate in proc sgplot

Anita_n — Tue, 29 Jun 2021 13:01:35 GMT

Dear all,

I have the following code to draw a bar chart and label each bar with the sum total of each bar or each year using an annotate dataset. The seglabels should be in percentage as shown in the graph, since I can't use datalabel here, I creatde an annotate dataset but can't get it working, I might probably be missing something somewhere. any help? My code is as below

data mydat;
infile datalines;
input year	level $12.	number_of_cases;
datalines;
2001	1	        90	
2001	2	        1	 
2001	3	        5	 
2001	4	        22	
2001	no cases	35	
2002	1	        99	
2002	2	        4	 
2002	3	        2	 
2002	4	        55	
2002	no cases	39	
2003	1	        200
2003	2	        45	
2003	3	        77	
2003	4	        39	
2003	no cases	23	
2004	1	        205
2004	2	        19	
2004	3	        29	
2004	4	        15	
2004	no cases	18	
2005	1	        159
2005	2	        323
2005	3	        99	
2005	4	        10	
2005	no cases	28	
2006	1	        88	
2006	2	        10	
2006	3	        16	
2006	4	        55	
2006	no cases	8	 
;
run;
*calculate sum of cases per year;
proc means data=mydat noprint nway;
	class year;
	var number_of_cases;
	output out=sum_mydat (drop=_TYPE_ _FREQ_) sum= / autoname;
run;
*calculate rel. perecntages;
data merge_mydat;
	merge mydat sum_mydat;
	by year;
	format rel_freq percent8.0; 
	rel_freq = divide(number_of_cases,number_of_cases_sum);
    drop number_of_cases_sum;
run;

*create annodat;
data anno;
length function color $8 style $12;                                                                                                  
retain function 'label' when 'a' size 1.25 position '2'                                                                              
          color 'black' xsys '2' ysys '2' style 'Albany AMT' weight 'bold';    
set sum_mydat;

midpoint=year;                                                                                                                       
   y=number_of_cases_Sum;                                                                                                                            
   text=put(number_of_cases_Sum,3.);  

run;

*create attrmap;
data myattrmap;
	infile datalines;
	input id $6. value $10. fillcolor $8. linecolor $9.;
	datalines;
level no cases  CXA9A9A9 CXA9A9A9
level 1         CX00FFFF CX00FFFF
level 2         CXFF7F50 CXFF7F50
level 3         CXFFF8DC CXFFF8DC
level 4         CX006400 CX006400
;
run;

*create chart with labels, the totals of each bar should be written on top of the bar;

proc sgplot data =merge_mydat dattrmap=myattrmap sganno=anno pad=(bottom=5%) noautolegend noborder;
xaxis grid type=discrete discreteorder=data;
vbarparm category= year  response=rel_freq /group=level groupdisplay=stack attrid=level seglabel seglabelattrs=(size=0.25 color="CXFDF5E6") 
seglabelfitpolicy=none barwidth=0.6 missing;
xaxis  display= (nolabel) valueattrs=(family="verdana" size=8pt) fitpolicy=rotate ValuesRotate=Vertical;
yaxis  display= (nolabel) valueattrs=(family="verdana" size=8pt);
keylegend / location=outside position=bottomright across=0 title=" " valueattrs=(size=2 family="verdana") outerpad=(top=0.2cm) noborder;
run;

Proc SGPANEL color by subject

csa — Tue, 22 Jun 2021 17:12:40 GMT

Hi All,

I am trying to generate a figure plot using the below code. I am getting the output as in the attached screen shot.

proc sgpanel data=have ;
panelby SUBJID / layout=PANEL novarname columns=4 rows=2;
series x=x_axis_var y=aval /group=trt01p_ attrid=myid lineattrs=(pattern=mediumdash thickness=2);
scatter x=x_axis_var y=aval /group=trt01p_ attrid=myid markerattrs = (symbol=circlefilled size=8);

rowaxis integer grid values=(0,50,75,100) label="test";
colaxis integer grid values=(1 to 10 by 1) discreteorder=data label="Day";

run;
quit;

My query is how can we change the Blue color plots 'Green', Red color plots to 'Blue' and the bottom Legend with Green and Blue colors?

Thanks.

DATTRMAP for data label text color

snoopy369 — Wed, 23 Jun 2021 22:13:20 GMT

Is it possible in SGPIE to use a dynamic data label text color? If so, what is the right variable name?

For example:

data attrmap;
  set sashelp.class;
  if sex='M' then textcolor='blue';
  else textcolor='pink';
  fillcolor='orange';
  id = 'mycolors';
  value = name;
  keep id value textcolor fillcolor;
run;
  
  
proc sgpie data=sashelp.class dattrmap=attrmap;
pie name/attrid=mycolors;
run;

I know that the attribute map is "working" for the rows since the fillcolor is changing, but the textcolor won't change, and I have tried everything I can - textcolor, datalabeltextcolor, datalabelattrs with 'color=blue' - but can't seem to find the option (or, it's not possible?). Don't really want to go the SGANNO route as suggested on a somewhat similar post, so perhaps it's just not possible... Thanks!

ODS LAYOUT ABSOLUTE with a BY statement in PROC SGPLOT

PaigeMiller — Wed, 23 Jun 2021 12:27:02 GMT

I can create an output with four plots on one page using ODS LAYOUT ABSOLUTE, as shown in code #1. But then I want to use a BY statement in PROC SGPLOT, such that the first I get four plots on one page, where all four plots are for the same level of the BY variable and the plots are always in the same position, see code #2.

Code #1, this works as expected, but I have not used a BY statement:

ods pdf file="absolute.pdf";
ods layout absolute ;
ods region x=.5in y=0.1in width=4in height=2.75in;
proc sgplot data=sashelp.cars;
    histogram msrp;
run;
ods region  x=5.1in y=.1in width=4in height=2.75in;
proc sgplot data=sashelp.cars;
    histogram invoice;
run;
ods region  x=.5in y=5.25in width=4in height=2.75in;
proc sgplot data=sashelp.cars;
    histogram mpg_city;
run;
ods region  x=5.1in y=5.25in width=4in height=2.75in;
proc sgplot data=sashelp.cars;
    histogram horsepower;
run;
ods layout end;
ods pdf close;

Now, code #2, with a BY Statement, it doesn't work at all, how can I get something like this to work? I want three output pages, the first page would be four plots for origin='Asia', with MSRP top left, INVOICE top right, MPG_CITY bottom left and HORSEPOWER bottom right; second page would be for origin='Europe', exact same arrangement of plots; third page would be for origin='USA', exact same arrangement of plots.

proc sort data=sashelp.cars out=cars;
    by origin;
run;
ods pdf file="absolute2.pdf";
ods layout absolute ;
ods region x=.5in y=0.1in width=4in height=2.75in;
proc sgplot data=cars;
    by origin;
    histogram msrp;
run;
ods region  x=5.1in y=.1in width=4in height=2.75in;
proc sgplot data=cars;
    by origin;
    histogram invoice;
run;
ods region  x=.5in y=5.25in width=4in height=2.75in;
proc sgplot data=cars;
    by origin;
    histogram mpg_city;
run;
ods region  x=5.1in y=5.25in width=4in height=2.75in;
proc sgplot data=cars;
    by origin;
    histogram horsepower;
run;
ods layout end;
ods pdf close;

Adding: I'm sure I could just wrap this in a macro and loop through all possible values of the BY variable and do it that way, but over the years I have found that SAS is very good at anticipating needs and creating features to do this, so really I'd like to do this without a macro if possible.

Proc SGPANEL color by subject

csa — Tue, 22 Jun 2021 17:41:31 GMT

Hi All,

I am trying to generate a figure plot using the below code. I am getting the output as in the attached screen shot.

proc sgpanel data=have ;
panelby SUBJID / layout=PANEL novarname columns=4 rows=2;
series x=x_axis_var y=aval /group=trt01p_ lineattrs=(pattern=mediumdash thickness=2);
scatter x=x_axis_var y=aval /group=trt01p_ markerattrs = (symbol=circlefilled size=8);

rowaxis integer grid values=(0,50,75,100) label="test";
colaxis integer grid values=(1 to 10 by 1) discreteorder=data label="Day";

run;
quit;

My query is how can we change the Blue color plots 'Green', Red color plots to 'Blue' and the bottom Legend with Green and Blue colors?

SAS GTL: Symbol (created by SymbolChar statement) moving on Y-Axis when increasing size

mborne — Thu, 10 Jun 2021 10:27:55 GMT

Hi all, I have an issue using SymbolChar statement to define MARKERATTRS symbol value for a scatterplot.

When I increase the size of the symbol in MARKERATTRS, then the symbol is moving on the Y-Axis. it looks like the symbol has something like POSITION=BOTTOM whereas I'm expecting MIDDLE or CENTER value. I can try to decrease the impact by adding vertical offset to my symbol. But if anybody has an idea of what is happening and how to avoid this it would be great for me (I'm using SAS 9.4 9.4 (TS1M5)).

Thank you for reading.
Mickael

data WORK.PLOT ;
infile datalines4 dsd missover dlm=',' ;
length DAY GEOMEAN LCLM UCLM 8. ;
input DAY GEOMEAN LCLM UCLM ;
datalines4 ;
1,5.851,5.721,5.983
15,10.103,8.531,11.965
29,17.723,14.523,21.628
57,25.449,20.171,32.108
71,76.898,62.926,93.972
85,119.914,101.658,141.447
;;;;
run ;

** Style to be used for the plot ** ;
/*this SAS style can be used: STYLES.Journal;*/
ods path (PREPEND) WORK.GRAPH ;

proc template ;
define statgraph TestSymbolChar / store=WORK.GRAPH ;
mvar GraphTitle EmDashSize ;
dynamic XVar YMean YLowBound YUpBound ;
begingraph ;
symbolchar name=EmDash char='2014'x / /*textattrs=(size=20pt)*/ hoffset=0 voffset=0 /*+0.1*/ ;
entrytitle textattrs=(size=10pt weight=bold) GraphTitle ;
layout overlay / xaxisopts=(label='Day')
yaxisopts=(griddisplay=on type=log logopts=(base=10 Thresholdmax=1 Thresholdmin=1 minorticks=true minortickcount=8 minorgrid=true)
label='Geometric Mean (95%CI)' griddisplay=on) ;
scatterplot x=XVar y=YMean / markerattrs=(symbol=EmDash size=EmDashSize weight=bold color=black) errorbarcapshape=serif
errorbarattrs=(thickness=1pt color=black) yerrorlower=YLowBound yerrorupper=YUpBound name="GeoMean" ;
endlayout ;
endgraph ;
end ;
quit ;
*** Path of the template catalogs *** ;
ods path SASUSER.TEMPLAT(UPDATE) SASHELP.TMPLMST(READ) WORK.MYTEMP(READ) WORK.GRAPH(READ) ;

** SAS SGRender ** ;
%let EmDashSize= ;
%let GraphTitle=Using &EmDashSize. as a size for Symbol Char ;
%let FigureRoot=C:\Temp ;
option center ;
ods listing close ;
ods graphics on / reset=index imagename="TestSymbolCharEmpty" height=10cm width=16cm noborder ;
ods html style=styles.journal2 image_dpi=150 gpath="&FigureRoot." path="&FigureRoot." (url=none) file="output.html" ;

*goptions reset=all ;
ods graphics / outputfmt=png ;

proc sgrender data=WORK.PLOT template=TestSymbolChar description='' ;
dynamic XVar='DAY' YMean='GEOMEAN' YLowBound='LCLM' YUpBound='UCLM' ;
run ;

ods html close ;
ods graphics off ;
ods listing ;

** SAS SGRender ** ;
%let EmDashSize=12pt ;
%let GraphTitle=Using &EmDashSize. as a size for Symbol Char ;
%let FigureRoot=C:\Temp ;
option center ;
ods listing close ;
ods graphics on / reset=index imagename="TestSymbolChar12pt" height=10cm width=16cm noborder ;
ods html style=styles.journal2 image_dpi=150 gpath="&FigureRoot." path="&FigureRoot." (url=none) file="output.html" ;

*goptions reset=all ;
ods graphics / outputfmt=png ;

proc sgrender data=WORK.PLOT template=TestSymbolChar description='' ;
dynamic XVar='DAY' YMean='GEOMEAN' YLowBound='LCLM' YUpBound='UCLM' ;
run ;

ods html close ;
ods graphics off ;
ods listing ;

%let EmDashSize=18pt ;
%let GraphTitle=Using &EmDashSize. as a size for Symbol Char ;
%let FigureRoot=C:\Temp ;
option center ;
ods listing close ;
ods graphics on / reset=index imagename="TestSymbolChar18pt" height=10cm width=16cm noborder ;
ods html style=styles.journal2 image_dpi=150 gpath="&FigureRoot." path="&FigureRoot." (url=none) file="output.html" ;

*goptions reset=all ;
ods graphics / outputfmt=png ;

proc sgrender data=WORK.PLOT template=TestSymbolChar description='' ;
dynamic XVar='DAY' YMean='GEOMEAN' YLowBound='LCLM' YUpBound='UCLM' ;
run ;

ods html close ;
ods graphics off ;
ods listing ;

%let EmDashSize=40 ;
%let GraphTitle=Using &EmDashSize. as a size for Symbol Char ;
%let FigureRoot=C:\Temp ;
option center ;
ods listing close ;
ods graphics on / reset=index imagename="TestSymbolChar4P0" height=10cm width=16cm noborder ;
ods html style=styles.journal2 image_dpi=150 gpath="&FigureRoot." path="&FigureRoot." (url=none) file="output.html" ;

*goptions reset=all ;
ods graphics / outputfmt=png ;

proc sgrender data=WORK.PLOT template=TestSymbolChar description='' ;
dynamic XVar='DAY' YMean='GEOMEAN' YLowBound='LCLM' YUpBound='UCLM' ;
run ;

ods html close ;
ods graphics off ;
ods listing ;

proc template/ proc sgrender graph

S_RAVI — Wed, 16 Jun 2021 23:42:44 GMT

I was trying for a Cumulative Incidence graph. Please correct me.


options nodate nonumber; 
data sj.data;
infile datalines dsd truncover;                                                                                                                                                          
input @1 SEQ 2. @4 ARMS $6. @10 armcode 2. @12 event 1. @14  mon 15. ;
datalines;					
01 Dose1 1 0 3.8107752956636
02 Dose1 1 0 2.66097240473062
03 Dose1 1 0 2.43101182654402
04 Dose1 1 0 2.20105124835742
05 Dose1 1 0 1.44546649145861
;

ods output survivalplot = kmdata;
proc lifetest data = test;
	time mon * event(1);
	strata arms;
run;

data kmdata2;
	set kmdata;
	survival2 = (1-survival);
	censor2 = (1-censored);
run;

proc template;
	define statgraph SurvivalPlot;
	begingraph / designwidth=8in designheight=5in;
	entrytitle 'Cumulative Incidence Graph'/ textattrs=(size=12);
	layout lattice;
	layout overlay / xaxisopts =(Label = "Months (Time)" type=linear labelattrs=(size=8pt weight=bold) tickvalueattrs=(size=9pt weight=bold) linearopts=(viewmin=0 viewmax=12 tickvaluelist=(0 2 4 6 8 10 12)))
		yaxisopts=(Label = "Cumulative Incidence" linearopts=(viewmin=0 viewmax=100) labelattrs=(size=8pt weight=bold) tickvalueattrs=(size=9pt weight=bold) linearopts=(viewmin=0 viewmax = 1.0 tickvaluelist=( 0.0 0.2 0.4 0.6 0.8 1.0)));
	stepplot x=mon y = survival2 / group=stratum  lineattrs=(thickness = 3.0 pattern=solid color = orange) name='s';
	scatterplot x=mon y = censor2 / GROUP=stratum name='c' markerattrs=(symbol=plus size=8);
	mergedlegend "c" "s" / location=inside halign=right valign=top order=columnmajor down=2 border=false valueattrs=(size=8pt) pad=(top=5 bottom=5);
	endlayout;
	layout overlay;
	entry halign=left 'N at Risk' /pad=(left=0 ) valign=bottom textattrs=(size=8pt);
	endlayout;
	blockplot x=tatrisk block=atrisk / class=stratum display=(values label) valuehalign=start valueattrs=(size=8) labelattrs=(size=8);
	endlayout;
	endgraph;
	end;
run;

proc sgrender data = kmdata2 template = survivalplot;
run;

Here are warnings which are shown when I run my codes.

376        options nodate nonumber;
377        ods graphics on / noborder;
378        ods output survivalplot = kmdata;

381        proc lifetest data = test;
382        	time prog_mon * prog_event(1);
383        	strata arms;
384        run;

NOTE: The LOGLOG transform is used to compute the confidence limits for the 
      quartiles of the survivor distribution. To suppress using this transform, 
      specify CONFTYPE=LINEAR in the PROC LIFETEST statement.
NOTE: Variable Stratum already exists on file KMDATA, using Stratum2 instead.
NOTE: The homogeneity tests across strata were omitted since there is only one 
      stratum.
NOTE: The data set KMDATA has 6 observations and 8 variables.
NOTE: The PROCEDURE LIFETEST printed page 1.
NOTE: PROCEDURE LIFETEST used (Total process time):
      real time           2.41 seconds
      cpu time            0.11 seconds
      

385        
386        data kmdata2;
387        	set kmdata;
388        	survival2 = (1-survival);
389        	censor2 = (1-censored);
390        run;

NOTE: Missing values were generated as a result of performing an operation on 
      missing values.
      Each place is given by: (Number of times) at (Line):(Column).
      6 at 389:14   
NOTE: There were 6 observations read from the data set KMDATA.
NOTE: The data set KMDATA2 has 6 observations and 10 variables.
NOTE: DATA statement used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds
      

391        
392        proc template;
393        	define statgraph SurvivalPlot;
394        	begingraph / designwidth=8in designheight=5in;
395        	entrytitle 'Cumulative Incidence Graph'/
395      ! textattrs=(size=12);
396        	layout lattice;
397        	layout overlay / xaxisopts =(Label = "Months (Time)"
397      ! type=linear labelattrs=(size=8pt weight=bold)
397      ! tickvalueattrs=(size=9pt weight=bold) linearopts=(viewmin=0
397      ! viewmax=12 tickvaluelist=(0 2 4 6 8 10 12)))
398        		yaxisopts=(Label = "Cumulative Incidence" linearopts=(viewmin=0
398      ! viewmax=100) labelattrs=(size=8pt weight=bold)
398      ! tickvalueattrs=(size=9pt weight=bold) linearopts=(viewmin=0 viewmax =
398      !  1.0 tickvaluelist=( 0.0 0.2 0.4 0.6 0.8 1.0)));
399        	stepplot x=mon y = survival2 / group=stratum
The SAS System

399      ! lineattrs=(thickness = 3.0 pattern=solid color = orange) name='s';
400        	scatterplot x=mon y = censor2 / GROUP=stratum name='c'
400      ! markerattrs=(symbol=plus size=8);
401        	mergedlegend "c" "s" / location=inside halign=right valign=top
401      ! order=columnmajor down=2 border=false valueattrs=(size=8pt)
401      ! pad=(top=5 bottom=5);
402        	endlayout;
403        	layout overlay;
404        	entry halign=left 'N at Risk' /pad=(left=0 ) valign=bottom
404      ! textattrs=(size=8pt);
405        	endlayout;
406        	blockplot x=tatrisk block=atrisk / class=stratum display=(values
406      ! label) valuehalign=start valueattrs=(size=8) labelattrs=(size=8);
407        	endlayout;
408        	endgraph;
409        	end;
NOTE: Overwriting existing template/link: SurvivalPlot
NOTE: STATGRAPH 'SurvivalPlot' has been saved to: SASUSER.TEMPLAT
410        run;
NOTE: PROCEDURE TEMPLATE used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds
      

411        
412        proc sgrender data = kmdata2 template = survivalplot;
413        run;

WARNING: The stepplot statement named 's' will not be drawn because one or more 
         of the required arguments were not supplied.
WARNING: The scatterplot statement named 'c' will not be drawn because one or 
         more of the required arguments were not supplied.
WARNING: The blockplot statement will not be drawn because one or more of the 
         required arguments were not supplied.

WARNING: The stepplot statement named 's' will not be drawn because one or more 
         of the required arguments were not supplied.
WARNING: The scatterplot statement named 'c' will not be drawn because one or 
         more of the required arguments were not supplied.
WARNING: The blockplot statement will not be drawn because one or more of the 
         required arguments were not supplied.
NOTE: There were 6 observations read from the data set KMDATA2.
NOTE: PROCEDURE SGRENDER used (Total process time):
      real time           0.19 seconds
      cpu time            0.06 seconds

Currently I am not getting anything in my graph.

Please let me know what exactly I am missing here

Can PROC GINSIDE test for points within overlapping polygons?

genemroz — Mon, 14 Jun 2021 17:46:02 GMT

Recently ballardw helped me with a question about how to test whether points lie within a polygon using PROC GINSIDE.

In the code below, I have modified the example code he provided to include a third polygon that overlaps the two he created for his example. Now, two points, #4 and #6, each lie within two polygons. Point #4 lies within polygons 1 and 3. Point # 5 lies within polygons 2 and 3. PROC GINISDE reports that these two points lies within polygon 1 and 2, respectively, but does not report that they are also within polygon 3.

How can I modify the code so that PROC GINISDE correctly reports all polygons occupied by a point? If possible, I would like to retain all the polygons within a single map file.

Thanks in advance for any help,

Gene

data examplemap;
   input id x y;
datalines;
1   0 0
1   0 10
1   5 15
1   8 10
1   8 0
1   0 0
2   12 12
2   12 20
2   18 20
2   18 12
2   12 12
3   3  9
3   3  19
3   16  19
3   16 9
3   3  9
;


data exampledata;
   input otherid x y;
datalines;
1  0  .5
2  1  4
3  0  20
4  4  12
5  10 10
6  13 14
7  20 1
;

proc ginside map =examplemap
             data=exampledata
             out =result
             includeborder
;
   id id;
run;

Regression Plot with Confidence 95% CI Bounds

ugoagwup — Thu, 10 Jun 2021 14:07:24 GMT

Hello!

I need some assistance. A simple linear regression script such as the following would give me a regression plot with 95% CI bounds on defaultRegression Plot.

Proc glm data= Tdata;

model LFTvalues = bckTrans_auc/ solution;
run;

Please how can I use Proc SGplot to replicate exact same plot and have full control of things display like removing the title from the plot?

Unicode symbol in data label attributes statement

sidello — Wed, 09 Jun 2021 10:01:46 GMT

Dear all,

In the vline statement below, I use character variable as datalabel.

This character string now contains "National average"; instead of these characters I would like to use x-bar symbol here.
What is the best way to do this?

I tried to use the unicode in this character variable, but as a result I just get the printed code (string of characters).

First code:

data leavetop;
set leavetop;
x2="National average= *"!!strip(round(x))!!"%";
run;

Attempt with unicode:

data leavetop;
set leavetop;
x2="(*ESC*){U+0304}"!!strip(round(x))!!"%";
run;

In the vline data label attribute statement I also changed the font to "GraphUnicodeText", but this didn't change much as well.

Code of the graph:

proc sgpanel data=leavetop dattrmap=attrmap noautolegend;
panelby _NAME_/
		onepanel
		columns=1
		rows=1
		spacing=7 
		novarname 
		noborder 
		noheaderborder 
		skipemptycells
		uniscale=all
		headerbackcolor=white
		headerattrs=(color=black size=8 weight=bold family=arial);
vbar hosp/ group=countryid groupdisplay=cluster response=col2 attrid=category fill nooutline datalabel=col2 
					datalabelattrs=(color=black size=5 family=arial) DATALABELFITPOLICY=none name="include";

vline hosp /response=x lineattrs=(color=darkgrey pattern=solid thickness=2) group=countryid groupdisplay=overlay splitchar="*"
		datalabel=test datalabelpos=topright datalabelattrs=(size=5 color=black family=GraphUnicodeText) name="dropper" ;

vbar hosp/ group=countryid groupdisplay=cluster response=col2 attrid=category fillattrs=(transparency=1) nooutline datalabel=col2 
					datalabelattrs=(color=black size=5 family=arial) DATALABELFITPOLICY=none name="dropper2";

colaxis display=(noticks noline) label="Hospital" fitpolicy=rotate discreteorder=data offsetmin=0.02 offsetmax=0.02
        valueattrs= (family=arial
	                color=black
	                size=7.5
	                style=normal)
		labelattrs= (family=arial
	                color=black
	                size=9
	                style=normal
	                weight=bold);
rowaxis display=(noticks noline) label="Percent" fitpolicy=thin
		values=(0 to 100 by 20)
		valuesdisplay=("0%" "20%" "40%" "60%" "80%" "100%")
        valueattrs= (family=arial
	                color=black
	                size=7.5
	                style=normal)
        labelattrs= (family=arial
	                color=black
	                size=9
	                style=normal
	                weight=bold);
keylegend "include" /noborder autooutline autoitemsize valueattrs=(size=9 color=black family=arial weight=bold) titleattrs=(color=white) sortorder=ascending;
format _NAME_ $edutop.;
run;

Any suggestions on how to work around this or alternatives are much appreciated.

Kind regards,

Forest plot

Manj — Mon, 07 Jun 2021 17:58:16 GMT

Hi,

I have created a forest plot as below. However I would like to make the below:

1) The CI bars in Green circle to be green in colour.

2) would need an Arrow under the text 'Favours Control Favours TMT ' as shown with the red arrow .

My plot and the code attached belo:

data ip;
length col1 $20 grp $25 adjmn $20;

col1= 'Non-Dialysis';grp= 'Non-Dialysis';grpcd=1;catcd=1;obsid=14;Output;
col1= 'Incident Dialysis';grp= 'Incident Dialysis';grpcd=2;catcd=2;obsid=9;Output;
col1= 'Dialysis';grp= 'Dialysis';grpcd=3;catcd=3;obsid=6;Output;

col1= 'Non-Dialysis';
grp='study001 AXXAA-ND, OL, QD'; grpcd=1;catcd=1;obsid=13;
N_val= 'TMT 1801 Control 1799';
adjmn='0.01 (-0.10,0.11)';
low=-0.1;
high=0.11 ;
mean=0.01 ;
output;

col1= 'Non-Dialysis';
grp=' study001 ESA users'; grpcd=2;catcd=1;obsid=12;
N_val= 'TMT 1001 Control 1013';
adjmn='0.28 (-0.05,0.54)';
low=-0.05;
high=0.54 ;
mean=0.28;
output;

col1= 'Non-Dialysis';
grp=' study001 ESA non-users'; grpcd=3;catcd=1;obsid=11;
N_val= 'TMT 480 Control 470';
adjmn='-0.06 (-0.17,0.26)';
low=-0.17;
high=0.26 ;
mean=-0.06;
output;

col1= 'Incident Dialysis';
grp='study002 AXXAA-ID, OL, QD'; grpcd=4;catcd=2;obsid=8;
N_val= 'TMT 1481 Control 1483';
adjmn='0.23 (-0.01,0.47)';
low=-0.01;
high=0.47;
mean=0.23;
output;

col1= 'Dialysis';
grp='study003 AXXAA-D, OL, QD'; grpcd=5;catcd=3;obsid=5;
N_val= 'TMT 1481 Control 1483';
adjmn='0.23 (-0.01,0.47)';
low=-0.01;
high=0.47;
mean=0.23;
output;

col1= 'Dialysis';
grp=' study003 HD';
N_val= 'TMT 1001 Control 1013';grpcd=6;catcd=3;obsid=4;
adjmn='0.28 (-0.05,0.54)';
low=-0.05;
high=0.54 ;
mean=0.28;
output;

col1= 'Dialysis';
grp=' study003 PD';
N_val= 'TMT 480 Control 470';grpcd=7;catcd=3;obsid=3;
adjmn='-0.06 (-0.17,0.26)';
low=-0.17;
high=0.26 ;
mean=-0.06;
output;

col1= 'Dialysis';
grp='study004 AXXAA-TD, DB, TIW'; grpcd=8;catcd=3; obsid=1;
N_val= 'TMT 204 Control 203';
adjmn='0.23 (-0.01,0.47)';
low=-0.01;
high=0.47;
mean=0.23;
output;

run;

proc sort;by catcd grpcd;run;

data ip;
set ip ;
if obsid in (1,5,6,8,9,13,14) then IndtWt=0;
else indtwt=1;

if obsid in (1,5,6,8,9,13,14 ) then id=1;
else id=2;
run;

data extra;
sid=16; low=-1.5; texta=' <------ Favours control '; high=0.5; textb=' Favours TMT ------>';output;
run;

data ip;
set ip extra;
run;

*** CREATE FIGURE ***;
ods listing close;
ods pdf file="&g_outfile..&g_fontsize" nobookmarkgen style=idsl dpi=600 ;

data attrmap;
length textweight $10;
id='text'; value='1'; textcolor='Black'; textsize=7; textweight='bold'; output;
id='text'; value='2'; textcolor='Green'; textsize=6; textweight='normal'; output;
id='obsid'; value='3'; color='Green'; output;
id='obsid'; value='3'; symbolcolor='Green'; output;
run;
proc sort data=ip out=allnew;by catcd grpcd grp ;run;

proc sgplot data=ip nowall noborder nocycleattrs noautolegend dattrmap=attrmap;
styleattrs axisextent=data ;
highlow y=obsid low=low high=high/highcap=serif lowcap=serif attrid=obsid;
scatter y=obsid x=mean / markerattrs=(symbol=circlefilled size=7) attrid=obsid;
scatter y=obsid x=mean / markerattrs=(size=0) x2axis ;
refline 0 / axis=x lineattrs=(pattern=1);
refline 1.2 / axis=x lineattrs=(pattern=3);

yaxis display=none colorbands=odd colorbandsattrs=(transparency=1) offsetmin=0.02;
x2axis display=(noline noticks novalues nolabel) labelattrs=(size=8);
xaxis values=(-2 -1 0 1 2) labelattrs=(size=8) labelattrs=(size=8);
*TEXT x=low y=sid text=texta / textattrs=(size=8 weight=normal) ;
* TEXT x=high y=sid text=textb / textattrs=(size=8 weight=normal) ;

x2axis label='<-------Favours Control Favours tmt -------->' display=(noline noticks novalues) labelattrs=(size=9 weight=bold);

yaxistable grp / location=inside position=left labelattrs=(size=9 weight=bold) valueattrs=(size=9 ) textgroup=id textgroupid=text indentweight=IndtWt label=' ' labelhalign=left;
yaxistable N_val / location=inside position=left labelattrs=(size=9 weight=bold) valueattrs=(size=9 ) label='N ' textgroup=id textgroupid=text labelhalign=center;
yaxistable adjmn / location=inside position=left pad=(left=15px right=15px) labelattrs=(size=9 weight=bold) textgroup=id textgroupid=text label='Adj. Mean Hgb difference (95% CI) ' labelhalign=left valueattrs=(size=9);

run;

ods pdf close;
ods listing;

Create a polygon and test whether point(s) lie within?

genemroz — Wed, 02 Jun 2021 20:06:23 GMT

I need to create a polygon and test whether a number of points lie within. I've created the polygon (from lat/long coordinates transformed to a Cartesian xy grid) using the POLYGON function available in the SG Annotation Function Dictionary. But I don't know how to test whether a x,y point lies within or without. Will PROC GINSIDE accept a polygon generated by the Polygon Function? It seems dubious given that the annotation function is about drawing and not about "mapping" as seems to be a requirement for PROC GINSIDE. Looking for some guidance on how best to proceed.

Thank you,

Gene

PROC template Graph

Adeoye — Thu, 03 Jun 2021 15:27:48 GMT

I have attached the data if it would be useful in your helping me achieve this task.

Data Selected; 
	Set mylibrary.Mean_Comparison_Proc;
	Period2 = input(Period, 2.);
	Drop Period;
	rename Period2=Period;
run;
Proc Sort Data=Selected;
	by resp_name Treatment Period;
run;

I created the above individual graphs using SGPLOT within a macro, with the below code:

Proc Sql noprint;
	select count(distinct resp_name) into :obs from Selected;
	select distinct resp_name into :resp1-:resp5 from Selected;
Quit;

%Macro TestSee;
	%do i = 1 %to &obs;
		%Let value1=(650 to 750 by 10); %Let value2=(1.0 to 2.0 by 0.1); %Let value3=(1 to 2.0 by 0.1); %Let value4=(2 to 5 by 0.1); %Let value5=(35 to 45 by 1);
		%Let ylabel1=BW, kg; %Let ylabel2=ECM/DMI, kg/kg; %Let ylabel3=Milk Fat Yield, kg; %Let ylabel4=Milk Fat Percentage, %; %Let ylabel5=Milk Yield, kg;
		%Let xlabel=Week of Trial;
		proc format;
		value $treatmentFmt
		"Bov" = "PRO-A"
		"BovC" = "PRO-B"
		"Control" = "CON" 
		;
		ods html style=listing; *In SAS 9.4 you need a style like LISTING to see different symbols;
		Proc SGPLOT data = Selected noautolegend;
			format treatment $treatmentFmt.;
			where resp_name="&&resp&i";
			styleattrs datasymbols=(CircleFilled triangleFilled SquareFilled) datalinepatterns=(ShortDash Solid MediumDash) datacontrastcolors=(darkblue darkred darkgreen); 
			Series x=Period y=estimate/ group=Treatment grouplc=Treatment groupmc=Treatment markers lineattrs=(thickness=1) name="Treatment";
			Scatter x=Period y=estimate/ group=Treatment name="a" YErrorUpper=Upper YErrorLower=Lower;*groupdisplay=cluster clusterwidth=0.6 datalabel=MSgroup datalabelattrs=(size=9pt color=gray);
			yaxis label="&&ylabel&i" Labelattrs=(Family=Arial Size=11 weight=Bold ) valueattrs=(size=10) values=&&value&i;
			xaxis discreteorder=data label="&xlabel." Labelattrs=(Family=Arial Size=11 weight=Bold ) valueattrs=(size=10) values=(1 to 12 by 1);
			keylegend "Treatment"/ title="Treatment:" location=inside position=topright across=1;
		run;
 	%end;
		Ods html CLOSE;
 %mend;
%TestSee;

What i however need is a grid of all 5 graphs with independent axis, all labelled as i did with the above graph. I did research on Proc template and it seems it could really do the job for me, but i still couldn't figure out how to make this work. I will really appreciate your help guys. Thank you.

Superscript in format for SGPLOT value labels

gabonzo — Tue, 01 Jun 2021 16:01:08 GMT

Hi,

I have defined a FORMAT that uses unicode superscripts:

	proc format;
		value agegroupfr 	1="1(*ESC*){super re} fois"
							2="2(*ESC*){super e} fois ou plus"
							3="18–24"
							4="25–34"
							5="35–44"
							6="45–54"
							7="55–64"
							8="65–74"
							9="75+"
							10="Toutes";
	run;

I want to use this format to be used to label values on the y-axes of a butterfly plot. Here is the portion of the code where I define the y-axes:

	 	yaxis reverse
		 	display=(noticks)
		 	valueattrs=(color=gray33 size=10pt family="Arial Unicode MS")
		 	labelposition=top
		 	colorbands=odd colorbandsattrs=(color=grayf3)
		 	label="&aglab."
		 	valuesformat=&agfmt.;
		y2axis reverse
			type=discrete
		 	display=(noticks)
		 	valueattrs=(color=gray33 size=10pt family="Arial Unicode MS")
		 	labelposition=top
		 	label="&aglab."
		 	valuesformat=&agfmt.;

The result is not what I expected:

I have tried several solutions proposed in this forum, but with no effect. What am I doing wrong?

Thanks.

SGPLOT change y axis label based on values of By variable.

James_Yu — Fri, 28 May 2021 06:38:18 GMT

Hi, I would like to use SGPLOT to create multi graphs using By-variable statement

- My code is like that:

Proc sgplot data=database;

by Symptom;

scatter x=VisitNum y=Mean / groupdisplay=cluster group= Treatment;

run;

- This is some of the graphs created:

My question is how could I change the label in the y axis to the name of variable (in this case: Symp1S Symp2S instead of Mean).

Thanks in advance

SGPLOT HEATMAP colormap labels/values size

marta25 — Wed, 26 May 2021 16:46:10 GMT

Hi,

I am using PROC SGPLOT with HEATMAP and I would like to change the size of the colormap values.

Here is my code and the resulting plot:

PROC SGPLOT data=  MARTA.mydata;
heatmap x=clocktime y=myvariable / colorresponse=aaa  colorstat=mean colormodel=(bigb red) xbinsize=0.25 ;
xaxis label='Clock time' values=(9 to 18 by 1) valueattrs=(size=14pt) labelattrs=(size=14pt);
yaxis label='My variable' valueattrs=(size=14pt) labelattrs=(size=14pt);
gradlegend / title='colormap' titleattrs=(size=14pt)  ;
run;

My problem seems to be simple but I can not find the solution. The values of the colormap are smaller compared to the other labels and values on the axes.

How to set such values to 14pt?

What is the explanation behind the limited graphical capabilities of SAS?

Malthe — Wed, 26 May 2021 08:29:56 GMT

I am a Ph.D. in medicine and have used SAS for a year and a half now. When starting my studies, I had to choose between SAS and R. The research-community in medicine in Denmark is split roughly 50/50 in SAS and R-users, with young, new users tending to choose R due to it being Open Source.

I chose SAS because most people at my department used it and it seemed like an easy way to get going. I have previous experience from Python, but wanted an easy start from the get-go. Now, SAS is very intuitive for data-management and most statistics (although the syntax for a statement like CONTRAST in Phreg is very unintuitive compared to the R-alternative).

I have grown to like SAS a lot, and therefore it both baffles and hurts me, that the graphical capabilities are so inferior to R. I hate having to export my datasets to be able to do what in R takes half an hour, but in SAS takes several hours and obscure code. I am sure that this is one of the reasons that SAS is loosing foothold in this battle.

How come a company with so big economic muscles and market-dominance doesn't prioritize this aspect of the program? Why can't you output vector-formats easily for editing in other programs? Why isn't there a wizard or some easy code, that recognizes all the typical objects you want to be able to manipulate, when having created a graph? I mean most of these functions are from a coding-perspective not that advanced compared to all the other mechanisms in SAS.

I hope that this side of the program will get some more attention - and a fresh new Googly/Apply-approach in the future!

Help with YAXIS Ranges= option

genemroz — Tue, 25 May 2021 14:59:04 GMT

Can someone please help me with the ranges option in the following code? It is not yielding the expected result of limiting the yaxis data to the range of 60 to 100%.

Thanks,

Gene

data data;
input LOHI $ k3_6 Group $;
datalines;
36/0 90.18 RoOPT
39/0 93.31 RoOPT
39/41 94.55 RoOPT
39/43 95.21 RoOPT
39/45 95.79 RoOPT
39/47 95.58 RoOPT
39/49 95.60 RoOPT
39/51 94.85 RoOPT
39/53 95.12 RoOPT
39/55 95.43 RoOPT
39/59 96.02 RoOPT
39/61 96.29 RoOPT
39/63 96.52 RoOPT
39/65 96.88 RoOPT
39/67 96.78 RoOPT
39/69 96.84 RoOPT
39/71 96.85 RoOPT
36/0 66.59 RoCOV
39/0 72.53 RoCOV
39/41 72.87 RoCOV
39/43 73.22 RoCOV
39/45 73.52 RoCOV
39/47 71.02 RoCOV
39/49 69.95 RoCOV
39/51 70.65 RoCOV
39/53 70.79 RoCOV
39/55 70.95 RoCOV
39/59 71.25 RoCOV
39/61 71.39 RoCOV
39/63 71.51 RoCOV
39/65 71.59 RoCOV
39/67 71.64 RoCOV
39/69 71.67 RoCOV
39/71 71.68 RoCOV
;
ods graphics / reset width=6.4in height=4.8in imagemap;

proc sgplot data=WORK.DATA;
title "Sensitivity of Optimization Results to Elevation Angle";
title2 "Region of Optimization (RoOPT): 440km by 440km";
title3 'Region of Coverage (RoCOV): All Gridpoints where k GE 1';

	vbar LOHI / response=k3_6 group=group groupdisplay=cluster;
	yaxis grid ranges=(60-100);
	xaxis label="LO/HI Elevation Angles";
	yaxis label="Gridpoints Covered by 3-6 Stations (%)";
run;

ods graphics / reset;

proc sgmap - removing white padding around edge of image

andersonh — Fri, 21 May 2021 23:05:35 GMT

I've been working with proc sgmap to generate some plots that are meant to be included in some printed material, and I'm trying to get rid of the white border around the edge of the plots. The images are being handed off to someone working in our communications department who's using a separate piece of software to format and typeset the printed materials, and I need to provide them with the images at a pre-defined resolution and DPI.

I can get the DPI of the overall image working just fine, and can turn off the titles/legends/borders without issue, but I've gotten completely stuck trying to make the whitespace go away. I got as far as experimenting with proc template and defining new styles, but I couldn't figure out how to make any of the pieces take effect. Currently I'm falling back on batch processing the images locally to trim the whitespace, but it would be much easier if I could just generate them without the whitespace in the first place, so I don't have to do some extra work to find how much is getting trimmed off (especially if I need to change the DPI, and the number of pixels being trimmed will change).

Unfortunately I can't share my data, but here's some code that illustrates the issue. It just draws a black box around Texas:

data work.ChoromapData;
	infile datalines delimiter=" ";
	input x y id $;
	datalines;
-94 25 Texas
-94 37 Texas
-106 37 Texas
-106 25 Texas
-94 25 Texas
;
run;

ods listing gpath="/home/my_user/" device=png;
ods graphics on / height=500 width=500 noborder;
proc sgmap mapdata=work.ChoromapData;
    openstreetmap; 
    choromap;
run;

The images I get out look like this. I've colored the whitespace around the edge a bright red, so it's easier to see against the forum background, but it's white in my program's output. Normally this would be the region between the map and the thin border that surrounds the whole image (but I've turned the border off with the ods graphics statement).

What I want is this (note the lack of whitespace around the edges):

I've been trying to get this sorted out for a few days now but have made zero progress. Anyone have any ideas/solutions that maybe I've just not come across yet?

Technical details: I'm working in Enterprise Guide 8.3, and SAS version 9.4M7 (full version number: 9.04.01M7P080520).

proc gmap : zoom part of map

Anita_n — Fri, 21 May 2021 16:37:41 GMT

I am tring to zoom this part of the map with data so that labels will fit in, also to enable proper visaulisation. I am using proc gproject but can't really get this working.

Proc gmap: displaying symbols for graph region

Anita_n — Fri, 21 May 2021 16:31:14 GMT

How can I display location symbols on the regions? some region names were not displayed, I think because of the tiny space. Is there any way to solve that?

How to add a symbol just for one group in sgplot for Hbox

Andrea_Peng — Fri, 21 May 2021 06:46:31 GMT

Hi,

I used sgplot to get hbox. I want to make one hbox (eg. the top bar) is different to show that is a specific case. Either to put a star aside the box or make the color different. Both I did not find how to make it. Hope some one can help.

Thanks,

Andrea

The graph looks like:

My code:

proc sgplot data=Caterpillar_graph noautolegend noborder;
hbox VAS_preop / category=surgeon_id nocaps extreme fillattrs=(color=CXE7AF19) WHISKERATTRS=(color=black) MEDIANATTRS=(color=CX2D3884 thickness=8) nomean;
refline 75 /axis=x lineattrs=(color=CX2D3884) LEGENDLABEL= "National Median" name='line';
refline '*'/axis=x;
yaxis label = "Individual Surgeon" display=(novalues noticks) labelattrs=(size=10pt FAMILY="Segoe UI Symbol") DISCRETEORDER=DATA;
xaxis label = "VAS" labelattrs=(size=10pt FAMILY="Segoe UI Symbol")  values = (0 to 100 by 10);
keylegend "line" ;
run;

proc gmap: color gradient shading

Anita_n — Fri, 21 May 2021 06:59:10 GMT

Dear all,

am tring to use color grading shades for the map I created but it's not showing any results. Am not sure if I have missing anything in my code. I will appreciate any help

proc template;
define style styles.grad1;
parent=styles.listing;
style threecolorramp /startcolor=cream endcolor=darkblue;
end;
run;

proc gmap data=mydata map=mapsgfk.germany all; 
id id;
choro cases / statistic=first discrete 
levels=4
coutline=gray anno=myanno
cdefault=CXFFFFFF
legend=legend1;
run;
quit;

Patient flow Sankey Diagram

Anita_n — Tue, 18 May 2021 09:06:57 GMT

Hi all,

I wanted to create this nice graph created by as shown on these links:

Outside-the-box: Circle link graph - Graphically Speaking (sas.com)

Outside-the-box: Directed circle link graphs - Graphically Speaking (sas.com)

I have a data of patients who visit a certain clinic, I will like to display this just like in the links above but vice-versa. This means in my case they are moving from their residents to the clinic.

I tried using the method @Jay54 used but I still need some help

data patientflow;
  input From $18. To $8. LinkCount ;
  datalines;
AAAAAAAAAAAAAAAAA   clinic	11
BBBBBBBBBBBBBBB	    clinic	5
DDDDDDDDDDDD	            clinic	30
EEEEEEEEEEEEEEEEE   clinic	2
FFFFFFFFFFFFF	            clinic	7
GGGGGGGGGGGGGG	   clinic	2
HHHHHHHHHHHHH	   clinic	11
IIIIIIIIIIIIIIII                           clinic	2
JJJJJJJJJJJJJJJJJ            clinic	25
KKKKKKKKKKKKKKK	   clinic	2
LLLLLLLLLLLLLLL	           clinic	 1
MMMMMMMMMMMM	   clinic	48
NNNNNNNNNNNNNN	   clinic	1
OOOOOO(clinic location)    clinic	500
PPPPPPPPPPPPP	           clinic	1
QQQQQQQQQQQQQQ	   clinic	3
RRRRRRRRRRRRRRR    clinic	1
SSSSSSSSSSSSSSSSS  clinic	1
TTTTTTTTTTTTTT	           clinic	1
UUUUUUUUUUUUUUU	   clinic	3
VVVVVVVVVVVVV	            clinic	7
WWWWWWWWWW	   clinic	7
XXXXXXXXXXXXXXX	   clinic	165
;
run;

please note that OOOOOO is the region where the clinic is located that is the reason why the number of cases is in this region very high. The names of the regions are quite long and the arrows should be coming from the different regions and ending at OOOOOO

I will appreciate any help or an alternative solution. Thanks

how to add tick marks with date variable

htst1003 — Mon, 17 May 2021 22:39:04 GMT

Hi, I am very beginner in SAS. I have date variable and 'daily COVID cases' variable.

When I make sgplot (xaxis=date and yaxis=# of cases), I only get a few tick marks on x axis. I want to add tick marks for each day. I have been searching to solve this problem for 4 days but none of article worked or I couldn't understand their explanations.....again I am very beginner, please explain basic details if you can.. Thank you so much.

This is my code.

What should be fixed?

proc sgplot data=after;
series x=date y=New_cases_fu;
xaxis interval=day label=' ' values=('08JUL2020'd to '16sep2020'd by 1) ;
yaxis label=" New cases";
RUN;

Remove bar chart tick labels but not ticks

RoddyJ — Fri, 14 May 2021 14:38:38 GMT

Hello,

I'm using SGPlot to plot a bar chart from a dataset that gets updated once a week with a new observation for a total number of events that week.

The x-axis is getting crowded so I'd like to display the text label for every other axis tick, but keep all the bars. I've tried using the values= option but that removes the bars as well as the tick labels. I've also tried using the FITPOLICY=thin option, which is closer to what I want as it removes tick labels but keeps the data. But that doesn't keep enough of the x axis ticks and removes the most recent tick label which is most important.

Can anyone help me out?

proc sgplot data=set5 noborder noautolegend;
yaxis label = 'Number' valueattrs=(size=10pt) LABELATTRS=(size=10pt) values=(0 to 100 by 5);
xaxis label = 'Week' valueattrs=(size=7pt) LABELATTRS=(size=6pt) values=(1 to &max_week by 2);
vbar week / datalabel datalabelfitpolicy=none datalabelattrs=(size=8) fillattrs=(color='CX0066CC');
format week donweek.;
title1 h=12pt j=left "Figure 2";
keylegend/ title=" ";
run;

3D Barchart

Manj — Thu, 13 May 2021 13:22:39 GMT

Hi,

This is the graph Im looking forward to create. Not really sure if this is one x axis and 2 yaxes Or 2 X axes and 1 Y axis. Can anyone show some light towards this? Any documentation shows a method to do this?

Thanks

Extracting Plots from procedure output

AlexeyS — Wed, 12 May 2021 17:42:24 GMT

Hi,

When i run the following code, i got in the output multiple plots.

proc ttest data=sashelp.class plots=all;
class sex; 
var age;  
title 'ttest of equality of means'; 
run; 
title;

I want to extract the plots and put some of them in the table two-by-two table. Can i do it?

I know that ODS OUTPUT destination enables to store any value that is produced by any SAS procedure(ods trace), but what about plots?

Thank you

Horizontal barchart with axis table, Labels won't align

andypandy_swe — Wed, 12 May 2021 06:43:23 GMT

I'm trying to create a horizontal barchart with a label above the y-axis values (left aligned) and also a label above the bars (left aligned).

In a simple GTL barchart I can get the "unit" label but couldn't get the "Some measurement" label above the bars (only at the bottom of the bars) so i figured I'd try and put the y-axis values and the "unit" label in an axistable. So I created a layout lattice with two layout overlays in it, one for the axistable (showing values and label) and one for the barchart (showing only the label). However as seen above the two labels won't align. The label for the axistable is higher than the label for the barchart.

So my question is:

Is it possible to get the x-axis label on top of the bars aligned with the y-axis label (and skip the axistable) or

keep the axistable and somehow align not only the axistable/barchart values but also the labels?

The template needs to accomodate diffrent number of bars (and thereby different graph heights) so using drawtext isn't really an option (as far as I can see)..

Fun With SAS ODS Graphics: MOM-WOW Mother's Day Card

tc — Sat, 08 May 2021 04:06:16 GMT

A quick SAS ODS Graphics text/scatter plot take on those MOM-WOW cards. Happy Mother's Day, all!

* Fun With SAS ODS Graphics: MOM-WOW Mother's Day Card
  Mother icon by DewDrops on @NounProject thenounproject.com/term/mother/851501
  Mother's Day color palette from codeofcolors.com/mothers-day-colors.html;

data mom;                                                                                           * Generate points and text for card;
retain x y 0 text 'MOM' x2 0 y2 .87;                                                                * Text to display and x/y points;
                                                                                                    * x,y="MOM" x2,y2=Mother icon;
ods graphics / reset height=5in width=5in antialias;                                                * 5"x5" image (good enough for Twitter!);
proc sgplot data=mom aspect=1 noborder noautolegend nowall;                                         * Generate MOM-WOW Mother's Day card (text/scatter plots);
styleattrs backcolor=cxFEE8DA;                                                                      * Specify background fill color;
symbolimage name=mom image='/folders/myfolders/motherchild.png';                                    * Mother icon; 
scatter x=x2 y=y2 / markerattrs=(symbol=mom size=34pt);                                             * Scatter plot w/mother icon marer; 
text x=x y=y text=text / textattrs=(color=cxE25B6A size=108pt weight=bold) position=top contributeoffsets=none; * MOM; 
refline 0 / axis=y lineattrs=(color=cx651B0E thickness=3pt);                                        * Line separator at 0;
text x=x y=y text=text / textattrs=(color=cxFFAEAA size=108pt weight=bold) rotate=180 position=top contributeoffsets=none; * WOW (MOM "reflection"); 
inset " HAPPY MOTHER'S DAY " / position=bottom textattrs=(color=cx651B0E size=22pt weight=bold);    * HAPPY MOTHER'S DAY;
xaxis display=none min=-1 max=1 offsetmin=0 offsetmax=0;                                            * Specify x/y axes ranges, hide axes;
yaxis display=none min=-1 max=1 offsetmin=0 offsetmax=0;

Bar plot, To present symbol (Star, Plus sign) on X-axis For specific bars (separate variable)

Sagarg — Fri, 07 May 2021 09:33:49 GMT

I am having simple Bar plot. See below code. In SCORE dataset for name "LMN" doesn't have any marks. I want to present star(*) at the x-axis below tickvalues (I can hide tickvalues).

Any thoughts How to present Star(*) at axis, indicating it is different from dataset row having no values.

data score;

input name $1-3 marks 5.;

datalines;

abc 50

bcd 70

LMN

kgh 90

uyt 10

;

run;

proc template;

define statgraph temp;

begingraph;

layout overlay;

barchart x=name y=score;

endlayout;

endgraph;

end;

run;

proc sgrender data=score template=temp;

run;

"Ask the Expert” webinar entitled How Do I Fix Common Problems in Output From ODS Graphics?

MarciaS — Wed, 05 May 2021 14:27:44 GMT

Hi SAS Graphics Programming Community,

I’m presenting a live “Ask the Expert” webinar on May 11^th, 10 – 11 AM ET entitled How Do I Fix Common Problems in Output From ODS Graphics? Join this webinar to learn the common pitfalls when using ODS Graphics procedures and how to avoid those issues in the future. I’ll also discuss ways to improve the overall appearance of a graph by using attribute maps or adding special symbols. I’ll host a live Q&A session at the end of the webinar.

You will learn:

How to get your desired colors and symbols in the graph output.
How to place annotation on a graph.
How to improve the overall look of your graph output. the advanced features of SAS survey procs.

Reserve your spot for the webinar.

Want more tips? Be sure to subscribe to the Ask the Expert board to receive follow up Q&A, slides and recordings from this and other SAS Ask the Expert webinars.

Can't join the live event? You can view this and other Ask the Expert sessions on-demand here.

Thanks,

Marcia

SGPlot annotation on a grouped VBAR plot

Chris_Gibby — Fri, 30 Apr 2021 10:35:02 GMT

Hello,

I'm having some trouble displaying annotated data on a grouped VBAR. Example code below.

I want to display an 'x' at the bottom of each Bar, i.e. the position of the 'x' is tied to the data value of each bar.

However when I run I get the following Warning and the 'x' is not displayed on the output:

WARNING: DrawText statement has missing/invalid value for position (X or Y). Draw statement discarded.

data data1;
 length genx genid $ 10;
  genx = 'GRP1';
  genid = 'HRD1'; d1=15; output;
  genid = 'HRD2'; d1=23; output;
  genid = 'HRD3'; d1=5; output;
  genid = 'HRD4'; d1=35; output;
  genx = 'GRP2';
  genid = 'NONHRD1'; d1=43; output;
  genid = 'NONHRD2'; d1=28; output;
  genid = 'NONHRD3'; d1=26; output;
  genid = 'NONHRD4'; d1=9; output;
run;

data anno_ds;
 set data1;
   function = 'text';
   x1space = 'datavalue';
   y1space = 'layoutpercent';
   textsize = 6;
   textweight = 'normal';
   width = 20;
   widthunit = 'percent';
   anchor = 'center';
   xc1 = genid;
   y1 = -1;
   label = 'x';
run;

  ods graphics on / reset height=6.8in width=8.5in imagename='test' noborder;
  
  proc sgplot data=data1 noautolegend sganno=anno_ds  pad=(bottom=10pct);
   vbar genx / group=genid response=d1 groupdisplay=cluster 
              dataskin=pressed ;
              
   xaxis label='Grouped data' display=(noticks) valueattrs=(size=7 family=Arial weight=bold) labelattrs=(size=7 family=Arial weight=bold);
   yaxis label='Score (%)' values=(0 to 100 by 10) valueattrs=(size=7 family=Arial weight=bold) labelattrs=(size=7 family=Arial weight=bold);
  run;

Thank you for looking,

Chris.

How to change the color in BARCHART

Coooooo_Lee — Tue, 04 May 2021 10:22:42 GMT

Hi there,

I just want to change the bar color in my GTL temple.

No group in my data, just want to show all bars in a same color, like the RED I specified here....
However, the bar color always black...

Anyone can help? Thanks!

proc template;
 define statgraph swimmer;
  begingraph / drawspace=datavalue;


  layout lattice / rows=1 columns=3 columnweights=(0.1 0.75 0.15) columngutter=0;

    layout overlay;
      entry textattrs=(color=white size=5) "blank";
    endlayout;

    layout overlay / walldisplay=none

        xaxisopts = (label="xxxxx"
                     labelattrs=(size=7)
                     tickvalueattrs=(size=7)
                     linearopts = (viewmin=0 viewmax=140
                                   tickvaluesequence = (start=0 end=140 increment=1)))


      yaxisopts = (label="xxxxxx"  
                     labelattrs=(size=7)
                     display=(label)
                     linearopts = (viewmin=1 viewmax=215
                                   tickvaluesequence = (start=1 end=215 increment=1)));


      barchart  x=obs y=dur /  orient=horizontal barwidth=0.1 fillattrs=(color=yellow);

        endlayout;

   
      layout overlay;
          entry textattrs=(color=white size=5) "blank";
        endlayout;   
    endlayout;
  endgraph;
  end;
run;

sgplot - change colors for each group in histogram and add statistics

AlexeyS — Fri, 30 Apr 2021 07:28:55 GMT

Hi,

I created a histogram and added a density plots for each group.

1. I want to change the colors of each group(the colors of histogram and density have to be the same for each group)

2. I want to add some basic statistics such mean, number of observation for each group. I tried to do it using inset statement, but my problem, how can i use inset statement for each group, i.e the first row statistics will be for the first group, the second row, will be for the second group

ods graphics / noborder;

proc sgplot data=test;
  histogram amount   / group=id transparency=0.5;       /* SAS 9.4m2 */
  density   amount   / type=normal group=id;      /* overlay density estimates */
  legenditem type=marker name="1"  / label="Try1" markerattrs=(color=red symbol=circle) labelattrs=( Family="Arial" Size=10 Style=italic Weight=bold);
  legenditem type=marker name="2"  / label="Try2"  markerattrs=(color=blue symbol=circle) labelattrs=(Family="Arial" Size=10 Style=italic Weight=bold);
  keylegend "1" "2"/ location=inside position=topright across=1 noborder AUTOITEMSIZE;
  inset "Average = &avg" "NObs = &nobs" / border title="Try1" TITLEATTRS=(Color=Green Family=Arial Size=8 Style=Italic Weight=Bold) position=topleft TEXTATTRS = (SIZE=6 COLOR=red Family=Arial Size=8 Style=Italic);
run;

Thank you

Creating boxplots

Anita_n — Thu, 29 Apr 2021 09:20:52 GMT

Hi all,

I will like to create a boxplot for the following 3 variables (var1, var2, var3) as in the table attached. Is there any way to achieve that? they should be displayed side by side in one diagramm

proc sgplot label y axis with another character variable

zimcom — Tue, 27 Apr 2021 18:06:11 GMT

Greetings!

proc sgplot data=example noautolegend;
styleattrs datasymbols=(circlefilled );
scatter x=b y=a ;
keylegend / location=outside position=NW across=1;
run;

For the above code, is it possible to label the y axis using variable "CAT" to make the y-axis more meaningful?

CAT	a	b
CAT1	1	2
CAT2	2	3
CAT3	3	5
CAT4	4	1
CAT5	5	6
CAT6	6	4

Change barwidth sgplot histogram

Anita_n — Mon, 26 Apr 2021 10:50:32 GMT

Hi all,

is there any way to change the barwidth

in sgplot histogram statement?

Categorical Data on X Axis

RaquelAR — Fri, 23 Apr 2021 17:30:49 GMT

Hi,

I am analyzing data of different strains with different genotypes and I am trying to format the x axis a certain way. My data is categorical and I would prefer not to change it numbers. I added a xaxis values statement and now my graph looks crazy. Any suggestions? I attached photos of before I used the xaxis value statement and after. This is my code:

*Creating figure plotting estiamted probs & CI**;

proc logistic data=estimated2 alpha=0.2 noprint;
class mosquitostrain (ref= 'NO') genotype / param=ref;
model Alive(Event='1')= mosquitostrain genotype;
/* 1. Use a procedure or DATA step to write Pred, Lower, and Upper limits */
output out=LogiOut predicted=estprob2 l=lower95 u=upper95;
run;
/* 2. Be sure to SORT! */
proc sort data=LogiOut;
by mosquitostrain genotype;
run;
/* 3. Use a BAND statement. If more that one band, use transparency */
title "Predicted Probabilities with 95% Confidence Limits";
title2 "Four Strains";
proc sgplot data=LogiOut;
band x=genotype lower=Lower95 upper=upper95 / group=mosquitostrain transparency=0.75;
series x=genotype y=estprob2 / group=mosquitostrain curvelabel;
xaxis grid;
xaxis values= ('WT,WT' 'WT,HET' 'HET,WT' 'HET,HET' 'MUT,WT' 'WT,MUT' 'HET,MUT' 'MUT,HET' 'MUT,MUT');
yaxis grid label="Predicted Probability of Survival" max=1;
keylegend "L" / location=inside position=NW title="Genotype" across=1 opaque;
run;

Multiple graphs in one panel

wj2 — Fri, 23 Apr 2021 15:38:14 GMT

Hello,

I am trying to combine 4 separate survival plots into a single, panel style graph. Please see below for my code for each of the individual survival plots. Any help with this would be greatly appreciated!

/***Graph #1***/
ods graphics on;
ods output survivalplot=OverallSample;
proc lifetest data=have plots=survival; 
time Days*censor(0);
run;
ods graphics off;

proc sgplot data=OverallSample;  
step y=survival x=time;
run;

/***Graph #2***/
ods graphics on;
ods output survivalplot=VariableA;
proc lifetest data=have plots=survival; 
strata varA;
time Days*censor(0);
run;
ods graphics off;

proc sgplot data=VariableA;  
step y=survival x=time/ group=stratum;
run;

/***Graph #3***/
ods graphics on;
ods output survivalplot=VariableB;
proc lifetest data=have plots=survival; 
strata varB;
time Days*censor(0);
run;
ods graphics off;

proc sgplot data=VariableB;  
step y=survival x=time/ group=stratum;
run;

/***Graph #4***/
ods graphics on;
ods output survivalplot=VariableC;
proc lifetest data=have plots=survival; 
strata varC;
time Days*censor(0);
run;
ods graphics off;

proc sgplot data=VariableC;  
step y=survival x=time/ group=stratum;
run;

Gaussian curves in SAS

Anita_n — Thu, 22 Apr 2021 06:14:36 GMT

Hi all,

any idea how to create Guassian curves in sas. Any code sample?

I guess it's easy but it's new for me. My client wants that for a set of diagnosis ages in a DB (Number of cases, Mean, Standard Deviation)

Thanks

Editing labels on Forest Plot in SGPlot and reordering variables

joesmama — Wed, 21 Apr 2021 20:45:30 GMT

Hello,

I haven't used SGPLOT aside from the very basics. I have a Forest plot for which I would like to edit the labels. I would also like to reorder the variables so they match the order in the odds ratio table. In my example below, I would like to remove the variables names "Race_new" and "Sex" but keep "Age". I would also like to put "Age" in top followed by "Sex" and then "Race". I've been scouring SAS resources but I can't seem to find anything to address these questions specifically. My code can be found below. Any help would be greatly appreciated.

proc sgplot data=orci_a;
scatter x=oddsratioest y=effect / xerrorlower=lowercl
xerrorupper=uppercl
markerattrs=or
(symbol=DiamondFilled size=8);
refline 1 / axis=x;

xaxis label="OR and 95% CI " min=0;
yaxis label="Covariates";

run;

sgplot scatter and series (supress one Legend)

Anita_n — Tue, 20 Apr 2021 20:22:12 GMT

Hi all,

I used sgplot to create a scatter and series plot. I get two legends as result. How can are supress one since they are both the same.

Adding Standard errors on graphs

Anita_n — Tue, 20 Apr 2021 13:10:16 GMT

Dear all,

I have a graph showing the mean ages in different diagnosis years. I wish to add standard errors to this plot. Is it right to do that with proc means like below and then deduct or add these values from/to the maen ages to get the upper and lower limits. Or do I say lower= mean_age- 2*stderr ,

and upper= mean_age+2*stderr. Why is this standard error important?

proc means data=mydata noprint;
var age;
class year sex;
output out = myoutdat mean=mean_age stderr=stderr ;
run;

how to label bins on x axis in histogram

fengyuwuzu — Tue, 20 Apr 2021 00:29:57 GMT

I am trying to plot histograms with proc sgpanel, but the default x-axis labeling is 0, 50, 100, 150.

I want to show 0, 10, 20, 30, 40, ... or 0, 20, 40, 60, 80,...

I tried with rdi=(0 to 200 by 10) but it does not work.

Any suggestions? Below is my code:

proc sgpanel data=exsum0 (where = (armcdn in (1, 2, 3)));
  panelby armcdn / rows=3 layout=rowlattice;
  histogram rdi /BINSTART=0 BINWIDTH = 10 rdi=(0 to 200 by 10);
run;

sgplot: using data response and datalabel

Anita_n — Mon, 19 Apr 2021 11:01:12 GMT

Hi all,

Is it possible to specify a response and a different datalabel in sgplot. I mean I wish to display the relative (percent) values on the yaxis but the absolute on that bars. I used the following code but the pecentages are display in both cases. Is it possible to do that anyway?

proc sgplot data =mergedat dattrmap=myattrimap pad=(bottom=5%) noautolegend noborder;
xaxis grid type=discrete discreteorder=data;
vbar year response=percentvar group=grup attrid=grup  datalabel = absolutevar datalabelattrs=(weight=bold) seglabel seglabelattrs=(size=0.25) 
seglabelfitpolicy=none barwidth=0.9 missing;
xaxis  display= (nolabel) valueattrs=(family="verdana" size=8pt) fitpolicy=rotate ValuesRotate=Vertical;
yaxis values=(0 to 60  by 10) display= (nolabel) valueattrs=(family="verdana" size=8pt);

run;

sgrender: alignemt of tick marks

Claus_Stenberg — Mon, 19 Apr 2021 11:27:50 GMT

Doing a adverse event plot were I have to put some stat numbers together with a scatter plot (see below). Now I want the stat numbers to be right aligned with the tickmark at the x axis (see red inserted line).

proc template;
   define statgraph demo;
      begingraph ;
            layout overlay;
                      scatterplot y=name x=w / markercharacter =weightc   MARKERCHARACTERPOSITION=left    xaxis=x2;
                       scatterplot y=name x=h / markercharacter =height   MARKERCHARACTERPOSITION=left    xaxis=x2;
             endlayout; 
      endgraph;
   end;
run;
 
data class2;
 set sashelp.class;
 weightc=put(weight,4.1);
 heightc=put(height,4.1);
 w="kg";
 h="cm";
run;

proc sgrender data=class2 template=demo;
run;

Have created a little sandbox example using the good old sashelp.class :-). Any ideas on how to solve this?

SAS KEEPS RUNNING after ERROR 79-322: Expecting a -.

Giampaolo — Sun, 18 Apr 2021 10:44:08 GMT

Dear SAS users,

After I run proc template sas keeps running and I receive the ERROR 79-322: Expecting a -.

Any idea, how can I stop SAS from processing the command and fix the issue?

Thank you!

Proc template;
define statgraph BrokenAxis1;
begingraph /datacolors=(CXD05B5B CX6F7EB3) datacontrastcolors=(red blue);
layout lattice / Columns=2 Rows = 2 columnweights= (.5 .5)
Rows=2 rowweights= (.47 .53) rowgutter = 5
rowdatarange=UNION;
/* Lattice header definitions */
column2headers ;
entry textattrs=(size=9 ) "ELECTIVE"/*/pad=(left=15 px)*/;
entry textattrs=(size=9 ) "EMERGENT"/*/pad=(left=15 px)*/;
endcolumn2headers;
/*define row 1 column 1*/
layout overlay / xaxisopts= (display= (line))
yaxisopts=(display=(line ticks tickvalues label)
label="Total Cost (X$10,000)"
griddisplay=on
linearopts=(includeranges=(min-170000 380000-384000)
tickvalueformat= bigvalues.
tickvaluelist=(20000 40000 60000 80000 100000 120000 140000 160000 380000)));

boxplot x=&GRP Y=eval(ifn(&strata=0,Total_cost, .))/
grouporder=ASCENDING
INCLUDEMISSINGGROUP=FALSE
group=&XVAR
groupdisplay=cluster
meanattrs = (SYMBOL = DIAMOND weight=bold )
medianattrs= (thickness = 2)
name ="L" ;
endlayout;
/*define row 1 column 2*/

layout overlay /xaxisopts= (display= (line))
yaxisopts=(display=(line ticks tickvalues)
griddisplay=on
linearopts=(includeranges=(min-170000 380000-max)
tickvalueformat= bigvalues.
tickvaluelist=(20000 40000 60000 80000 100000 120000 140000 160000 380000)));

boxplot x=&GRP Y=eval(ifn(&strata=1,Total_Cost, .))/
grouporder=ASCENDING
INCLUDEMISSINGGROUP=FALSE
group=&XVAR
groupdisplay=cluster
meanattrs = (SYMBOL = DIAMOND weight=bold )
medianattrs= (thickness = 2)
name ="L" ;
endlayout;
/*define row 2 column 2*/
layout overlay /xaxisopts=(display=(ticks tickvalues line))
yaxisopts=(display=(ticks tickvalues line label)label="Postop LOS (days)"
LABELFITPOLICY=SPLITALWAYS LABELSPLITCHAR='#' LABELSPLITCHARDROP=TRUE
griddisplay=on
linearopts=(tickvalueFormat= best12.
includeranges=(min-42 115-118)
tickvaluelist=(5 10 15 20 25 30 35 40 115)));

boxplot x= &GRP Y=eval(ifn(&strata=0,PostopLOS, .)) /grouporder=ASCENDING
INCLUDEMISSINGGROUP=FALSE
group=&XVAR
groupdisplay=cluster

meanattrs = (SYMBOL = DIAMOND weight=bold )
medianattrs= (thickness = 2)name ="L" ;

endlayout;

/*define row 2 column 2*/

layout overlay /xaxisopts=(display=(ticks tickvalues line))
yaxisopts=(display=(ticks tickvalues line)
griddisplay=on
linearopts=(tickvalueFormat= best12.
includeranges=(min-42 115-max)
tickvaluelist=(5 10 15 20 25 30 35 40 115)));

boxplot x= &GRP Y=eval(ifn(&strata=1,PostopLOS, .)) /
grouporder=ASCENDING
INCLUDEMISSINGGROUP=FALSE
group=&XVAR
groupdisplay=cluster

meanattrs = (SYMBOL = DIAMOND weight=bold )
medianattrs= (thickness = 2)name ="L" ;

endlayout;
sidebar / align = bottom; discretelegend "L" / halign=center across= 2 border=true VALUEATTRS= (size= 9);
endsidebar;
endlayout;
endgraph;
end;
run;
quit;

Relative Risk plot

Manj — Wed, 14 Apr 2021 18:57:11 GMT

The Expected plot is below:

I have the following code

ods path show;
ods path(prepend) work.templat(update);

ods pdf file="&g_outfile..PDF" nogtitle /*nogfootnote*/ bookmarkgen=no;
ods graphics on / reset=all width=1000px height=500px border=off ;

proc template;
define statgraph AEbyRelativeRisk;
dynamic _thk _grid;
begingraph;
layout lattice / columns=2 rowdatarange=union columngutter=1;

/*--Row block to get common external row axes--*/
rowaxes;
rowaxis / griddisplay=_grid display=(tickvalues) tickvalueattrs=(size=5);
endrowaxes;

/*--Column headers with filled background--*/
column2headers;
layout overlay;
entry " ";
endlayout;
layout overlay ;
entry "<--- Tmt1 Better Tmt2 Better --->";
endlayout;
endcolumn2headers;

/*--Left side cell with proportional values--*/
layout overlay / xaxisopts=(label="Percent" display=(ticks tickvalues label) labelattrs=(size=8) tickvalueattrs=(size=7));
referenceline y=teaex / lineattrs=(thickness=_thk) datatransparency=0.9;
scatterplot y=ae x=tmt1Pct / markerattrs=graphdata2(symbol=circlefilled size=10)
name='a' legendlabel="tmt1 (N=&NA)";
scatterplot y=ae x=tmt2Pct / markerattrs=graphdata1(symbol=trianglefilled size=10)
name='b' legendlabel="Tmt2 (N=&NB) ";
discretelegend 'a' 'b' / border=true halign=left valign=bottom valueattrs=(size=8) displayclipped=true;
endlayout;

/*--Right side cell with Relative Risk values--*/
layout overlay / xaxisopts=(label="Relative Risk with 0.95 CI (*ESC*){unicode '000a'x} tmt1 vs tmt2"
labelattrs=(size=8) tickvalueattrs=(size=7));
referenceline y=teaex / lineattrs=(thickness=_thk) datatransparency=0.9;
scatterplot y=ae x=RelRisk / xerrorlower=LowerCL xerrorupper=UpperCL
markerattrs=(symbol=circlefilled size=10);
referenceline x=1 / lineattrs=graphdatadefault (pattern=solid color=grey);
endlayout;

endlayout;
endgraph;
end;
run;

ods graphics / reset width=6in height=4in imagename='AEbyRelativeRisk';

proc sgrender data=ip template=AEbyRelativeRisk;
dynamic _thk='0' _grid='on';
run;

ods pdf close;

The figure I get is as below:

The issues are:

1) how to split each of the y axis tick values into two or more lines.

2) the figure looks very pale ,how to improve the texts and lines to be a bit more bold or darker so that it will be readable.

Creating text in sas and exporting it to the XL

JaneS — Mon, 12 Apr 2021 18:53:55 GMT

I created a couple of data sets in sas that export to excel. I need to create a description tab that describe my data sets. I was tying to find out a code sample for SAS, where I could put some text in that code and it would create a tab in XL file with that text in it.

Thank you!

How to prevent slice labels from overlapping pie chart (gchart)

Anita_n — Sat, 10 Apr 2021 18:02:43 GMT

Hi all,

can anyone help me out with this issue. The labels of my created pie chart overlaps. Is there any way to avoid that in a gchart?

Gchart clustered bar chart

JaneS — Fri, 09 Apr 2021 18:14:39 GMT

I have a question how to create a clustered graph is SAS.

The data is the following:

Year	County	Statewide
2015	20.00%	50.00%
2016	40.00%	32.00%
2017	30.00%	12.00%
2018	70.00%	0.00%
2019	90.00%	80.00%
2020	100.00%	90.00%
2021	100.00%	100.00%

And the graph I am trying to create in SAS is the following:

Could you please help with the code for it? I was thinking to use proc gchart but then it looks like there should be two vbar arguments and I am not sure how to combine them.

Thank you!

Be "that one guy"!

RobertAllison_SAS — Thu, 08 Apr 2021 15:44:26 GMT

Too good not to share... 🙂

proc sgpie

Anita_n — Thu, 08 Apr 2021 11:14:52 GMT

Hi all,

I saw this pie chart (using proc sgpie) which I really like. But it seems not to work in sas base(32 bits). Is there any way do create similar charts in sas base (32bits)? using gcharts or alternatives

Creating Custom legend

ewillisva — Wed, 07 Apr 2021 18:13:17 GMT

I need a way to produce only a legend which can be used on a panel in a template I have created. It's nothing fancy, just a title, 4 colored squares (with Hex or color name specified) with some text next to each. I've seen some SAS programming use gplot to produce only a legend, but that isn't working for my purposes. I came across Lex Jansen's page on it, but wasn't able to use that either. Anyone have some canned SAS code that will do this? I just need a way to call it in PROC GREPLAY so that I can add it to a template.

Values less than or equal to percent-of-total shouldn't be summarised as Other gchart(pie)

Anita_n — Sun, 04 Apr 2021 15:17:02 GMT

Hello all,

is there any option to add to this code, so that small values are not summarised as other in my pie chart?

legend1 label=none
        shape=bar(4,1.5)
        position=(top left)
        offset=(5,)
        across=1
        mode=share;

proc gchart data=mydat;
    pie var1 / sumvar=cases
	            noheading
                     coutline=black
				value= inside
				slice=arrow
				percent=outside
				legend=legend1;
			
run;
quit;

Fun With SAS ODS Graphics: Ellipses Easter Egg

tc — Sat, 03 Apr 2021 22:53:48 GMT

Ellipses Easter Egg

Drawing an egg could involve a lot of math...or you can cut a few corners and join the top of an ellipse with the bottom of a circle. Happy Easter, all!

* Fun With SAS ODS Graphics: Ellipses Easter Egg
  Bunny icon courtesy of Marie-Louise Merser, DK thenounproject.com/search/?q=easter+bunny&i=960660;

data egg;                                                            * Generate points and text for Easter Egg;
retain id 0 xt1 0 yt1 .195 text1 'HAPPY' xt2 0 yt2 0 text2 'Easter' xb 0 yb -.375; * Text to display and x/y points;
x=-1; y=0; output; x=1; output; y=-1; output; x=-1; output;          * Rectangle points (used to mask bottom half of top ellipse);

ods graphics / reset antialias width=5in height=5in;
proc sgplot data=egg aspect=1 noborder noautolegend nowall;          * Generate Easter Egg (half-ellipse + masking polygon + circle + text + scatter plots);
symbolimage name=bunny image='/folders/myfolders/easterbunny.png';   * Rabbit icon; 
ellipseparm semimajor=.75 semiminor=1 / nooutline fill fillattrs=(color=purple);* Top ellipse of Easter Egg;
polygon x=x y=y id=id / nooutline fill fillattrs=(color=white);      * Mask off bottom half of top ellipse;
ellipseparm semimajor=.75 semiminor=.75  / nooutline fill fillattrs=(color=purple);* Bottom ellipse of Easter Egg (circle);
text x=xt1 y=yt1 text=text1 / textattrs=(color=white size=18pt);     * Happy; 
text x=xt2 y=yt2 text=text2 / textattrs=(color=white size=63pt);     * Easter; 
scatter x=xb y=yb / markerattrs=(symbol=bunny size=48pt);            * Easter Bunny image; 
xaxis display=none min=-.925 max=.925 offsetmax=.001 offsetmin=.001; * Limit x/y axis bounds;
yaxis display=none min=-.80 max=1.05 offsetmax=.001 offsetmin=.001;

Ellipses Easter Egg "Wireframe"

Request for help with %sganno and %sgrectangle

genemroz — Sat, 03 Apr 2021 16:36:09 GMT

I'm trying to use %SGANNO and %SGRECTANGLE to overlay a square centered on the intersection of the reference lines of the following graph. Seems simple enough, but no square appears. Must be something simple I'm overlooking and would greatly appreciate your insights. The code I'm using appears below the graph.

Thanks,

Gene

/* create annotation data set */
%sganno;
 data anno;
 %SGRECTANGLE(drawspace="datavalue",x1=0,y1=0,width=200,height=200);
 output;
 run;
/*Create HeatMap*/
ods graphics / reset width=6.4in height=6.4in imagemap;
proc sgplot data=work.optcamer_k_coverage sganno=anno aspect=1;
where tz=80;
title height=12pt "OPTModel Camera Coverage @100km for 200x200km grid at 80km";
	heatmap x=Tx y=Ty / name='HeatMap' discretex discretey 
		colorresponse=k_coverage colormodel=(cxfeebe2 cxfbb4b9 cxf768a1 cxc51b8a cx7a0177)
		transparency=0;
	gradlegend /title='Number of Cameras';
	xaxis label='KM' valuesrotate=vertical grid values=(-300 to 300 by 10);
	yaxis label='KM' grid values=(-300 to 300 by 10);
	refline "0" / axis=x lineattrs=(thickness=2 color=blue) label 
		labelattrs=(color=blue);
	refline "0" / axis=y lineattrs=(thickness=2 color=blue) label 
		labelattrs=(color=blue);
run;

How to fill the peaks and valleys of a series plot?

Quentin — Sat, 03 Apr 2021 13:27:08 GMT

There was an interesting question on LinkedIn, where the desire was to fill the peaks of a series plot with one color, and the valleys as another color.

https://www.linkedin.com/feed/update/urn:li:activity:6783694808239128576

I'm wondering how folks would approach this.

Graphing guru Kriss Harriss had the good idea of using bandplot, and posted the approach in his excellent Daily Coding series:

https://www.linkedin.com/posts/krissharris_dailycoding-day91-sgf-activity-6783513455380041728-vO9n

When I played with this, I started with a sine wave and the goal of coloring the peaks and valleys different colors, without specifying the x-values to define each peak. I got close to what I wanted, but unfortunately there is no /BREAK option for the band statement, so SGPLOT connects the peaks.

My code:

data foo;
  do x = 0 to 3.14*4 by 0.01;
    y = sin(x);   
    output;
  end;
run;

data foo2 ; 
  set foo ;
  if y > .5 then do ;
    group=1 ;
    limit=.5 ;
  end ;
  else if y < -.75 then do ;
    group= 2 ;
    limit= -.75 ;
  end ;
run ;

proc sgplot data=foo2 ;
  series x=x y=y ;
  band x = x lower = limit upper = y   / group=group nomissinggroup;
run ;

My output, within the unwanted connecting bands highlighted:

I thought this was an interesting question, so was curious about other approaches.

Since the OP has reference lines for each band, I think my approach might actually 'work' because you could make the reference lines cover the unwanted band lines.

datalabelfitpolicy= rotate and xaxis fitpolicy= rotate

Anita_n — Fri, 02 Apr 2021 11:04:44 GMT

Dear all,

I have a problem using sgplot fitpolicy. Nothing seems to work.

I used datalabel and want the values to rotate so that they are horizontally displayed but it's not working (see myfilename1.pdf). Please not that this only seen when the results is outputted as pdf.

Is there any way to suppress 0 values?

My second problem is the xaxis fitpolicy =rotate which also doesn't work. Attached are all my test data.

I will appreciate any help.

Here is my code

ods graphics/ imagefmt = pdf width= 3in height= 2in noborder push;
goptions device=pdfc gunit=pct rotate=landscape ftext= "verdana" ;
ods listing gpath="mypath/myfolder" ;
ods graphics on / reset=index imagename="myfilename";

footnote justify=center "Cases per Agegroup per year" ;
proc sgplot data=test dattrmap=myattrmap pad=(bottom=5%) noautolegend noborder; 
format sex $sex.;
xaxis grid type=discrete discreteorder=data;
vbar AGEGROUP /response=CASES_PER_AGEGROUP group=sex datalabel datalabelfitpolicy= rotate
datalabelattrs=(size=0.25 color=black family= "verdana") groupdisplay=cluster attrid=sex
barwidth=0.9 missing ; 
xaxis  display= (nolabel) valueattrs=(family="verdana" size=8pt) fitpolicy=rotate ValuesRotate=Vertical;
yaxis values=(0 to 250 by 5) display=(nolabel) valueattrs=(family="verdana" size=8pt) grid offsetmin=0.0 offsetmax=0.0;
keylegend / location=inside position=topleft across=0 title="" valueattrs=(size=10 family="verdana") outerpad=(top=0.2cm) noborder;
run;
footnote " ";

How do I add specific sections in figures ?

_sba_ — Thu, 01 Apr 2021 14:51:50 GMT

Hi, any idea how can I do to program the below in SAS ? I mean, the sections in yellow and red that are based on 'y' values of 54 mg/dL, 70 mg/dL, 180 mg/dL and 250 mg/dL ?

Any advice is welcome!

Thanks!

sgplot barchart with n on y-axis and % floating over the bars

telligent — Wed, 31 Mar 2021 20:13:10 GMT

How do I make a bar chart with N on the y-axis, a category on the x-axis and percent as a label on the bars without creating an output data set with proc freq or something similiar. Just using raw data. . .

Thank you.

How do I have forecast and actual times series on a same graph

thuphammm — Mon, 29 Mar 2021 08:20:17 GMT

Hi, I've just done estimating and forecasting using ARIMA, below is my code. I want to do the Hold-out method where I have my test set and my training set separated, then compare the test set with what I forecasted, I want to have both of them on the same graph now. Please help!

ods html close;
ods html;

data tesla;
infile "/folders/myfolders/TSLA (1).csv" delimiter=',' firstobs=2;
input date :anydtdte. open high low close adj_close vol;
format date MONYY7.;
informat date MONYY7.;
proc print data=tesla;

ods graphics on;
proc sgplot data=tesla;
series x=date y=open;
Xaxis values=('01Jun2010'd to '01Jun2019'd by year) valuesformat=year4.;
title 'TSLA Opening Price';
run;

proc arima data=tesla;
identify var=open(1);
estimate p=3 q=2 method=ml noconstant; 
forecast lead=12 out=teslaforecast; 
proc print data=teslaforecast;
run;

Proc sgplot in SAS university edition

uzma03505621 — Fri, 26 Mar 2021 04:37:47 GMT

Hello everyone,

I want to create a bar chart as shown in the picture. Will I be able to do it using SAS. If yes, can someone help me with the syntax for the same?

Thank you for your time and consideration.

I used the following syntax but did not get the desired output as I have 2 response variables (mental&physical), also I want to know how I can add different colors?

PROC SGPLOT DATA = testrun;
VBAR category / RESPONSE = mental stat=mean barwidth=0.1 ;
TITLE 'Mental activity score';
RUN;

Blue: Category1 Red=category2 green=category3

Only want the year to appear on my x-axis

thuphammm — Fri, 26 Mar 2021 02:05:36 GMT

Hi everyone, first-time poster here!

So I have this data set:

This is my data

It produced a series plot like this:

and this is the plot I got

I want my graph to look like this with the x-axis to consist only year, obviously i failed using the SAS code below because the x-axis looks right but there's no graph.

I want the plot to have the x-axis like this (obviously I failed, the graph does not appear)this is what I used

Please help 😞

Change font and font size in proc sgpanel

mcmaxwell — Fri, 26 Mar 2021 11:29:39 GMT

Hello, I'm new to SAS and am using SAS studio. Thanks in advance for the help!

I am using the following code to make a graph that represents some Likert-scale data I collected:

ods graphics / reset width=6in height=3in imagename='Likert_Panel_';  
title 'Student perceptions of the social aspects of cPLTL';  
proc sgpanel data=all;  
	styleattrs datacolors=(deeppink hotpink pink MistyRose lightSalmon salmon Tomato)   
		datacontrastcolors=(black);  
	panelby question / layout=panel columns=1 onepanel novarname noborder nowall;  
	highlow y=cPLTL low=low high=high / group=group type=bar nooutline 
		lowlabel=sumdisagree highlabel=sumagree dataskin=matte;  
    rowaxis display=(nolabel noticks ) fitpolicy=none;  
colaxis display=(nolabel);    
    keylegend / noborder;  
    where category = "social"; 
run;

The resulting image looks like this:

Is there a way to change the font style and font size of the whole graph to make it easier to read? Also, is there a way to get rid of the boxes around each panel? Thanks!

Graficar

Uriel26 — Wed, 24 Mar 2021 17:12:34 GMT

Alguien podría ayudarme, quiero hacer una gráfica horizontal que me permita manejar dos conjuntos de datos en la misma gráfica,por ejemplo quiero tener alturas, diámetros y que estén en función de número de sitio

Waterfall

SASSLICK001 — Wed, 24 Mar 2021 00:06:50 GMT

How can i draw following graph using sas? Any ideas would be appreciated. 001, 002 are patient numbers, the colored lines are a criterial presenting for each patient.

Stacked bar graph of column percentages for subgroups

Jason2020 — Tue, 23 Mar 2021 23:27:53 GMT

Greetings,

I have a data set of a survey, and the program below generates a stacked bar graph of the percentage of DISEASE status (3 categories) by LOCATION. However, the graph uses the percentages of the total sample, whereas I want the graph to show the column percentages, and only for diseased individuals.

I tried using WHERE statement to include only the diseased groups, but the percentages becomes wrong. Is there a way around that?

ods graphics on;

proc surveyfreq data=moroc6 order=formatted;

tables disease * location /col nowt nototal plots= freqplot(groupby=column twoway=stacked scale=percent);

where disease in ("Acute" , "Chron" , "None");

cluster province;

strata region;

weight sampling_weight;

run;

ods graphics off;

line shift on y-axis tick label

Claus_Stenberg — Mon, 22 Mar 2021 12:32:52 GMT

I am doing an template for an AE dotplot for clinical trials. The text string for a preferred term is shown on the y axis. Often the text string is very long why I want have a line break so the term will be on two lines (the y axis tick mark label). I have tried to get force a line shift in the text string by setting in the unicode for a line break - (*ESC*){unicode '000a'x} - but that does not work.

I am doing the plot in GTL

proc template;
   define statgraph aedotplot;
   
      begingraph ;

    
         layout lattice /columns=5 
                         columnweights = (0.4 0.1 0.1 0.2 0.2) 
                         rowdatarange=unionall 
                         columngutter=0 
                         ;
                         

            layout overlay / walldisplay=none 
                             border=false 
                             x2axisopts=(display=(tickvalues ) 
                             displaysecondary=(line)) 
                             yaxisopts=(linearopts=(tickvaluesequence=(start=1 end=&count increment=1))
                                       display=none
                                       displaysecondary=none
                                       );

      /*--Draw alternating bands using referenceline--*/
               referenceline y = display2nd/ lineattrs=(thickness=&RefThickness.PCT) datatransparency=0.9;
/*             left colum with text box */
               scatterplot y=display x=display_label / markercharacter=display MARKERCHARACTERPOSITION=right yaxis=Y xaxis=X2 markerattrs=(size=0) includemissinggroup=true;

            endlayout;
......

Any ideas ?

Remove x-axis labels proc tree

Lollieboo — Sun, 21 Mar 2021 18:05:55 GMT

I have some data and code and I am producing a dendrogram using the proc tree statement as I have been instructed to do so.

However, my dataset has multiple thousands of observations, and printing thousands of labels for each observation is visually impractical. I would like to remove the labels from the x-axis of the dendrogram by amending the proc tree code, but I can't figure out how. The only thing I could find online to try was haxis=none, but that did not work.

Any advice?

Keylegend Elements SGPLOT

Anita_n — Sat, 20 Mar 2021 14:51:17 GMT

Dear all,

I have a plot that looks like this:

I wish to insert the percentages in the graph and the keylegend should just be one or two lines outside the graph. I used across=2 but it isn't working. Any help?

proc sgplot data = mydata dattrmap= myattrmap pad=(bottom=5%) noautolegend noborder;
xaxis grid type=discrete discreteorder=data;
vbar year / response=frequencies group = mygroup attrid=mygroup barwidth=0.9 missing ;
xaxis  display= (nolabel) valueattrs=(family="verdana" size=8pt) fitpolicy=rotate ValuesRotate=Vertical;
yaxis  display= (nolabel) valueattrs=(family="verdana" size=8pt);
keylegend / location=outside position=bottomright across=2 title="" valueattrs=(size=10 family="verdana") outerpad=(top=0.2cm) noborder;
run;

Use SGPLOT to produce an animated .svg via the HTML or HTML5 ODS destinations

DaveShea — Sat, 20 Mar 2021 08:57:44 GMT

Hi,

I would like to produce a chart animation in a single .html file that contains a nice and simple Scaleable Vector Graphic (.svg) file that is the result of running Proc SGPLOT.

I seem to be able to do this with Proc GCHART, but I definitely want to use Proc SGPLOT.

I only seem to be able to get SGPLOT and the .svg parts working if I use the ODS PRINTER destination and that produces a .svg file, rather than an .html file.

I really do want a single .html file that, when opened in a browser, shows my nice animation, and crucially, displays the pause/refresh controls to allow the user to manage the .svg animation.

My demo code via the ODS PRINTER destination, using SASHELP.CARS, so you can all join in, looks like the following.

You will probably need to amend the path used in the ODS PRINTER File option.

GOptions Reset=ALL;
Options  PrinterPath=SVG 
         Animate=START 
         AnimDuration=0.75 
         PaperSize=("45cm", "15cm")
         SVGTitle="Daves Animation"
         ;

**************************************************************;
*Close ALL ODS destinations before we start.                  ;
**************************************************************;
ODS _ALL_ CLOSE;
**************************************************************;
*Open just the ODS PRINTER destination.                       ;
**************************************************************;
ODS PRINTER File="C:\Temp\SGPlotDemo03.svg";
**************************************************************;
*Use each SGPLOT to produce a frame of the animation.         ;
**************************************************************;
title1 "Plot of MPG City versus Horsepower - Cars from Asia" ;  
proc sgplot data=sashelp.cars; 
Where Origin="Asia";
  bubble x=horsepower y=mpg_city size=cylinders / fillattrs=(color=green transparency=.5);
run;

title1 "Plot of MPG City versus Horsepower - Cars from Europe" ;  
proc sgplot data=sashelp.cars; 
Where Origin="Europe";
  bubble x=horsepower y=mpg_city size=cylinders / fillattrs=(color=yellow transparency=.5);
run;

title1 "Plot of MPG City versus Horsepower - Cars from USA" ;  
proc sgplot data=sashelp.cars; 
Where Origin="USA";
  bubble x=horsepower y=mpg_city size=cylinders / fillattrs=(color=red transparency=.5);
run;
Title;

**************************************************************;
*Stop the animation.                                          ;
**************************************************************;
Options Animation=STOP;
**************************************************************;
*Close the ODS PRINTER destination.                           ;
**************************************************************;
ODS PRINTER CLOSE;

I hope that someone can put me out of misery here.

Am I chasing ghosts trying to get an .svg animation into a single .html file using Proc SGPLOT or is there something fundamental that I am missing here ? I do not understand why Proc GCHART plays ball writing out to the ODS HTML destination, but Proc SGPLOT does not.

Many thanks,

Downunder Dave

Proc SGPLOT - Legend

A_M_ — Fri, 19 Mar 2021 11:48:13 GMT

Hello,

I am using sgplot to do 2 histograms in a figure, comparing the distribution of 2 groups (one in blue, one in red). I chose to keep the histogram bars empty (nofill option). The legend therefore shows only a colored line, one red and one blue. I would however like to have boxes in the legend, not lines. To match the appearance of my histogram bars, I would like to therefore have empty/white boxes, with colored outlines matching my groups. I can't seem to figure out how to do this. I can either have just lines, or fully colored red and blue boxes. Any experience/ideas with that? Thanks a lot!

Changing the bar color in high low plot

Leo9 — Tue, 16 Mar 2021 23:51:11 GMT

Hi All,

Is there a way to change the bar color in the high-low plot without using attrid and dattrmap ? I am trying to plot high-low where their are three bars and I want to sign unique color to each bar based on treatment. Your help is appreciated.

How do I change the datalabel title in SGPLOT VBAR?

marta25 — Tue, 16 Mar 2021 15:28:37 GMT

Hello,

I am trying to show data labels on the bar chart like in the image.

My code looks like:

proc sgplot data=mydata
	styleattrs datacolors=(blue red) datalinepatterns=(solid);
	vbar season_ / response=quantity
        group=condition limitstat=stderr numstd=1 stat=mean
	datalabel=quantity datalabelattrs=(size=12pt) datalabelpos=bottom
	limitattrs=(color=black thickness=1) dataskin=pressed attrid=condition groupdisplay=cluster;
.....

I would like to change the label "Ev (mean)" to another name I choose (for example just "mean"), does someone know if there is a way to do it?

Thanks.

create forest plot using output from multiple models

jdr1 — Mon, 15 Mar 2021 19:39:38 GMT

Dear SAS Community,

I've been using SAS for several years and hate to admit that I typically cut and paste model results and graph or create tables using a different program. But I've decided to take my SAS graphing skills farther.

Here is the problem. I'm looking at temperature sensitive (exposure) in a severe maternal morbidity (outcome) across 5 critical windows of pregnancy: (1) pre-conception, (2) Trimester 1, (3) Trimester 2, (4) Trimester 3, and (5) last month of pregnancy. For each critical period, I'm modeling relative risk for 4 separate exposure outcomes : very cold vs moderate, cold vs moderate, very hot vs moderate, and hot vs moderate across three temperature variables Average temperature, maximum temperature, and minimum temperature.

For example, results examining average temperature exposure and SMM relative risk across 5 critical periods will produce results from 20 models. How do I create a forest plot for maximum temperature and illustrate Relative Risk (RR) for 4 separate exposure outcomes: very cold vs moderate, cold vs moderate, very hot vs moderate, and hot vs moderate across 5 critical periods of pregnancy?

Here is example code that produces RR for SMM in response to 4 exposure categories of maximum temperature:

**LMvcold_mod LMvhot_mod LMcold_mod LMhot_mod;

proc genmod data= SMMheat.Tmaxdup;

*by preterm race_ethn;

class RFA_ID ;

model SMM20 (Event='1')= LMvcold_mod / dist=poisson type3 link = log dscale;

repeated subject= RFA_ID / type=exch modelse;

estimate 'Beta LMvcold_mod' LMvcold_mod 1 / exp;

run;

proc genmod data= SMMheat.Tmaxdup;

*by preterm race_ethn;

class RFA_ID ;

model SMM20 (Event='1')= LMvhot_mod / dist=poisson type3 link = log dscale;

repeated subject= RFA_ID / type=exch modelse;

estimate 'Beta LMvhot_mod' LMvhot_mod 1 / exp;

run;

proc genmod data= SMMheat.Tmaxdup;

*by preterm race_ethn;

class RFA_ID ;

model SMM20 (Event='1')= LMcold_mod / dist=poisson type3 link = log dscale;

repeated subject= RFA_ID / type=exch modelse;

estimate 'Beta LMcold_mod' LMcold_mod 1 / exp;

run;

proc genmod data= SMMheat.Tmaxdup;

*by preterm race_ethn;

class RFA_ID ;

model SMM20 (Event='1')= LMhot_mod / dist=poisson type3 link = log dscale;

repeated subject= RFA_ID / type=exch modelse;

estimate 'Beta LMhot_mod' LMhot_mod 1 / exp;

run;

CDF plots (cumulative density) using proc univariate

Manj — Mon, 15 Mar 2021 17:02:50 GMT

Hi, I using the following code to build the plot attached. I would need to change the attributes of the lines to pink solid line and green dashed line. At present it is showing blue and red solid lines be default. Anyone suggest a method to change the attributes?

ods listing close;
ods pdf file="&g_outfile..&g_fontsize" nobookmarkgen style=idsl dpi=600 ;

ods noproctitle ;

title1 j=l height=10pt "Protocol: %upcase(&g_study_id)" j=r "Page 1 of 1";
title2 j=l height=10pt "Population: &g_poplbl";
title3 j=center "Figure &g_dsplynum";
title4 j=center "&g_title1 ";
title5 %if &g_title2 = %then %do; color = white "." %end; %else %do; j=center "&g_title2." %end;;

*footnote1 j=left h=9pt "&g_foot1";
Footnote1 j=left h=9pt j=left "&g_userid: &g_pgmpth &sysdate9 &systime";

proc univariate data=sf36_obs noprint ;
class trtppn ;
var chg;
cdfplot /overlay odstitle=none ;
run;

ods pdf close;
ods listing;

SGPanel with different x-axis range

grezek_tcfbank_com — Sat, 13 Mar 2021 14:37:11 GMT

I am trying to get SGPanel to produce a scatter plot graphs in 4 panels with different x-axis scales. I divide the data into 4 groups of x-values--the by statement divides the data such that each panel displays the entire common range (see attach). Is it possible to have each panel with its own range such that the detail within the range is clearer.

The code I use:

ods rtf body="K:\MRM\Model ...\graphs\AVM vs. True Scatterplot by &name..rtf";
proc sgpanel data=&lib..reg ;
panelby &by. / /*uniscale=row*/ sort =data ;
* scatter x=lastsaleprice y=avm_value / ;
reg x=lastsaleprice y=avm_value / cli /*nomarkers*/ lineattrs=(color=red);
format avm_value cpred dollar12.0;
*colaxis min=0 max=1750000;
* rowaxis min=0 max=1750000;
TITLE "AVM vs. True Scatterplot by &name";
run;
ods rtf close ;

The image I'm getting is attached. Any suggestions are welcome. Thanks for your help. -- George

sgplot stack bar by descending order while keeping group order inside each bar

garfieldyue — Fri, 12 Mar 2021 20:42:26 GMT

Hi,

I’m having trouble in the following bar plot.

My goals are:
- Descending order by the total count of each key word
- The same order of two teams inside each bar

The code I’m using now is as below:

proc sort data=rct; 
  by descending totalct team;
run;

proc sgplot data = rct;  
                                    
  VBAR keyword /RESPONSE = count  /*y-axis: count per keyword per team*/  
                  group = team    /*stack each bar with two teams*/
                  categoryorder = respdesc; /*order by count of each team per keyword, not the total count of each keyword*/
       
run;

Part of the data is like below:

I know the above code is not 100% correct in terms of achieving my goals, but I don’t know what code (proc or options) to use in this case.

Please advise.

Thanks!

Alicia

Waterfall Plot: x-axis tick mark vertical labels

Banoo — Fri, 12 Mar 2021 15:43:20 GMT

Hello,

In my waterfall, I am plotting change in tumor size for 145 subjects. As a result the subject numbers are not displayed clearly and x-axis is not readable. Is there a way to make the subject numbers (x-axis) tick marks vertical? Attached to this post is an image of x-axis from the plot. Any help is much appreciated. Thanks!!

Space out clustered bar chart more

lrackley — Thu, 11 Mar 2021 18:57:32 GMT

Can anyone tell me which options I have to change to make these bars bigger and more space between them? Just looks too on top of each other as it is currently. I tried barwidth and clusterwidth and it's not producing what I want.

proc sgpanel data=aedata;
panelby sev_abbr;
hbar SOC/response=count group=aesevn groupdisplay=cluster datalabel=labelvar clusterwidth=0.85 barwidth=0.85;
colaxis label="Number of Events";
rowaxislabel="System Organ Class";
run;

How to change color scale of a Heatmap in PROC SGPLOT?

gmgrosse — Wed, 10 Mar 2021 07:46:52 GMT

Dear all,

I try to create a Heatmap of inter-biomarker correlations using PROC SGPLOT. However, I repeatedly failed to adjust the color scale. I would like to define white as the "0"-level of Pearson's r, and blue for negative and red for positive values, as I think that this would be the intuitive Interpretation. See the code below. I tried to define colormodel(STARTCOLOR "-1") etc. but this was not helpful at all.

I am grateful for any suggestions.

Thank you and kind regards

proc sgplot data=CorrLong noautolegend;

heatmap x=Variable y=CorrelationID / colorresponse=Correlation name="nope1" discretex discretey x2axis colormodel=ThreeColorRamp;

text x=Variable y=CorrelationID text=p_value / textattrs=(size=5pt) x2axis name='nope2';

label correlation='Pearson Correlation';

yaxis reverse display=(nolabel);

x2axis display=(nolabel);

gradlegend;

run;

Plot multiple x-y relations with 3D-tilt

Sinistrum — Tue, 09 Mar 2021 14:34:30 GMT

Dear community,

I would please like to ask for your help concerning the following issue.

In a nutshell, I want to create a graph akin to this:

That is, I have, as a data set, a discrete approximation of n probability-density functions (here, n=3).

I want to plot them such that they do not appear "on top of each other", but with a 3D-tilt.

Plotting "on top of each other" I achieved via proc gplot. However, I do not know how to tackle my current issue.

It would be great, if you could please provide help with respect to how achieve my goal.

Sample-Code is embedded in the spoiler;

data	in1;
	input	x	t1 t2 t3;
	datalines;
		0.05	0.1	0.05	0.001
		0.1	0.2	0.1	0.001
		0.15	0.2	0.1	0.001
		0.2	0.3	0.4	0.01
		0.25	0.4	0.4	0.02
		0.3	0.5	0.5	0.05
		0.35	1	1	0.1
		0.4	2	2	3
		0.45	2.5	3	4
		0.5	3	4.5	6
		0.55	2.5	3	4
		0.6	2	2	3
		0.65	1	1	0.1
		0.7	0.5	0.5	0.05
		0.75	0.4	0.4	0.02
		0.8	0.3	0.4	0.01
		0.85	0.2	0.1	0.001
		0.9	0.2	0.1	0.001
		0.95	0.1	0.05	0.001
	;
run;

	symbol1
		interpol		=			join 
		value			=			DiamondFilled 
		color			=			black 
		height			=			0.25;      
	symbol2
		interpol		=			join 
		value			=			DiamondFilled 
		color			=			blue 
		height			=			0.25;                                                                         
	symbol3
		interpol		=			join 
		value			=			DiamondFilled 
		color			=			orange 
		height			=			0.25;                                                                         
 	axis1
		label			=			("x") 
		minor			=			none 
		value			=			(color=black angle=45)
		;      
	axis2
		label			=			(angle=90 "\phi(x)")
		minor			=			(n=1);  
	proc	gplot
		data	=	in1
					;
		plot (t1 t2 t3)* x
		/
		overlay
			haxis	=	axis1
			vaxis	=	axis2
		;
	run; quit;

Yours sincerely,
Sinistrum

Histogram options - Specify Color bar & Title/Description

IChatterji — Fri, 05 Mar 2021 18:26:12 GMT

Hello,

I am using the following code to create Histograms that I will ultimately include in a presentation package. I need to format the titles and color bars so that it matches the rest of the presentation template. My 2 needs are :

- What is the code to specify the color that is filled in the histogram bar?

- What is the code to display the BY variables in the Plot Description? Currently the histogram description is "Distribution of PnL". I want to be able to say "Distribution of PnL - <Symbol>_<RiskProfile>"

Any suggestions would be greatly appreciated. If this question should be posted in Graphics community, please let me know.

ods select Histogram;
proc univariate data=Mstr_rule4a_analysis noprint;
by symbol riskprofile;
var PnL;
histogram/barlabel=count ;
where win_loss ne 'FEE' and year ge 2015 ;
inset n mean std="Std Dev" / pos = ne format = 5.2;
title 'IB One-5yrs-Histogram PnL';
run;
title;

Ask the Expert Webinar: How Can I Create Graphs Using SAS®

wreeves — Wed, 24 Feb 2021 21:04:04 GMT

Hi Graphics Programming Community,

SAS customers Richann Watson and Kriss Harris are presenting a live “Ask the Expert” webinar on March 4^th, 11 AM – Noon ET entitled How Can I Create Graphs Using SAS^®? Learn to harness the power of SAS to create meaningful graphs using Statistical Graphics procedures and Graph Template Language (GTL). Richann and Kriss will share their expertise in creating SAS graphs and show you how to customize graphs for your audience. They’ll also host a live Q&A session at the end of the webinar.

You will learn:

How to use Statistical Graphics procedures.
Ways to adjust your graph font, color, and more.
Ways to use ODS output objects to get data from another SAS procedure.
Techniques using Graph Template Language (GTL) to create a custom graph.

Click here to register for the webinar.

Want more tips? Be sure to subscribe to the Ask the Expert board to receive follow up Q&A, slides and recordings from this and other SAS Ask the Expert webinars.

Can't join the live event? You can view this and other Ask the Expert sessions on-demand here.

Proc Gchart + Values dont Fit inside Stack Bar

pmksasbi — Sat, 20 Feb 2021 19:22:25 GMT

Hi Team,

i am working on creating Stacked Vertical Bar chart. and when my values are small they don't fit inside the bar. can you please help on this.

The code run good to this values:

data sold;

input city $ Sales Type$;

if type = 'customer' then color = '1';

else if type = 'employer' then color = '3';

else if type = 'manufact' then color = '2';

else color = '4';

datalines;

current 40 customer

previos 25 manufact

previos 10 customer

previos 5 employer

;

run;

In the below code employer values are not shown inside the bar. Also i would like to have the dollar values inside the midddle of the bar instead on the top.

Here is the code :

data sold;

input city $ Sales Type$;

if type = 'customer' then color = '1';

else if type = 'employer' then color = '3';

else if type = 'manufact' then color = '2';

else color = '4';

datalines;

current 20 customer

previos 12 manufact

previos 7 customer

previos 1 employer

;

run;

goptions reset = all border cback = white hsize=6in htitle =2.0 htext =5.0 imagestyle = vsize=3in;

goptions device=jpeg gsfname=out;

goptions xpixels=550 ypixels=650;

ods graphics on / width =2in;

ods graphics / height =20in;

ods listing close;

ods html path = 'mypath' body="test.htm" style =htmlblue gpath = '/mypath/test.jpg';

pattern1 v =solid color='#33ccff';

pattern2 v =solid color='#ff9900';

pattern3 v =solid color='#bfbfbf';

proc gchart data=sold;

vbar city/discrete sumvar=sales nozero group = city subgroup = color type = sum

width = 8 space =15 NOFRAME noaxis noloegend coutline = gray99 woutline =1 inside=sum;

format sales dollar8.

run;

quit;

ODS HTML CLOSE;

ODS LISTING;

Change Error Bar Line Style by Group

moreka — Tue, 23 Feb 2021 17:10:10 GMT

Is there anyway to alter the style of error bars by group in regards to line style or thickness?

I'm working on a graph where the error bars overlap, so I want them to appear distinct on the graph.

I've attempted with SGPLOT and SGDESIGN, but have only been able to differentiate the groups by color.

I do not want to stagger the bars using the cluster option.

These are examples of the code I'm working with.


proc template;
define statgraph sgdesign;
dynamic __XCLAS __PREDICTED __UCLM __LCLM __GROUP2;
begingraph / dataskin=none attrpriority=Auto DataColors=(CX000000 CX848284) DataContrastColors=(CX000000 CXA4A5A4) DataLinePatterns=(4 1) DataSymbols=(circlefilled trianglefilled);
   layout lattice / rowdatarange=data columndatarange=data rowgutter=10 columngutter=10;
      scatterplot x=__XCLAS y=__PREDICTED / group=__GROUP2 yerrorupper=__UCLM yerrorlower=__LCLM name='scatter' groupdisplay=Overlay errorbarattrs=(thickness=1);
   endlayout;
endgraph;
end;
run;

proc sgrender data=B.EFFECTPLOTRAWDATA1 template=sgdesign;
dynamic __XCLAS="'_XCLAS'n" __PREDICTED="'_PREDICTED'n" __UCLM="'_UCLM'n" __LCLM="'_LCLM'n" __GROUP2="'_GROUP'n" ;
run;

proc sgplot data=b.effectPlotRawData1;
scatter x = _XCLAS y = _PREDICTED / yerrorupper=_UCLM yerrorlower=_LCLM group=_INDEX name='scatter' ERRORBARATTRS=(THICKNESS = 1);
styleattrs datacontrastcolors=(black darkgray) datalinepatterns=(ShortDash LongDash);
;run;

How to use a line chart instead of a barchart

Scooby3g — Tue, 23 Feb 2021 16:51:38 GMT

Hi,

I'm new to using charts in sas. I want to have a line chart and a table below the chart. I used code from the following reference from the link: https://support.sas.com/kb/49/697.html but instead of using a bar chart. I want to use a line chart so that each value for one of the variables ("Slept") is graph.

Can i replace "barchart" with something else to make this into a line graph? Thanks in advance for your thoughts and suggestions.

DATA Slept;
INPUT QUARTER SLEPT PERCENT C_Count;
DATALINES;
2019Q1	NA	5	6
2019Q2	NO	2	3
2019Q3	UKNOWN	6	8
2019Q4	YES	5	7
2020Q1	NA	5	6
2020Q2	NO	9	11
2020Q3	UKNOWN	5	6
2020Q4	YES	5	7
2019Q1	NA	0	
2019Q2	NO	17	
2019Q3	UKNOWN	0	
2019Q4	YES	83	
2020Q1	NA	0	
2020Q2	NO	14	
2020Q3	UKNOWN	14	
2020Q4	YES	71	
2019Q1	NA	0	
2019Q2	NO	100	
2019Q3	UKNOWN	0	
2019Q4	YES	0	
2020Q1	NA	0	
2020Q2	NO	0	
2020Q3	UKNOWN	0	
2020Q4	YES	100	
2019Q1	NA	0	
2019Q2	NO	0	
2019Q3	UKNOWN	0	
2019Q4	YES	100	
2020Q1	NA	0	
2020Q2	NO	0	
2020Q3	UKNOWN	0	
2020Q4	YES	100	
;
RUN;

proc template;
   define statgraph GRAPHY;
   begingraph;
      entrytitle 'SLEPT ';
      layout lattice / rowweights=(.65 .35);
         layout overlay;
	    barchart X=Quarter y=percent / group=SLEPT barlabel=true dataskin=gloss;
	 endlayout;
	 layout overlay / xaxisopts=(type=discrete label='Number per CY Quarter' 
                          display=(label)) walldisplay=none;
	    blockplot x=Quarter block=C_Count / class=Quarter display=(outline values label) 
                                              repeatedvalues=true;
	 endlayout;
      endlayout;
   endgraph;
   end;
run;
proc sgrender data=SLEPT template=GRAPHY;
run;