r - 使用带有直方图数据的 R 绘制 ECDF

Question

我有形式的直方图数据

Key  |  #occurences_of_key
--------------------------
 -10 | 1200
   0 | 1000
  10 | 700
  33 | 500
  67 | 200
  89 | 134
--------------------------

实现它的代码：

structure(c(-10, 0, 10, 33, 67, 89, 1200, 1000, 700, 500, 200, 134), .Dim = c(6L, 2L))

我想用这些数据绘制一个Empirical Cumulative Distribution Chart (percentile chart)使用。R我是新手R，所以我很感激任何指示。我阅读了有关ecdf可用功能的信息，R但我很难理解。

Answer 1

我能想到的一种方法是rep用来重建原始数据并ecdf在其上使用。

mat <- structure(c(-10, 0, 10, 33, 67, 89, 1200, 1000, 700, 500, 200, 134), .Dim = c(6L, 2L))

original <- unlist(apply(mat, 1, function(x) rep(x[1], x[2])))

original_ecdf <- ecdf(original)

plot(original_ecdf)

Answer 2

如果您的数据很大（这就是为什么您在加载到 R 之前预先将其制成表格的原因），您不想再次生成一些“虚拟”数据。您可以破解执行ecdf以接受表格数据：

tab_ecdf <- function (xs, counts) 
{
  n <- sum(counts)
  if (n < 1) 
    stop("'x' must have 1 or more non-missing values")
  rval <- approxfun(xs, cumsum(counts) / n, 
                    method = "constant", yleft = 0, yright = 1, f = 0, ties = "ordered")
  class(rval) <- c("ecdf", "stepfun", class(rval))
  assign("nobs", n, envir = environment(rval))
  attr(rval, "call") <- sys.call()
  rval
}

然后用它代替原来的ecdf()功能。