如何在scipy genextreme分佈中設置「尾部概率」？

有沒有人知道如何在scipy的'genextreme'分佈中設置'q'參數（控制較低或較高的尾部概率）？如何在scipy genextreme分佈中設置「尾部概率」？

＃的/ usr/bin中/ env的蟒

進口numpy的作爲NP
進口pylab作爲PLT
從scipy.stats導入genextreme

畝，標準差= 10，4
＃問題：如何在genextreme中設置q，較低或較高的尾部概率？
RV = genextreme（C = [0]，在上述=畝，標度=西格馬）
X = np.linspace（0，40，50）＃生成從0-40

y_pdf = RV 50個值。 PDF（X）
plt.plot（X，y_pdf， 'RO-'，標記= 'PDF'）

y_cdf = rv.cdf（X）
plt.plot（X，y_cdf，「b * - 」，標記= 'CDF'）

plt.legend（）
plt.show（）

@Warren：謝謝你提供的信息。我是新來的scipy統計。我怎樣才能移動分佈，使得pdf在高端更加重要？我嘗試將'c'參數改爲[0.9]，但輸出非常不穩定，隨着曲棍棒的上升，然後劇烈的幾乎直線下降到零。

來源

2014-01-27 TrailDreaming

我在回答中添加了更多信息。 –

q不是scipy的參數genextreme distribution。它是方法ppf（cdf的倒數）和isf（存活函數的倒數）的參數。如果您查看其他發行版的文檔 - 例如。 norm,gamma - 你會發現所有的班級文件列表在其「參數」中列出了q，但該文件是對所有方法的參數的概述。 genextreme具有標準的位置和比例參數，以及一個形狀參數c。

例子：

>>> genextreme.cdf(genextreme.ppf(0.95, c=0.5), c=0.5) 
0.94999999999999996

你可以找到更多關於generalized extreme distribution on wikipedia，但要注意，在使用SciPy的形狀參數的符號相反的維基百科文章中的形狀參數。例如，下面的生成維基百科情節：

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

# Create the wikipedia plot. 
plt.figure(1) 
x = np.linspace(-4, 4, 201) 
plt.plot(x, genextreme.pdf(x, 0.5), color='#00FF00', label='c = 0.5') 
plt.plot(x, genextreme.pdf(x, 0.0), 'r', label='c = 0') 
plt.plot(x, genextreme.pdf(x, -0.5), 'b', label='c = -0.5') 
plt.ylim(-0.01, 0.51) 
plt.legend(loc='upper left') 
plt.xlabel('x') 
plt.ylabel('Density') 
plt.title('Generalized extreme value densities') 
plt.show()

結果：

genextreme densities

c時0 <，支撐分佈的左端是loc + scale/c。下面創建與c = -0.25，和loc = -scale/c分佈的曲線圖（這樣所述支撐件的左端爲0）：

c = -0.25 
scale = 5 
loc = -scale/c 
x = np.linspace(0, 80, 401) 
pdf = genextreme.pdf(x, c, loc=loc, scale=scale) 
plt.plot(x, pdf) 
plt.xlabel('x') 
plt.ylabel('Density')

簡介：

genextreme example

的上stats tutorial scipy文檔站點有關於分發類的更多信息，包括關於shifting and scaling a distribution的部分。

來源

2014-01-27 03:50:26

如何在scipy genextreme分佈中設置「尾部概率」？

回答

相關問題