Hi all,
Could you please help to code the 'Standard Error' for the gsplot. The plot has to be as below. What should be a standard of the dataset? What are the options for proc sgplot?
Try something like:
proc sql;
create table classGraph as
select
sex,
age,
mean(weight) as meanWeight,
mean(weight) - std(weight) as lowStdWeight,
mean(weight) + std(weight) as highStdWeight
from sashelp.class
group by sex, age;
quit;
proc sgplot data=classGraph;
scatter x=age y=meanWeight /
yerrorlower=lowStdWeight yerrorupper=highStdWeight group=sex;
series x=age y=meanWeight / lineattrs=(pattern=solid) group=sex;
run;
The main options you want to use on the VLINE statement are the following:
1. LIMITSTAT = CLM | STDERR | STDDEV -- for your case, use STDERR
2. NUMSTD = <number> -- This is a multiplier for STDERR or STDDEV. The default is 1, which mean you would get one standard error in your case.
3. LIMITS = UPPER | LOWER | BOTH -- the default is BOTH, just like in your picture.
Here is a simple example:
proc sgplot data=sashelp.class;
vline age / response=weight stat=mean limitstat=stderr numstd=2;
run;
You can also modify the visual attributes of the error bars using the ERRORATTRS option.
Hope this helps!
Dan
Thank you for the reply! Could you please say how to modify your code to draw two lines: one for 'male', another for 'female'.
And another question: if I have SE calculated in previous steps (this is an obligate rule) how to put SE on the plot in this case?
Try something like:
proc sql;
create table classGraph as
select
sex,
age,
mean(weight) as meanWeight,
mean(weight) - std(weight) as lowStdWeight,
mean(weight) + std(weight) as highStdWeight
from sashelp.class
group by sex, age;
quit;
proc sgplot data=classGraph;
scatter x=age y=meanWeight /
yerrorlower=lowStdWeight yerrorupper=highStdWeight group=sex;
series x=age y=meanWeight / lineattrs=(pattern=solid) group=sex;
run;
What does your data look like?
Perhaps something like this example?
Thank you for the last link. I used it:
invalue defines an INFORMAT. You need a VALUE statement to define a format in proc format.
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 16. Read more here about why you should contribute and what is in it for you!
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.
Select SAS Training centers are offering in-person courses. View upcoming courses for: