[R] "normal probability plot" with a percentile scale?

[Peter Dalgaard BSA]

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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._