data forest;
input Study $1-16 grp OddsRatio LowerCL UpperCL Weight;
format weight percent5. Q1 Q3 4.2 oddsratio lowercl uppercl 5.3;
ObsId=_N_;
OR='HR'; LCL='LCL'; UCL='UCL'; WT='Weight';
if grp=1 then do;
Q1=OddsRatio-OddsRatio*weight;
Q3=OddsRatio+OddsRatio*weight;
lcl2=lowercl;
ucl2=uppercl;
end;
else study2=study;
datalines;
Outcome1 . . . . .
AC N=60787 1 0.62 0.53 0.73 1
AX N=33680 1 0.60 0.49 0.76 1
AV N=19766 1 0.61 0.49 0.75 2
AD N=7420 1 0.62 0.51 0.74 1
Outcome2 . . . . .
AC N=60738 1 0.51 0.40 0.65 1
AX N=33830 1 0.48 0.35 0.67 1
AV N=19736 1 0.55 0.39 0.76 1
AB N=7240 1 0.50 0.38 0.66 1
;
run;
proc sort data=forest out=forest2;
by descending obsid;
run;
/* Add sequence numbers to each observation */
data forest3;
set forest2 end=last;
retain fmtname 'Study' type 'n';
studyvalue=_n_;
if study2='Overall' then study2value=1;
else study2value = .;
/* Output values and formatted strings to data set */
label=study;
start=studyvalue;
end=studyvalue;
output;
if last then do;
hlo='O';
label='Other';
end;
run;
/* Create the format from the data set */
proc format library=work cntlin=forest3;
run;
/* Apply the format to the study values and remove Overall from Study column. */
/* Compute the width of the box proportional to weight in log scale. */
data forest4;
format studyvalue study2value study.;
drop fmtname type label start end hlo pct;
set forest3 (where=(studyvalue > 0)) nobs=nobs;
if studyvalue=1 then studyvalue=.;
/* Compute marker width */
x1=oddsratio / (1 ** (weight/2));
x2=oddsratio * (1 ** (weight/2));
/* Compute top and bottom offsets */
if _n_ = nobs then do;
pct=0.3/nobs;
call symputx("pct", pct);
call symputx("pct2", 2*pct);
call symputx("count", nobs);
end;
run;
ods listing close;
ods html image_dpi=100 path="." file='sgplotforest.html';
ods graphics / reset width=600px height=400px imagename="Forest_Plot_Vector" imagefmt=gif;
title "";
title2 h=8pt 'Hazard Ratio and 95% CL';
proc sgplot data=forest4 noautolegend;
scatter y=study2value x=oddsratio / markerattrs=graphdata2(symbol=diamondfilled size=10);
scatter y=studyvalue x=oddsratio / xerrorupper=ucl2 xerrorlower=lcl2 markerattrs=graphdata1(symbol=squarefilled size=0);
vector x=x2 y=studyvalue / xorigin=x1 yorigin=studyvalue lineattrs=graphdata1(thickness=8) noarrowheads;
scatter y=studyvalue x=or / markerchar=oddsratio x2axis;
scatter y=studyvalue x=lcl / markerchar=lowercl x2axis;
scatter y=studyvalue x=ucl / markerchar=uppercl x2axis;
refline 1 2 / axis=x;
refline 1 1 / axis=x lineattrs=(pattern=shortdash) transparency=0.5;
xaxis type=log offsetmin=0 offsetmax=0.35 min=0.3 max=1.2 minor display=(nolabel) ;
x2axis offsetmin=0.7 display=(noticks nolabel);
yaxis display=(noticks nolabel) offsetmin=0.2 offsetmax=0.1 values=(1 to &count by 1);
run;
ods html close;
ods listing;Hi everyone,
I am trying to modify the above SAS code but running into a few errors.
- not all data points are displayed.
-I am trying to spell the subgroup completely but adjusting the group length to 1-35 creates issues with the graph
-For the list of groups I am trying to make them aligned to the left and increase the space for the subgroup names with the listed N. Please find a screen shot below
Accepted Solutions
SAS INNOVATE 2024
Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.
If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.
Back in the Classroom!
Select SAS Training centers are offering in-person courses. View upcoming courses for:
-
462 views
-
1 like
-
-
