%newsas(smsargn);
Proc format;
     value region 1='N East'   2='N Central'   3='South'    4='West';
Title 'Crime Rate by REGION in the SMSA data - Notched Boxplots';
Data smsa;
     label crimrate ='Crime rate per 1000 population';
     format REGION region. ;
     set psy303.smsa;
     CRIMRATE = CRIME / TOTPOP;
 
Proc sort data=smsa;
     by REGION;
 
* Use Proc univariate to get median and HSPR values for notched
  boxplot calculations;
 
Proc univariate data=SMSA noprint;
     by  REGION;
     var CRIMRATE;
     output out=SUMRY median=MEDIAN QRANGE=HSPR  n=N;
 
data sumry;
     set sumry;
     NOTCHWID = 1.58 * HSPR / sqrt(N);
     LONOTCH = MEDIAN - NOTCHWID;
     HINOTCH = MEDIAN + NOTCHWID;
     LOGMED =  log10(MEDIAN);
     LOGHSPR = log10(HSPR);
	  run;

proc print;
   format notchwid lonotch hinotch logmed loghspr 7.3;
*options ls=72 ps=50;
*include macros(splot);
options mprint;
%splot(data=smsa, class=REGION, var=CRIMRATE);
 
*-- Determine whether there is a need to transform the data
    to equalize variability [the plot says NO];
 
proc plot data=sumry;
     Plot LOGHSPR * LOGMED = REGION;
	  run;