使用Chisqaure表计算似然比检验

时间:2019-06-05 10:27:04

标签: scipy

我正试图提出一个简短的程序,该程序可以从两个嵌套模型的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进行比较?

0 个答案:

没有答案