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