r - 如何在对数刻度上绘制 CCDF 图？

Question

我想在对数轴上为我的一些模拟幂律尾部数据绘制 CCDF 图，下面是我在法线轴上绘制 CCDF 图的 R 代码，我使用了链接上的代码: ( How to plot a CCDF gragh? )

> load("fakedata500.Rda")
> x<-fakedata500
> f<-ecdf(x)
> f
Empirical CDF 
Call: ecdf(x)
 x[1:500] = 0.50174, 0.50307, 0.50383,  ..., 81.674, 140.63
> plot(f)

下面是 ECDF 图:

> plot(sort(x), 1-f(sort(x)), type="s", lwd=1)

此命令为我提供了 CCDF 图:

但是，我想在对数轴上绘制我的 CCDF 图以生成如下图所示的图形:(图表来自“最小化识别生物体的 Lévy 飞行行为的错误。”)

有没有办法在 R 中做到这一点？

如果是这样，如何对 CCDF 图进行线性回归？我试过使用下面的命令，但这对我不起作用。

a<-plot(sort(x), 1-f(sort(x)), type="s", lwd=1)
> a
NULL
> res=lm(a)
Error in terms.formula(formula, data = data) : 
  argument is not a valid model

非常感谢。

---

更新:

我使用了@BondedDust给出的代码，成功生成了一张CCDF图:

(plot(sort(x) , 1-ecdf(x)(sort(x) ), log="xy"))

下面是我如何生成数据集的代码:

u<-runif(500)
fakedata500<-((2*(1-u))^(-1))

Answer 1

您没有提供任何数据，因此不会尝试回归输出，但构建绘图应该非常容易:

plot(sort(x) , 1-ecdf(x)(sort(x) ), log="xy")

您确实会收到警告，因为最极端点的 ecdf 为 1，而 1-1 不能有对数值，因为它被省略了。使用绘图看起来像的对数正态变量的随机样本。轴的替代标记可以通过在绘图调用中抑制轴刻度然后使用 axis 函数添加它们来完成。这是通过使用默认参数的对数正态分布的 500 个实现样本完成的: