假设尺度参数是时间的函数,使用 Cunnane Plotting Position 和 Cramer von-Mises 拟合优度标准拟合 Gumbel 分布。
我的数据格式如下:
年 | X |
---|---|
1945年 | 10700 |
1946年 | 5260 |
1947年 | 424 |
1948年 | 10724 |
1949年 | 22100 |
1950 | 2520 |
1951年 | 589 |
1952年 | 3320 |
1953年 | 4970 |
1954年 | 5690 |
#Random values
X <- c(10700,5260,424,161,8435,1656,6516,6863,1561,1681,1556,651,6849,1533)
#Cunnane (alpha = 0.4)
alpha <- 0.4
cvm.test(X, null = "pgumbel", alpha=alpha, scale)
我感到困惑的是如何使比例参数成为时间的函数?