r - 如何使用 ECDF 绘制 "Reverse"累积频率图

Question

绘制以下累积 频率 图没有问题
像这样。

     library(Hmisc)
     pre.test <- rnorm(100,50,10)
     post.test <- rnorm(100,55,10)
     x <- c(pre.test, post.test)
     g <- c(rep('Pre',length(pre.test)),rep('Post',length(post.test)))
     Ecdf(x, group=g, what="f", xlab='Test Results', label.curves=list(keys=1:2))

但我想以值 > x 的“反向”累积频率的形式显示图形。
(即相当于 what="1-f"的东西)。

有没有办法做到这一点？

除了使用 Hmisc 之外，R 中的其他建议也非常受欢迎。

Answer 1

使用 Musa 建议:

pre.ecdf <- ecdf(pre.test)
post.ecdf <- ecdf(post.test)

r <- range(pre.test,post.test)
curve(1-pre.ecdf(x), from=r[1], to=r[2], col="red", xlim=r)
curve(1-post.ecdf(x), from=r[1], to=r[2], col="blue", add=TRUE)

您可以设置一些参数，如标题、图例等。

如果你想要频率而不是比例，简单的解决方案是:

pre.ecdf <- ecdf(pre.test)
post.ecdf <- ecdf(post.test)

rx <- range(pre.test,post.test)
ry <- max(length(pre.test),length(post.test))
curve(length(pre.test)*(1-pre.ecdf(x)), from=rx[1], to=rx[2], col="red", xlim=rx, ylim=c(0,ry))
curve(length(post.test)*(1-post.ecdf(x)), from=rx[1], to=rx[2], col="blue", add=TRUE)