[R] Extracting values from a ecdf (empirical cumulative distribution function) curve
Duncan Mackay
dulcalma at bigpond.com
Fri Nov 1 03:40:14 CET 2013
Hi
There is a print method for ecdf
So print(f) should give you an idea of what is going on
See ?ecdf
HTH
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mackay at northnet.com.au
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of Manoranjan Muthusamy
Sent: Friday, 1 November 2013 11:18
To: Rui Barradas
Cc: r-help at r-project.org
Subject: Re: [R] Extracting values from a ecdf (empirical cumulative
distribution function) curve
Thank you, Barradas. It works when finding y, but when I tried to find x
using interpolation for a known y it gives 'NA' (for whatever y value). I
couldn't find out the reason. Any help is really appreciated.
Thanks,
Mano
On Thu, Oct 31, 2013 at 10:53 PM, Rui Barradas <ruipbarradas at sapo.pt> wrote:
> Hello,
>
> As for the problem of finding y given the ecdf and x, it's very easy,
> just use the ecdf:
>
> f <- ecdf(rnorm(100))
>
> x <- rnorm(10)
> y <- f(x)
>
> If you want to get the x corresponding to given y, use linear
> interpolation.
>
> inv_ecdf <- function(f){
> x <- environment(f)$x
> y <- environment(f)$y
> approxfun(y, x)
> }
>
> g <- inv_ecdf(f)
> g(0.5)
>
>
> Hope this helps,
>
> Rui Barradas
>
> Em 31-10-2013 12:25, Manoranjan Muthusamy escreveu:
>
>> Hi R users,
>>
>> I am a new user, still learning basics of R. Is there anyway to
>> extract y (or x) value for a known x (or y) value from ecdf
>> (empirical cumulative distribution function) curve?
>>
>> Thanks in advance.
>> Mano.
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________**________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/m
>> ailman/listinfo/r-help> PLEASE do read the posting guide
>> http://www.R-project.org/** posting-guide.html
>> <http://www.R-project.org/posting-guide.html>
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
[[alternative HTML version deleted]]
______________________________________________
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