r - 如何从 Ecdf() 返回值中提取 Ecdf 值？

Question

This 的答案这里的问题提出了一种应用ecdf的方法。

但是，我使用 Hmisc 包中的 Ecdf()，因为它提供了一种方便的方法来绘制 ccdf(互补累积分布函数)图。 (通过将 what 选项设置为“1-F”)

默认情况下，Ecdf() 进行绘图并返回包含 x 和 y 的嵌套列表。

如何提取特定x值的y值？然后将其绘制在原始绘图上？

仅供引用:

 > str(Ecdf(rnorm(20), lwd = 2))
    List of 2
     $ x: num [1:21] -1.46 -1.46 -1.18 -1.17 -1.16 ...
     $ y: num [1:21] 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 ...
     - attr(*, "N")=List of 2
      ..$ n: num 20
      ..$ m: num 0

首先，我正在考虑使用 R List to Data Frame 建议的方法将此列表转换为 data.frame ，但我的数据很大，而且 rbind 看起来真的很慢。

Answer 1

Ecdf 返回一个列表，而 ecdf 返回一个函数。使用 R-stats 函数 ecdf 比使用类似以下的东西要容易得多:Ecdf(.)$y[ min(which(Ecdf(.)$x>val) ))]。如果您想要来自 Ecdf-object 的 ecdf(x=0) 值，那么这应该可行:

 ecdf( Ecdf(rnorm(20), lwd = 2)$x ) (v=0)
[1] 0.5238095

(事实证明，ecdf 返回的函数的形式参数是“v”。)但是，如果您想要不太优雅的方法，并且您已经将结果分配给名为 ' 的对象， oneEcdf':

oneEcdf <- Ecdf(rnorm(20), lwd = 2)
oneEcdf$y[ min( which(oneEcdf$x > 0 ))]
[1] 0.6