R help - May 2002 - "normal probability plot" with a percentile scale?

I'd like to generate some plots like you'd see on the old "normal
probability graph paper", like the first plot in:

  <http://www.itl.nist.gov/div898/handbook/eda/section3/normprpl.htm>

except the horizontal scale would have 1%, 5%, 25%, 50%, 75%, 95%,
99%, or similar quantiles, with associated tick/grid lines.  [still
hunting around for a good example...] something like these, but with
the probability along the horizontal axis:

  <http://www.fw.umn.edu/biochr/assoc/dho/107/notes/minitab/NORMPROB.HTM>
  <http://www.survo.fi/gallery/002.html>
  <http://www.unistat.com/534.htm>
 
<http://www.mathworks.com/access/helpdesk/help/toolbox/stats/normplot.shtml>

I dug around in the mailing list archives, but didn't find anything.

Thanks!

-- 
Russell Senior         ``The two chiefs turned to each other.        
seniorr at aracnet.com      Bellison uncorked a flood of horrible       
                         profanity, which, translated meant, `This is
                         extremely unusual.' ''
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

Russell Senior <seniorr at aracnet.com> writes:
> I'd like to generate some plots like you'd see on the old
"normal
> probability graph paper", like the first plot in:
> 
>   <http://www.itl.nist.gov/div898/handbook/eda/section3/normprpl.htm>
> 
> except the horizontal scale would have 1%, 5%, 25%, 50%, 75%, 95%,
> 99%, or similar quantiles, with associated tick/grid lines.  [still
> hunting around for a good example...] something like these, but with
> the probability along the horizontal axis:
This'll be close:

x <- rnorm(100)
pp <- c(.01, .05, .25, .50, .75, .95,.99)
qq <- qnorm(pp)
qqnorm(x,axes=F)
axis(2)
axis(1,at=qq,labels=pp)
box()

For smaller samples, you may need to throw in an xlim=range(qq)
argument. 

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._