numpy.random.chisquare#
- random.chisquare(df, size=None)#
從卡方分佈中抽取樣本。
當 df 個獨立的隨機變數(每個都具有標準常態分佈,平均值為 0,變異數為 1)被平方並加總時,結果分佈為卡方分佈(請參閱「註解」)。此分佈常用於假設檢定。
- 參數:
- dffloat 或 float 的類陣列 (array_like)
自由度,必須 > 0。
- size整數 (int) 或整數元組 (tuple of ints),選填
輸出形狀。如果給定的形狀為,例如
(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