计算和生成 R 中超出概率的图？

Question

我正在尝试在 R 中计算 probability of exceedance。这里是 link，其中包含计算超出概率的详细信息和公式。我尝试使用以下示例代码复制该过程，但是，我离所附示例绘图很远。有一次，我想出了如何做到这一点——然后我想在一个有多个variables 的data.frame 上应用这个过程。

library(tidyverse)

A = sample(0:5000, 2500)
A = A[order(A, decreasing = TRUE)]
Rank = 1:2500

DF = data.frame(cbind(A,Rank))
DF$Prob =  100*(DF$Rank/(length(DF$Prob+1)))

ggplot(data = DF, aes(x=Prob, y=A))+
  geom_line() + scale_y_continuous(trans = "log10")

这是我想要制作的情节样本

Answer 1

计算Prob的公式在问题中没有很好地实现。它将 1 添加到 DF$Prob，然后将其 length 添加 1 到 length。

DF <- data.frame(A, Rank)
DF$Prob <-  DF$Rank/(length(DF$Rank) + 1)

ggplot(data = DF, aes(x = Prob, y = A)) +
  geom_line() + 
  scale_x_continuous(breaks = seq(0, 1, by = 0.20),
                     labels = percent) +
  scale_y_continuous(trans = "log10")

数据创建代码。

我已更改数据集示例。
下面的代码是可重现的，因为使用了set.seed，而sample 的概率会递减。

set.seed(1234)
A <- sample(0:5000, 2500, prob = exp(seq(10, 0, length.out = 5001)))
A <- A[order(A, decreasing = TRUE)]
Rank <- 1:2500