%newsas(draftbox);
%include goptions;
%include data(draftusa);

goptions hsize=6in vsize=7in htext=1.8 htitle=2.2;

data draftusa;
   set draftusa;
	birthday = intck('day', mdy(1,1,96), mdy(month,day,96));
   label birthday = 'Brithday (day of year)'
         priority = 'Draft Priority value';

*-- plot priority vs. birthday;
proc gplot data=draftusa;
   plot priority * birthday /vaxis=axis1;
	symbol1 v=square h=1.5;
	axis1 label=(a=90 r=0);
run;
%gskip;

*-- boxplots by month show the pattern more clearly;
axis1 order=(1 to 12) 
	value=(t=1 'Jan' t=2 ' ' t=3 'Mar' t=4 ' ' t=5 'May' t=6 ' '
	t=7 'Jul' t=8 ' ' t=9 'Sep' t=10 ' ' t=11 'Nov' t=12 ' ') offset=(4)
	minor=none;
%boxplot(data=draftusa, class=Month, var=priority, id=day,
	classfmt=mon., varlab=Draft Priority value, connect=3, cnotch=red,
	haxis=axis1);
%gskip;

*-- ... as does a smoothed curve through the data points;
%lowess(data=draftusa, x=birthday, y=priority, gplot=YES, symbol=square,
	htext=2, step=5,silent=y);