options linesize=80;
* filename rain 'e:\Stat210\raindata.txt' ;
filename rain 'h:\HAL\Courses\Stat210\raindata.txt' ;
*;
* read in data, transform variables;
*;
data rain2;
    set rain;
    lograin = log(rainfall);
    sqrtrain = sqrt(rainfall);
**                        **;
**  sort file by group    **;
**                        **;
proc sort data=rain2;
   by treat;
**                         **; 
**  obtain normal scores   **;
**  and save output (orig  **;
**  data and normal scores)**;
**  in data file "norm"    **;
**                         **;
proc rank normal=blom data=rain2 out=norm;
  var rainfall;
  by treat;
  ranks nrm;
**                         **;
**  and produce normal     **;
**  probability plot       **;
**  for each of two groups **;
**                         **;
proc gplot data=norm;
  plot rainfall*nrm lograin*nrm sqrtrain*nrm;
  by treat;
proc ttest data=norm;
  class treat;
  var rainfall lograin sqrtrain;
run;