我试图在python中使用内置的最大似然法而不是在MATLAB中获得GEV函数的参数估计,但MATLAB(gevfit_rev)中估计的参数与使用scipy.stats.genextreme.fit找到的参数不同。 。使用数据的示例
data = np.array([ 22.20379411, 22.99151292, 24.27032696, 24.82180626,
25.23163221, 25.39987272, 25.54514567, 28.56710007,
29.7575898 , 30.15641696, 30.79168255, 30.88147532,
31.0236419 , 31.17380647, 31.61932755, 32.23452568,
32.76262978, 33.39430032, 33.81080069, 33.90625861,
33.99142006, 35.45748368, 37.0342621 , 37.14768791,
38.14350221, 42.72699534, 44.16449992, 48.77736737,
49.80441736, 50.57488779])
Scipy提供参数(形状,位置,比例)
genextreme.fit(data)
(0.46794425159397518, 29.321600290365303, 4.7986571014800141)
MATLAB提供(形状,比例,位置)
gevfit_rev(data)
0.0144 5.8983 29.7621
有人能说出为什么这些估计有所不同吗?