[R] Finding percentile of a value from an empirical distribution
Robert A'gata
rhelpacc at gmail.com
Wed Jan 11 06:43:37 CET 2012
Thank you. That's easier than I thought.
On Wed, Jan 11, 2012 at 12:06 AM, Jorge I Velez
<jorgeivanvelez at gmail.com> wrote:
> Hi Robert,
>
> Try
>
> set.seed(123)
> x <- seq(100)
> x <- sample(x, 1000, replace = TRUE)
> f <- ecdf(x)
> f(10)
> # [1] 0.099
> f(71)
> # [1] 0.716
>
> See ?ecdf for more information.
>
> HTH,
> Jorge.-
>
>
> On Tue, Jan 10, 2012 at 11:52 PM, Robert A'gata <> wrote:
>>
>> Hello,
>>
>> I am not sure how to do this in R. Any suggestion would be
>> appreciated. I have a vector of values from where I build an empirical
>> CDF. For example:
>>
>> > x <- seq(1,100)
>> > x <- sample(x,1000,replace=T)
>> > quantile(x,probs=seq(0,1,.05))
>> 0% 5% 10% 15% 20% 25% 30% 35% 40% 45%
>> 50% 55%
>> 1.00 5.00 10.00 16.00 20.00 25.00 31.00 36.00 41.00 45.55
>> 50.00 56.00
>> 60% 65% 70% 75% 80% 85% 90% 95% 100%
>> 60.00 65.00 70.00 74.00 80.00 85.00 91.00 95.05 100.00
>>
>> I would like to write a function that takes in a number z and vector x
>> (i.e. the raw vector).It returns percentile of z wrt x. E.g.
>>
>> > f(71,x)
>>
>> Should return something around 0.708 or 0.709. I am wondering if there
>> is any pre-packaged functions that do this in R? If not, how can I
>> write such a function. Any suggestion would be appreciated.
>>
>> Robert
>>
>> ______________________________________________
>> 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.
>
>
More information about the R-help
mailing list