使用广义极值分布 (GEV) 计算返回值

Question

我有 30 年的年度温度数据，想要使用 50 年和 100 年的GEV 分布计算该数据的返回值返回周期。

我的30年数据：

data=[28.01,29.07,28.67,21.57,21.66,24.62,21.45,28.51,22.65,21.57,20.89,20.96,21.05,22.29,20.81,21.08,20.77,23.18,22.98,21.88,21.07,20.74,22.69,22.42,31.81,25.78,29.09,28.11,22.18,21.6]

如何使用GEV找到返回值？

Answer 1

要估计给定回报期 T 的回报水平，首先估计广义极值分布的参数，然后计算拟合分布的 1/T 处的生存函数的反函数。 （生存函数 SF(x) 只是 1 - CDF(x)。如果您阅读有关计算回报水平的内容，您通常会看到解决 CDF(x) = 1 - 1/T 的问题。这是相同的求解 SF(x) = 1/T。）

这是一个使用

scipy.stats.genextreme

来估计数据在多个返回周期的返回水平的脚本。方法

genextreme.isf()

是生存函数的逆函数。

import numpy as np
from scipy.stats import genextreme


data = np.array([28.01, 29.07, 28.67, 21.57, 21.66, 24.62, 21.45, 28.51,
                 22.65, 21.57, 20.89, 20.96, 21.05, 22.29, 20.81, 21.08,
                 20.77, 23.18, 22.98, 21.88, 21.07, 20.74, 22.69, 22.42,
                 31.81, 25.78, 29.09, 28.11, 22.18, 21.6])

# Fit the generalized extreme value distribution to the data.
shape, loc, scale = genextreme.fit(data)
print("Fit parameters:")
print(f"  shape: {shape:.4f}")
print(f"  loc:   {loc:.4f}")
print(f"  scale: {scale:.4f}")
print()

# Compute the return levels for several return periods.
return_periods = np.array([5, 10, 20, 50, 100])
return_levels = genextreme.isf(1/return_periods, shape, loc, scale)

print("Return levels:")
print()
print("Period    Level")
print("(years)   (temp)")

for period, level in zip(return_periods, return_levels):
    print(f'{period:4.0f}  {level:9.2f}')

输出：

Fit parameters:
  shape: -0.9609
  loc:   21.5205
  scale: 1.0533

Return levels:

Period    Level
(years)   (temp)
   5      25.06
  10      29.95
  20      39.45
  50      67.00
 100     111.53

Answer 2

希望你能收到我的消息。我可以应用上述步骤使用 GEV 分布计算极端降水的重现期吗？我将不胜感激您的回复。谢谢！