From: Peter Dalgaard BSA (p.dalgaard@biostat.ku.dk)
Date: Wed 01 May 2002 - 18:11:46 EST
Message-id: <x2sn5cjni5.fsf@blueberry.kubism.ku.dk>
Russell Senior <seniorr@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@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@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
This archive was generated by hypermail 2.1.3 : Wed 16 Oct 2002 - 11:57:14 EST