[Rd] Quantile issue

Therneau, Terry M., Ph.D. therneau at mayo.edu
Wed Apr 30 19:00:33 CEST 2014


This is likely yet another instance of round off error, but it caught me by surprise.

tmt% R --vanilla
  (headers skipped, version 3.0.2 on Linux)

> load('qtest.rda')
> length(temp)
[1] 3622
> max(temp) >= quantile(temp, .98)
   98%
FALSE

I can send the file to anyone who would like to understand this more deeply.
The top 3% of the vector is a single repeated value.

Terry Therneau



More information about the R-devel mailing list