使用evd包拟合形状参数>0的GEV模型报错求助

报错原因 #

报错'x' must be a non-empty numeric vector是因为fgev()函数要求输入数值向量，但你传入的是整个数据框df，而非数据框中的单列数值数据。

修正步骤及代码 #

1. 检查数据结构，确认目标列名称
   读取数据后先查看数据结构，明确你要拟合的列名（示例假设列名为W）：

   library(evd)  # R区分大小写，需用小写library()
   df <- read.csv("W.csv", header = TRUE)
   str(df)  # 查看数据结构，确认目标列的准确名称
   

2. 提取单列向量并拟合约束形状参数的GEV模型
   从数据框中提取目标列作为数值向量，同时通过lower参数约束形状参数大于0：

   # 提取目标列的数值向量，替换W为你的实际列名
   x <- df$W
   # 拟合形状参数>0的GEV模型，lower参数的第三个值设为0，限制形状参数下限
   fit <- fgev(x, lower = c(-Inf, -Inf, 0))
   

3. 获取所需结果
   查看参数估计值、标准误及对数似然值：

   # 输出参数估计值和标准误
   summary(fit)
   # 提取模型对数似然值
   logLik(fit)
   

额外注意事项 #

• 提前检查目标列数据：用is.numeric(x)确认列类型为数值型，sum(is.na(x))查看是否有缺失值，必要时用na.omit(x)清理缺失值。

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