Я пытаюсь придумать короткую программу, которая может вычислить тест отношения правдоподобия, начиная с logLikelihood двух вложенных моделей.
Затем мне нужно вычислить отношение правдоподобия, а затем для этого отношения,Я должен сравнить его со значениями хи-квадрат на моем значимом уровне (0,05)
Пока я написал эту программу, но кажется, что она не работает должным образом:
from scipy.stats.distributions import chi2
def likelihood_ratio(llmin, llmax):
return(2*(llmax-llmin))
LR = likelihood_ratio(-45523.66,-46458.58)
import scipy as sp
p = sp.stats.chi2.cdf(LR, 7)
print('p: %.30f' % p)
результатэтой программы 0,0000 Но в соответствии с моим примером книги это должно быть 0,05, а критическое значение должно быть 14,067 (взято из таблицы хи-квадрат с уровнем достоверности 0,05)
Как я могу получить критическое значение изтаблицу хи-квадрат на желаемом уровне достоверности и сравните ее с моим LR?