[R] ecdf

David Winsemius dwinsemius at comcast.net
Mon Apr 19 16:24:33 CEST 2010


On Apr 19, 2010, at 9:12 AM, Downey, Patrick wrote:

> Hi Thierry,
>
> That worked perfectly. Thanks for the suggestion.
>
> For reference, in the documentation, it

What is "it"?

> never lists {Hmisc}'s function as
> starting with E instead of e.

Every instance of the documentation in the r help system using ??ecdf  
that I can find for Hmisc's version has it properly capitalized. I  
have Hmisc_3.7-0

> I don't know who's in charge of
> documentation

It would be the package maintainer if there were a problem.

> , but that

What is "that"?

-- 
David

> should probably be corrected.
>
> Thanks again.
>
> -Mitch
>
> -----Original Message-----
> From: ONKELINX, Thierry [mailto:Thierry.ONKELINX at inbo.be]
> Sent: Monday, April 19, 2010 9:08 AM
> To: Downey, Patrick; R help
> Subject: RE: [R] ecdf
>
> R is case sensitive. ecdf() is in the stats package, Ecdf() is in  
> Hmisc.
> So you want Ecdf(x,what='1-F')
>
> Thierry
> ------------------------------------------------------------------------
> ----
> ir. Thierry Onkelinx
> Instituut voor natuur- en bosonderzoek
> team Biometrie & Kwaliteitszorg
> Gaverstraat 4
> 9500 Geraardsbergen
> Belgium
>
> Research Institute for Nature and Forest
> team Biometrics & Quality Assurance
> Gaverstraat 4
> 9500 Geraardsbergen
> Belgium
>
> tel. + 32 54/436 185
> Thierry.Onkelinx at inbo.be
> www.inbo.be
>
> To call in the statistician after the experiment is done may be no  
> more
> than asking him to perform a post-mortem examination: he may be able  
> to
> say what the experiment died of.
> ~ Sir Ronald Aylmer Fisher
>
> The plural of anecdote is not data.
> ~ Roger Brinner
>
> The combination of some data and an aching desire for an answer does  
> not
> ensure that a reasonable answer can be extracted from a given body of
> data.
> ~ John Tukey
>
>
>> -----Oorspronkelijk bericht-----
>> Van: r-help-bounces at r-project.org
>> [mailto:r-help-bounces at r-project.org] Namens Downey, Patrick
>> Verzonden: maandag 19 april 2010 15:04
>> Aan: R help
>> Onderwerp: [R] ecdf
>>
>> 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
>>
>> ______________________________________________
>> 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.
>>
>
> Druk dit bericht a.u.b. niet onnodig af.
> Please do not print this message unnecessarily.
>
> Dit bericht en eventuele bijlagen geven enkel de visie van de  
> schrijver
> weer
> en binden het INBO onder geen enkel beding, zolang dit bericht niet
> bevestigd is
> door een geldig ondertekend document. The views expressed in  this  
> message
> and any annex are purely those of the writer and may not be regarded  
> as
> stating
> an official position of INBO, as long as the message is not  
> confirmed by a
> duly
> signed document.
>
> ______________________________________________
> 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.

David Winsemius, MD
West Hartford, CT



More information about the R-help mailing list