如何在R中拟合Gumbel分布？含ismev、extRemes包参数固定方法问询

一、使用已有的GEV包实现 #

1. ismev包

ismev的gev.fit()函数支持固定特定参数，只需通过fixed参数将形状参数xi固定为0，同时传入位置参数mu和尺度参数sigma的初始值即可：

library(ismev)
# 替换为你的数据向量
data_vec <- rgev(100, loc = 5, scale = 2, shape = 0)
# 为mu和sigma设置初始值（通常用样本均值和标准差）
init_params <- c(mean(data_vec), sd(data_vec))
# 拟合Gumbel分布（固定xi=0）
gumbel_fit <- gev.fit(data_vec, init = init_params, fixed = list(xi = 0))
# 查看拟合结果
summary(gumbel_fit)

2. extRemes包

extRemes提供了两种更直观的方式拟合Gumbel分布：

• 方式一：直接指定分布类型为"gumbel"

library(extRemes)
data_vec <- rgev(100, loc = 5, scale = 2, shape = 0)
# 直接拟合Gumbel分布
gumbel_fit1 <- fevd(data_vec, type = "gumbel")
summary(gumbel_fit1)

• 方式二：在GEV模型中固定xi=0

# 通过fixed参数固定形状参数为0
gumbel_fit2 <- fevd(data_vec, type = "gev", fixed = list(xi = 0))
summary(gumbel_fit2)

两种方式的拟合结果完全一致，第一种更简洁。

二、其他可用于Gumbel分布推断的R包 #

1. fitdistrplus包

这是一个专注于分布拟合的工具包，支持极大似然、矩估计等多种方法，直接指定distr = "gumbel"即可拟合：

library(fitdistrplus)
data_vec <- rgev(100, loc = 5, scale = 2, shape = 0)
gumbel_fit <- fitdist(data_vec, distr = "gumbel", method = "mle")
summary(gumbel_fit)

该包还提供了拟合诊断图（如QQ图、密度拟合图），方便验证拟合效果。

2. MASS包（基础生态）

MASS包中的fitdistr()函数支持Gumbel分布的极大似然估计，无需额外安装基础相关包：

library(MASS)
data_vec <- rgev(100, loc = 5, scale = 2, shape = 0)
gumbel_fit <- fitdistr(data_vec, "gumbel")
print(gumbel_fit)

注意：不同包对Gumbel分布的参数化可能存在细微差异（比如部分包可能使用尺度参数的倒数），拟合前建议查看对应函数的帮助文档确认参数定义。

内容的提问来源于stack exchange，提问作者oliver