使用 ggplot2 在同一图上使用 ECDF

Question

我有一个数据框，在应用融化函数后，它看起来类似于：

 var       val
1 a 0.6133426
2 a 0.9736237
3 b 0.6201497
4 b 0.3482745
5 c 0.3693730
6 c 0.3564962

......

初始数据框有 3 列，列名称为 a、b、c以及它们的相关值。 我需要在同一张图上绘制，使用 ggplot 为每一列关联 ecdf (ecdf(a)、ecdf(b)、ecdf(c))，但我未能做到这一点。我尝试过：

p<-ggplot(melt_exp,aes(melt_exp$val,ecdf,colour=melt_exp$var))
pg<-p+geom_step()

但我收到错误：参数意味着不同的行数：34415, 0。

有谁知道如何做到这一点？该图应类似于plot(ecdf(x)) 返回的图，而不是阶梯状的图。

谢谢你！

Answer 1

我的第一个想法是尝试使用 stat_function ，但由于 ecdf 返回一个函数，我无法快速使其工作。相反，这里有一个解决方案，要求您首先将计算值附加到数据框（使用 Ramnath 的示例数据）：

library(plyr) # function ddply()
mydf_m <- ddply(mydf_m, .(variable), transform, ecd = ecdf(value)(value))

ggplot(mydf_m,aes(x = value, y = ecd)) + 
    geom_line(aes(group = variable, colour = variable))

如果您想要平滑地估计 ECDF，您还可以将 geom_smooth 与 spline 包中的函数 ns() 一起使用

library(splines) # function ns()
ggplot(mydf_m, aes(x = value, y = ecd, group = variable, colour = variable)) + 
    geom_smooth(se = FALSE, formula = y ~ ns(x, 3), method = "lm")

<图片src="https://i.sstatic.net/WWbuv.png" alt="在此处输入图像描述">

正如上面评论中所述，从版本 0.9.2.1 开始，ggplot2 有一个用于此目的的特定统计数据： < a href="http://docs.ggplot2.org/0.9.2.1/stat_ecdf.html" rel="nofollow noreferrer">stat_ecdf。使用它，我们只需做这样的事情：

ggplot(mydf_m,aes(x = value)) + stat_ecdf(aes(colour = variable))

Answer 2

基于 Ramnath，采用上述方法，您可以通过执行以下操作从 ggplot2 获取 ecdf：

require(ggplot2)
mydf = data.frame(
   a = rnorm(100, 0, 1),
   b = rnorm(100, 2, 1),
   c = rnorm(100, -2, 0.5)
)

mydf_m = melt(mydf)

p0 = ggplot(mydf_m, aes(x = value)) + 
   stat_ecdf(aes(group = variable, colour = variable)) 
print(p0)

Answer 3

这是一种方法

require(ggplot2)
mydf = data.frame(
  a = rnorm(100, 0, 1),
  b = rnorm(100, 2, 1),
  c = rnorm(100, -2, 0.5)
)

mydf_m = melt(mydf)

p0 = ggplot(mydf_m, aes(x = value)) + 
  geom_density(aes(group = variable, colour = variable)) +
  opts(legend.position = c(0.85, 0.85))