SAS and R Scripts to Compute Normal Scores

1. R script to compute normal scores for sample of size n
   using qnorm (normal quantiles) function.

p = seq(0.1, 0.9, 0.1)
print(qqnorm(nsco))

2. SAS script to compute normal scores for sample of size n
   using probit (inverse normal) function.

data nsco;
   keep p nsco;
   do p = 0.1 to 0.9 by 0.1;
      nsco = probit(p);
      output;
   end;
   
proc print;

3. SAS script to compute normal scores for a sample of size n
   using proc ranks.  Use Van der Waerden's method of obtaining 
   normal scores.
   
data orig;
   keep p;
   do p = 0.1 to 0.9 by 0.1;
      output;
   end;
   
proc print;

proc rank data=orig out=nscores normal=vw;
   var p;
   ranks nsco;