rjvbertin@hotmail.com
2002-Jun-12 17:54 UTC
Fwd: Re: Tcl8.0.4 missing function?! (PR#1640)
>From: Peter Dalgaard BSA <p.dalgaard@biostat.ku.dk> >To: "RenE J.V. Bertin" <rjvbertin@hotmail.com> >Subject: Re: Tcl8.0.4 missing function?! (PR#1640) >Date: 12 Jun 2002 19:37:21 +0200 > > > > From help(quantile): > > > > ... > > > > `quantile(x,p)' as a function of `p' linearly interpolates the > > points ( (i-1)/(n-1), ox[i] ), where `ox <- order(x)' (the ``order > > statistics'') and `n <- length(x)'. > > > > This gives `quantile(x, p) == (1-f)*ox[i] + f*ox[i+1]', where `r > > <- 1 + (n-1)*p', `i <- floor(r)', `f <- r - i' and `ox[n+1] :> > ox[n]'. > > > > I think that what is meant is 'ox <- x[order(x)]'?! When I do that, > > the values from the above formulae correspond to what quantile() gives > > me...!RenE J.V. Bertin College de France/LPPA 11, place Marcelin Berthelot 75005 Paris, France _________________________________________________________________ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel 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-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._