Rqchisq関数は、p 値と自由度を対応するカイ 2 乗値に変換します。同等の Python ライブラリはありますか?
私は何も見つけずにSciPyを見回しました。
5172 次
2 に答える
20
それはscipy.stats.chi2.ppf- パーセントポイント関数 (cdf の逆) です。たとえば、R では:
> qchisq(0.05,5)
[1] 1.145476
Python で:
In [8]: scipy.stats.chi2.ppf(0.05, 5)
Out[8]: 1.1454762260617695
于 2013-08-06T01:59:26.440 に答える
3
@VadimKhotilovich が彼の回答で指摘しているように、使用できますscipy.stats.chi2.ppf。の関数chdtriを使用することもできますがscipy.special、引数として 1-p を使用します。
R:
> qchisq(0.01, 7)
[1] 1.239042
> qchisq(0.05, 7)
[1] 2.16735
スキピー:
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
chdtrioverを使用する唯一の利点は、はるかに高速scipy.stats.chi2.ppfであることです。
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
于 2013-08-06T02:05:31.610 に答える