1

我现在正在使用该extremes包来拟合广义极值 (GEV) 分布，并且我想使用 Kolmogorov-Smirnov 检验来估计拟合优度，但出现以下错误：

library(extRemes)
library(eva)
data("PORTw", package = "extRemes")
fit1 <- fevd(TMX1, PORTw, units = "deg C")
ks.test(PORTw$TMX1,"pgev",fit1$results$par[[1]],fit1$results$par[[2]],shape=fit1$results$par[[3]])

`Warning message:
In ks.test(PORTw$TMX1, "pgev", fit1$results$par[[1]], fit1$results$par[[2]],  :
  ties should not be present for the Kolmogorov-Smirnov test`

所以，我的问题是，如何对 GEV 进行 Kolmogorov-Smirnov 检验以符合关系？或者，是否有任何其他拟合优度测试可用于拟合 R 中可用的分布？非常感谢。

4

1 回答 1

1

我推荐“EnvStats”包。您将有更多的适合度测试的多功能性：

library(EnvStats)
# For a data set called X
X <- rgevd(500)
# Generalized Extreme Value (EnvStats)
egevd(X, method = "mle")# Maximum likelihood
# Goodness of fit test
gofTest(X, distribution = "gev",test = "ks")#Kolmogorov-Smirnov
gofTest(X, distribution = "gev",test = "chisq")#Chi-Squared
于 2019-04-30T17:36:48.440 回答