numpy.random.RandomState.chisquare#
方法
- random.RandomState.chisquare(df, size=None)#
從卡方分佈中抽取樣本。
當 df 個獨立隨機變數(每個都具有標準常態分佈,平均值為 0,變異數為 1)平方並加總時,結果分佈為卡方分佈(請參閱「註解」)。此分佈常用於假設檢定。
- 參數:
- dffloat 或 float 的 array_like
自由度數量,必須 > 0。
- sizeint 或 int 元組,選用
輸出形狀。如果給定的形狀為例如
(m, n, k)
,則會抽取m * n * k
個樣本。如果 size 為None
(預設),如果df
是純量,則會傳回單一值。否則,會抽取np.array(df).size
個樣本。
- 傳回:
- outndarray 或純量
從參數化的卡方分佈中抽取的樣本。
- 引發:
- ValueError
當 df <= 0 或給定不適當的
size
(例如size=-1
) 時。
另請參閱
random.Generator.chisquare
新程式碼應使用此方法。
註解
透過加總 df 個獨立的、標準常態分佈隨機變數的平方所獲得的變數
\[Q = \sum_{i=1}^{\mathtt{df}} X^2_i\]是卡方分佈,表示為
\[Q \sim \chi^2_k.\]卡方分佈的機率密度函數為
\[p(x) = \frac{(1/2)^{k/2}}{\Gamma(k/2)} x^{k/2 - 1} e^{-x/2},\]其中 \(\Gamma\) 是伽瑪函數,
\[\Gamma(x) = \int_0^{-\infty} t^{x - 1} e^{-t} dt.\]參考文獻
[1]NIST “工程統計手冊” https://www.itl.nist.gov/div898/handbook/eda/section3/eda3666.htm
範例
>>> np.random.chisquare(2,4) array([ 1.89920014, 9.00867716, 3.13710533, 5.62318272]) # random