/* gatherered temp data from Weather Underground */
/* historical data like this: */
/* http://www.wunderground.com/history/airport/KRDU/2009/1/17/DailyHistory.html?req_city=NA&req_state=NA&req_statename=NA */
footnote "Data source: http://www.wunderground.com";
data parties;
length date 8 hi 8 lo 8;
format date year4.;
informat date date9.;
label lo = "Low temp"
      hi = "High temp"
	  date = "Party year";
infile datalines dsd;
input date hi lo;
datalines;
14JAN1995,73,61
13JAN1996,45,24
11JAN1997,44,23
10JAN1998,55,30
09JAN1999,62,35
15JAN2000,44,19
13JAN2001,55,33
12JAN2002,51,24
18JAN2003,32,16
17JAN2004,51,22
15JAN2005,42,33
14JAN2006,59,35
13JAN2007,69,48
12JAN2008,58,41
17JAN2009,36,11
;

/* REQUIRES SAS 9.2 for ODS graphics and SGPLOT */
ods graphics /width=600 height=350;
proc sgplot data=parties;
vbar date /response=hi transparency=.8 fillattrs=graphdata2 ;
vbar date /response=lo barwidth=.5 fillattrs=graphdata1;
yaxis label="Temps on party day";
refline 30 /axis=y label="Jan avg low";
run;

data parties;
set parties;
lo_c = (lo-32) * (5/9);
run;

ods graphics /width=500 height=400;
proc sgplot data=parties(where=(year(date)<>1995));
vbar date /response=lo_c barwidth=.5 fillattrs=graphdata1;
yaxis label="Temps on party day (C)" min=-10 max=30;
inset "Avg Jan low: 30 degrees (F)";
run;
quit;