绘制步行持续时间的 CCDF

Question

我已经按照maximum plot points in R? 帖子中提到的问题部分绘制了 CCDF，以获得带有此代码的 plot(image1)：

ccdf<-function(duration,density=FALSE)
{
freqs = table(duration)
  X = rev(as.numeric(names(freqs)))
  Y =cumsum(rev(as.list(freqs)));
  data.frame(x=X,count=Y)
}
qplot(x,count,data=ccdf(duration),log='xy')

现在，根据teucer 对Howto Plot “Reverse” Cumulative Frequency Graph With ECDF 的回答，我尝试使用以下命令绘制CCDF：

f <- ecdf(duration)
plot(1-f(duration),duration)

我得到了一个类似 image2 的情节。
此外，我在Plotting CDF of a dataset in R? 中的一个答案中从 cmets 中读到，因为 CCDF 只不过是 1-ECDF。
我对如何获取数据的 CCDF 感到非常困惑。

图片1

Image2

Answer 1

我使用 ggplot 来获得我想要的数据的 ccdf 图，如下所示：

>>ecdf_x <- ecdf(x) 
>>dfr <- data.frame( ecdf = ecdf_x(x), 
>>ccdf = 1 - ecdf_x(x) ) 
>>p_ccdf <- ggplot(dfr, aes(x, ccdf)) + geom_line() + scale_x_log10() 
>>p_ccdf

抱歉这么晚才发布。 谢谢大家！

Answer 2

生成一些数据并找到 ecdf 函数。

x <- rlnorm(1e5, 5)
ecdf_x <- ecdf(x)

在 x 的范围内定期生成向量。 （编辑：在这种情况下，您希望它们在对数刻度上均匀分布；如果您有负值，则在线性刻度上使用样本。）

xx <- seq(min(x), max(x), length.out = 1e4)
#or
log_x <- log(x)
xx <- exp(seq(min(log_x), max(log_x), length.out = 1e3))

为绘图创建带有 x 和 y 坐标的数据。

dfr <- data.frame(
  x = xx,
  ecdf = ecdf_x(xx),
  ccdf = 1 - ecdf_x(xx)
)

画图。

p_ccdf <- ggplot(dfr, aes(x, ccdf)) + 
  geom_line() +
  scale_x_log10()
p_ccdf

（也可以看看aes(x, ecdf)。）