如何找到R中GEV分布中给定值的累积概率？

我已经将我的数据拟合到GEV分布，我想知道如何根据fevd的帮助页面，部分

详细信息
，找到p（x的概率：

GEV df由以下公式给出：

PrX
library（EnvStats）
ams非常感谢您的帮助。但是，当我检查
pevd
的参数时，我发现它需要
q
，这是一个分位数的数字向量。因此，我们可以在这里使用
40
（因为40是一个值，而不是分位数）？@YangYang是的，在文档中它调用了参数
q
，但我已经从上面的代码中得到了一个
x
，仅此而已。
40
是一个分位数。很抱歉，我有点困惑。在我的例子中，40是一个数据点，而不是分位数。例如，40代表40毫米降雨量。@YangYang，这是你想要的名称问题你想要
P（X好的，谢谢你的解释。我在
https://stackoverflow.com/questions/48434236/how-to-plot-a-fevd-object-using-ggplot2-in-r
，不知您是否可以看一看？谢谢您抽出时间。
library(extRemes)
ams <- c(44.5,43.2,38.1,39.1,32.3,25.4,33.0,32.5,48.5,34.3,45.7,35.3,76.7,34.0,86.6,48.5,59.4,53.3,30.5,42.7,83.3,59.2,37.3,67.3,38.4,47.0,38.1,72.4,40.9,47.0,36.3,85.3,35.6,55.9,44.2,45.2,51.6,59.4,47.8,55.4,42.4,40.1,36.6,47.0,48.8,51.3,39.4,45.7)
fit_mle <- fevd(x=ams, method = "MLE", type="GEV",period.basis = "year")

location <- fit_mle$results$par[1]
scale <- fit_mle$results$par[2]
shape <- fit_mle$results$par[3]
x <- 40
exp(-(1 + shape*(x - location)/scale)^(-1/shape))
#    shape 
#0.3381735

pevd(x, location, scale, shape)
#[1] 0.3381735

library(EnvStats)
ams <- c(44.5,43.2,38.1,39.1,32.3,25.4,33.0,32.5,48.5,34.3,45.7,35.3,76.7,34.0,86.6,48.5,59.4,53.3,30.5,42.7,83.3,59.2,37.3,67.3,38.4,47.0,38.1,72.4,40.9,47.0,36.3,85.3,35.6,55.9,44.2,45.2,51.6,59.4,47.8,55.4,42.4,40.1,36.6,47.0,48.8,51.3,39.4,45.7)
fit_gev <- egevd(ams, method = "mle")# Parameters estimation
pgevd(40, location = fit_gev$parameters[[1]], scale = fit_gev$parameters[[2]],
  shape = fit_gev$parameters[[3]])

0.3381751