我可以在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”

说这是不可能的，这是不言而喻的。