R qchisq
関数は、p値と自由度の数を対応するカイ二乗値に変換します。同等のPythonライブラリがありますか?Rのqchisq関数に相当するPythonはありますか?
私は何も見つけずにSciPyを見回しました。
R qchisq
関数は、p値と自由度の数を対応するカイ二乗値に変換します。同等のPythonライブラリがありますか?Rのqchisq関数に相当するPythonはありますか?
私は何も見つけずにSciPyを見回しました。
これはscipy.stats.chi2.ppf
- パーセントポイント関数(逆関数)です。 Rで 例えば、:Pythonで
> qchisq(0.05,5)
[1] 1.145476
:
In [8]: scipy.stats.chi2.ppf(0.05, 5)
Out[8]: 1.1454762260617695
@VadimKhotilovichが彼の答えで指摘しているように、あなたがscipy.stats.chi2.ppf
を使用することができます。 scipy.special
の関数chdtri
を使用することもできますが、1-pを引数として使用します。
R:
> qchisq(0.01, 7)
[1] 1.239042
> qchisq(0.05, 7)
[1] 2.16735
scipyのダウンロード:
In [16]: from scipy.special import chdtri
In [17]: chdtri(7, 1 - 0.01)
Out[17]: 1.2390423055679316
In [18]: chdtri(7, 1 - 0.05)
Out[18]: 2.1673499092980579
scipy.stats.chi2.ppf
上chdtri
を使用する唯一の利点は、それが速くくらいであるということである。
In [30]: from scipy.stats import chi2
In [31]: %timeit chi2.ppf(0.05, 7)
10000 loops, best of 3: 135 us per loop
In [32]: %timeit chdtri(7, 1 - 0.05)
100000 loops, best of 3: 3.67 us per loop
@mnel:その質問カイ2乗値からp値を計算することに関するue。私は他の方向に進もうとしています。確かに – jveldridge
。閉じる投票が取り消されました(もっと慎重に読む必要があります) – mnel
http://stackoverflow.com/questions/11725115/p-value-from-chi-sq-test-statistic-in-python –