GEV 不合适

Question

我有每年的最大值，我尝试对其进行 GEV 分析。然而，当尝试时，pdf 根本不适合数据。我尝试最小化类似问题中提出的建议，但这效果不佳。代码如下：

%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import genextreme

yearly_max = [3465, 3461, 3318, 3304, 3211, 3307, 3224, 3440, 3317, 3393, 3360, 3695,
            3265, 3328, 3326, 3272, 3479, 3275, 3355, 3316, 3324, 3411, 3298, 3424,
            3311, 3387, 3374, 3320, 3307, 3343, 3440]

shape_gev, location_gev, scale_gev = genextreme.fit(yearly_max)

input_WL = np.linspace(2000, 6000, 1000)
pdf_gev = genextreme.pdf(input_WL, shape_gev, loc=location_gev, scale=scale_gev)
plt.hist(yearly_max, bins=40, range=(3000, 4000), density = True, label="normed occurence")
plt.plot(input_WL, pdf_gev, linewidth=4, label='pdf')

有人可以帮我解决这个问题吗？