我可以在extRemes或R中的其他极值包中定义概率函数吗？

Question

我想知道是否可以在extRemes包（或其他EV包）中定义一个函数。例如，要在fitdistrplus中定义Gumbel，我们可以执行以下操作：

dgumbel <- function(x, a, b) 1/b*exp((a-x)/b)*exp(-exp((a-x)/b))
pgumbel <- function(q, a, b) exp(-exp((a-q)/b))
qgumbel <- function(p, a, b) a-b*log(-log(p))

但是当我尝试使用fevd（来自extRemes）时，它不起作用。

fit1 <- fevd(as.vector(x), method = "MLE", type="NewDistr")

谢谢

Answer 1

据我了解，您想为type参数指定一个非标准值。但是，相应的错误

# Error in match.arg(type) : 
#   'arg' should be one of “GEV”, “GP”, “PP”, “Gumbel”, “Exponential”

说这不可能是很不言自明的。