[R] ecdf

Downey, Patrick PDowney at urban.org
Mon Apr 19 15:04:10 CEST 2010


Hello,

I'd like to plot an empirical cumulative distribution function, except
instead of the fraction of values < x, I'd like the fraction of values > x.


I think this can be done using the ecdf function in {Hmisc}. I installed
the package and loaded it. However, when following the example given in the
documentation, I get an error:

x <- rnorm(100)
ecdf(x,what='1-F')
Error in ecdf(x, what = "1-F") : unused argument(s) (what = "1-F")

I believe that this is because R is attempting to access the ecdf function
in base R, which does not have the what option. Am I correct, and if so,
how can I change that?

Note: I also tried to "do it myself" without the {Hmisc} ecdf function, and
couldn't figure out a way. 

x2 <- 1-ecdf(x)

doesn't work, and neither does

x2 <- rep(0,times=100)
for(i in 1:100){
  x2[i] <- 1-ecdf(x)[i]
}

Both result in errors.

Thanks in advance for any suggestions you can offer.

-Mitch



More information about the R-help mailing list