因为需要将gumbel分布拟合效果和其他几个分布进行比较,而fitdistrplus里原有的分布里没有gumbel分布,需要自定义。
R语言 fitdistrplus包 fitdist()函数中文帮助文档(中英文对照)从这里找到了机翻的fitdist()函数帮助文档,其中实例(4)就是自定义gumbel函数。由于原链接里只给了机翻,这里把实例(4)拿出来专门翻译注释一下。
(4) defining your own distribution functions, here for the Gumbel distribution #定义您自己的分布函数,这里以Gumbel分布示例
#dedicated to probability distributions #致力于概率分布
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-blog(-log§) #分布函数
fitgumbel <- fitdist(serving,“gumbel”,start=list(a=10,b=10)) #gumbel函数的fitdist类对象
summary(fitgumbel) #拟合分布的参数估计
plot(fitgumbel)