应用广义极值(GEV)理论后,我有一个形状,位置和比例参数来描述我的分布。现在我试图在Python中使用这三个参数绘制CDF。在matlab中有一个"cdf" function可以做到这一点。我找不到用scipy做的方法吗?
答案 0 :(得分:2)
使用scipy.stats.genextreme
的cdf
方法。
例如,以下ipython会话......
In [40]: from scipy.stats import genextreme
In [41]: shape = -0.5
In [42]: loc = 0.0
In [43]: scale = 2.5
In [44]: x = np.linspace(scale/shape, 20, 200)
In [45]: y = genextreme.cdf(x, shape, loc, scale)
In [46]: import matplotlib.pyplot as plt
In [47]: plt.plot(x, y, label="GEV CDF")
Out[47]: [<matplotlib.lines.Line2D at 0x10e135790>]
In [48]: plt.legend(loc='best')
Out[48]: <matplotlib.legend.Legend at 0x10de4cc50>
生成此图:
请注意,scipy代码中的形状参数c
与wikipedia article on the generalized extreme value distribution中使用的形状参数ξ符号相反。