Hi,
I have a gplot with x var* y var, to which I need to insert few dynamic data dependent vertical lines. Need to show some events happened at different time points with a vertical dashed lines. Any idea how to get those?
Thanks
The option you are looking for is to add HREF= to your Plot statement in the options to set the locations and LHREF= to set the line type (anything other than 1 will be some sore to dot or dash line), and CHREF to set color. LHREF and CHREF also support lists so you can set different line type and color for each HREF entry if so desired.
HREF= (5, 8, 12) would set vertical reference lines at the listed values, the syntax is basically the same as for an HAXIS statement.
The fun part, since you say dynamic is where you get the values from. Easiest would be create a data set with the desired values, your choice of methods and then build a macro string such as:
Proc sql noprint;
select values into : hreflist separated by ','
from hrefdata;
quit;
and then the href code would look like
HREF=(&hreflist)
Most people use SQL or a DATA _NULL_ step to put the reference values into macro variables that you can then use to draw the reference lines. For example, the following code writes the Q1 and Q3 values of a variable to a data set and then creates macro variables with the values:
/* write quantiles to data set */
proc univariate data=sashelp.class;
var Weight;
ods output quantiles=Q;
run;
proc print data=Q;
run;
/* create macro variables for interquartile range */
data _null_;
set Q;
if Quantile="75% Q3" then call symput("Q3", Estimate);
if Quantile="25% Q1" then call symput("Q1", Estimate);
run;
%put &=Q1;
%put &=Q3;
You can then use these values in SAS procedures. For example, in PROC SGPLOT I would write:
proc sgplot data=sashelp.class;
histogram weight;
refline &Q1 &Q3 / axis=x;
run;
The option you are looking for is to add HREF= to your Plot statement in the options to set the locations and LHREF= to set the line type (anything other than 1 will be some sore to dot or dash line), and CHREF to set color. LHREF and CHREF also support lists so you can set different line type and color for each HREF entry if so desired.
HREF= (5, 8, 12) would set vertical reference lines at the listed values, the syntax is basically the same as for an HAXIS statement.
The fun part, since you say dynamic is where you get the values from. Easiest would be create a data set with the desired values, your choice of methods and then build a macro string such as:
Proc sql noprint;
select values into : hreflist separated by ','
from hrefdata;
quit;
and then the href code would look like
HREF=(&hreflist)
As an additional SGPLOT alternative, you can have the REFLINE reference a column of values for the reference line instead of using macor variables.
data airrefs;
ref='01jan1950'd;
output;
ref='01jan1960'd;
output;
run;
data air;
merge sashelp.air airrefs;
run;
proc sgplot data=air;
series x=date y=air;
refline ref / axis=x lineattrs=(pattern=dash);
run;
If you are using a BY variable for the subject, you can match-merge the reference line data to the original data by subject (if the ref data is not already in the original data) and reference the reference line data from the REFLINE statement.
Dynamic data dependent vertical lines
Not sure this is what the op wanted.
It pre-supposes the op wanted to enter some values.
see
https://goo.gl/l7zMKo
https://communities.sas.com/t5/SAS-GRAPH-and-ODS-Graphics/dynamic-data-dependent-vertical-lines/m-p/344098
HAVE
====
Up to 40 obs from sashelp.class total obs=19
Obs HEIGHT WEIGHT
1 69.0 112.5
2 56.5 84.0
3 65.3 98.0
4 62.8 102.5
5 63.5 102.5
6 57.3 83.0
7 59.8 84.5
8 62.5 112.5
9 62.5 84.0
10 59.0 99.5
11 51.3 50.5
12 64.3 90.0
13 56.3 77.0
14 66.5 112.0
15 72.0 150.0
16 64.8 128.0
17 67.0 133.0
18 57.5 85.0
19 66.5 112.0
WANT
====
Pop up window
+--------------------------+
| |
| Enter q1: _84<enter>__ |
| Enter q3: _112<enter>_ |
| |
+--------------------------+
WEIGHT
150 + A
|
|
| A
| A
125 +
|
| A A
|------------------B--------------- 84
| AA
100 + A A
|
| A
|-------AB--A--A------------------- 112
|
75 + A
|
|
|
|
50 + A
-+----------+----------+----------+
50 60 70 80
HEIGHT
FULL SOLUTION
=============
%symdel q1 q3 / nowarn;
options missing=' ';
data _null_;
length q1 q3 $10;
set mty;
if _n_=1 then do;
window start irow=1 rows=12 color=white
#3 @10 'Enter q1:' +1 q1 $10. attr = underline
#5 @10 'Enter q3:' +1 q3 $10. attr = underline;
display start;
end;
call symputx('q1',q1);
call symputx('q3',q3);
rc=dosubl('
options ls=64 ps=32;
proc plot data=sashelp.class(keep=height weight);
plot weight*height /vref=&q1 &q3;
run;quit;
');
run;quit;
1885 %symdel q1 q3 / nowarn;
1886 options missing=' ';
1887 data _null_;
1888 length q1 q3 $10;
1889 set mty;
1890 if _n_=1 then do;
1891 window start irow=1 rows=12 color=white
1892 #3 @10 'Enter q1:' +1 q1 $10. attr =
1892! underline
1893 #5 @10 'Enter q3:' +1 q3 $10. attr =
1893! underline;
1894 display start;
1895 end;
1896 call symputx('q1',q1);
1897 call symputx('q3',q3);
1898 rc=dosubl('
1899 options ls=64 ps=32;
1900 proc plot data=sashelp.class(keep=height weight);
1901 plot weight*height /vref=&q1 &q3;
1902 run;quit;
1903 ');
1904 run;
SYMBOLGEN: Macro variable Q1 resolves to 80
SYMBOLGEN: Macro variable Q3 resolves to 112
NOTE: There were 19 observations read from the data set
SASHELP.CLASS.
NOTE: PROCEDURE PLOT used (Total process time):
real time 0.03 seconds
user cpu time 0.00 seconds
system cpu time 0.03 seconds
memory 967.65k
OS Memory 16876.00k
Timestamp 03/24/2017 01:26:34 PM
Step Count 170 Switch Count 0
NOTE: There were 1 observations read from the data set WORK.MTY.
NOTE: DATA statement used (Total process time):
real time 5.90 seconds
user cpu time 0.04 seconds
system cpu time 0.09 seconds
memory 967.65k
OS Memory 16876.00k
Timestamp 03/24/2017 01:26:34 PM
Step Count 170 Switch Count 9
1904! quit;
SAS Innovate 2025 is scheduled for May 6-9 in Orlando, FL. Sign up to be first to learn about the agenda and registration!
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.
Ready to level-up your skills? Choose your own adventure.