* Pearson Example 
  source code file normal-plot.sas
  
  Create the normal plot with proc ranks
  and plot actual data points vs. normal
  scores with sgplot;

options ls=64 nodate pageno=1;
data heights;
   infile 'c:/datasets/pearson.txt' firstobs=2;
   input father son;
   label father='Height of Father (in)'
         son='Height of Son (in)';

* Use proc rank to create normal plot
  of father's heights;
proc rank data=heights out=out normal=blom;
   var father;
   ranks nsco;

proc sgplot;
  title 'Normal Plot of Father''s Heights  Method=blom';
  scatter x=nsco y=father;

run;
quit;