You are not using the inv_ecdf function that Rui sent. His was
inv_ecdf_orig <-
function (f)
{
x <- environment(f)$x
y <- environment(f)$y
approxfun(y, x)
}
(There is no 'xnew' in the environment of f.)
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
From: Manoranjan Muthusamy [mailto:ranjanmano...@gmail.com]
Sent: Friday, November 01, 2013 4:38 AM
To: William Dunlap; dulca...@bigpond.com
Cc: Rui Barradas; r-help@r-project.org
Subject: Re: [R] Extracting values from a ecdf (empirical cumulative
distribution function) curve
Thanks, Bill & Duncan. Actually I tried values which are inside the defined
region. please find below the extracted script
> xnew<-rlnorm(seq(0,4000000,10000), meanlog=9.7280055, sdlog=2.0443945)
> f <- ecdf(xnew)
> y <- f(x)
> y1<-f(2000000) ## finding y for a given xnew value of 2000000
> y1
[1] 0.9950125 ## It works.
> inv_ecdf <- function(f){
+ xnew <- environment(f)$xnew
+ y <- environment(f)$y
+ approxfun(y, xnew)
+ }
## Interpolation to find xnew for a known y value.
> g <- inv_ecdf(f)
> g(0.9950125)
[1] NA
> g(0.99) ## It doesn't
[1] NA
> g(0.5)
[1] NA ## again
> g(0.2)
[1] NA ## and again
I am stuck here. Any help is appreciated.
Mano.
On Fri, Nov 1, 2013 at 2:48 AM, William Dunlap
<wdun...@tibco.com<mailto:wdun...@tibco.com>> wrote:
> it gives 'NA' (for whatever y value).
What 'y' values were you using? inf_f maps probabilities (in [0,1]) to
values in the range of the orginal data, x, but it will have problems for
a probability below 1/length(x) because the original data didn't tell
you anything about the ecdf in that region.
> X <- c(101, 103, 107, 111)
> f <- ecdf(X)
> inv_f <- inv_ecdf(f)
> inv_f(seq(0, 1, by=1/8))
[1] NA NA 101 102 103 105 107 109 111
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com<http://tibco.com>
> -----Original Message-----
> From: r-help-boun...@r-project.org<mailto:r-help-boun...@r-project.org>
> [mailto:r-help-boun...@r-project.org<mailto:r-help-boun...@r-project.org>] On
> Behalf
> Of Manoranjan Muthusamy
> Sent: Thursday, October 31, 2013 6:18 PM
> To: Rui Barradas
> Cc: r-help@r-project.org<mailto:r-help@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
> <ruipbarra...@sapo.pt<mailto:ruipbarra...@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@r-project.org<mailto:R-help@r-project.org> mailing list
> >> https://stat.ethz.ch/mailman/**listinfo/r-
> help<https://stat.ethz.ch/mailman/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@r-project.org<mailto:R-help@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.
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
