Chi2cdf верхний питон

Question

Итак, я перевожу Matlab в код Python и наткнулся на эту строку

Pv = chi2cdf(Devk,Md,'upper');

Где Devk - трехмерная матрица, а Md - df = 10.

Я искалэквивалент в Python.Я нашел это:

from scipy.stats import chi2
Pv = chi2.cdf(Devk,Md)

Однако, это не дает аналогичных результатов, я думаю, из-за «верхнего» аргумента, который имеет Matlab.У кого-нибудь есть идеи, как реализовать это в Python?

Кто-то попросил сделать некоторые игрушечные данные, чтобы показать проблему:

Matlab:

Devk = [1,2,3,4,5]
Md = 10
Pv = chi2cdf(Devk,Md,'upper')
Pv =
0.9998    0.9963    0.9814    0.9473    0.8912

Python:

from scipy.stats import chi2
Devk = np.array([1,2,3,4,5])
Md = 10
Pv = chi2.cdf(Devk,Md)
Pv
Out[125]: array([0.00017212, 0.00365985, 0.01857594, 0.05265302, 0.10882198])

Answer 1

Я считаю, что справедливо следующее:

chi2cdf(..., 'upper') == (1 .- chi2cdf(...))

Итак, в вашем коде Python, скорее всего, вам нужно сделать:

from scipy.stats import chi2

Pv = 1 - chi2.cdf(Devk, Md)

РЕДАКТИРОВАТЬ:

Еще лучше использовать:

Pv = chi2.sf(Devk, Md)

, потому что, согласно документации 1013 *, функция выживания sf() (также определенная как 1 - cdf()) иногда более точна.