- 积分
- 766
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-9-23
- 最后登录
- 1970-1-1
|
发表于 2021-6-25 21:38:35
|
显示全部楼层
# (4) defining your own distribution functions, here for the Gumbel distribution
# for other distributions, see the CRAN task view
# 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-b*log(-log(p))
fitgumbel <- fitdist(serving, "gumbel", start=list(a=10, b=10))
summary(fitgumbel)
plot(fitgumbel)
来自R语言fitdist函数(fitdistrplus包)帮助文档,希望对你有所启发。 |
|