获取与ggplot + stat_ecdf（）相关的数据

Question

我喜欢ggplot2包的 stat_ecdf（）功能部分，我觉得这对于探索数据系列非常有用。然而，这只是视觉的，我想知道它是否可行 - 如果是的话 - 如何获得相关的表格？

请查看以下可重现的示例

MartialArtsCheckInAppPurchaseKey

Answer 1

我们可以重新创建数据：

#Recreate ecdf data
dat_ecdf <- 
  data.frame(x=unique(iris$Sepal.Length),
             y=ecdf(iris$Sepal.Length)(unique(iris$Sepal.Length))*length(iris$Sepal.Length))
#rescale y to 0,1 range
dat_ecdf$y <- 
  scale(dat_ecdf$y,center=min(dat_ecdf$y),scale=diff(range(dat_ecdf$y)))

下面2个图表看起来应该是一样的：

#plot using new data
ggplot(dat_ecdf,aes(x,y)) +
  geom_step() +
  xlim(4,8)

#plot with built-in stat_ecdf
ggplot(iris, aes_string(x = "Sepal.Length")) +
  stat_ecdf() +
  xlim(4,8)

Answer 2

正如@krfurlong在this question中向我展示的那样，ggplot2中的layer_data函数可以为您提供所需的确切信息，而无需重新创建数据。

p <- ggplot(iris, aes_string(x = "Sepal.Length")) + stat_ecdf()
p.data <- layer_data(p)

p.data中的第一列“ y”包含ecdf值。 “ x”是绘图中x轴上的Sepal.Length值。