Re: [R] how to compute the inverse percentile of a given observation w.r.t. a reference distribution

From: David Winsemius <dwinsemius_at_comcast.net>
Date: Wed, 25 May 2011 18:00:37 -0400

On May 25, 2011, at 5:50 PM, Dennis Murphy wrote:

> Hi:
>
> On Wed, May 25, 2011 at 12:42 PM, rudi <rudi.strasser@gmail.com>
> wrote:
>> Hi,
>>
>> can anyone help me to figure out how to compute the percentile of an
>> individual observation with respect to a reference distribution.
>>
>> What I mean is. Let's assume I have a vector consisting of 10 numbers
>> {3,5,8,1,9,5,4,3,5.5,7} and I want figure out what percentile the
>> number 4.9 corresponds to. I failed to find any reference to such a
>> function, although I would assume this must frequently be necessary.
>
> The simple answer is, I believe,
>
> x <- c(3,5,8,1,9,5,4,3,5.5,7)

Try instead:

ecdf(x)(4.9)
> [1] 0.4

ecdf returns a function, so why not use it as such? It is also linked from the quantile help page where it is called the "inverse of quantile".

-- 
David.

> plot(ecdf(x))
> sum(x <= 4.9)/length(x)
> [1] 0.4
>
(Somewhat more complicated than necessary.)
> This would correspond to the empirical cumulative distribution
> function (ecdf) to which Jorge alluded.
>
> HTH,
> Dennis
>>
>> Thanks in advance for you help.
David Winsemius, MD West Hartford, CT ______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Received on Wed 25 May 2011 - 22:03:35 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 25 May 2011 - 22:50:09 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive