广告合作
  • 今日头条

    今日头条

  • 百度一下

    百度一下,你就知道

  • 新浪网

    新浪网 - 提供新闻线索,重大新闻爆料

  • 搜狐

    搜狐

  • 豆瓣

    豆瓣

  • 百度贴吧

    百度贴吧——全球领先的中文社区

  • 首页 尚未审核订阅工具 订阅

    如何用python进行似然比检验

    来源:网络收集  点击:  时间:2024-03-04
    【导读】:
    统计学中,似然比检验城用来对同一问题的两种模型的拟合优度的检验,这篇经验将为大家介绍如何利用python进行似然比检验。工具/原料morepython方法/步骤1/6分步阅读

    在似然比检验中,我们设参数空间为θ我们的原假设与备选假设为:

    原假设:参数θ属于参数空间的θ的一个子集θ0

    备选假设:参数θ属于θ0在参数空间的θ下的补集

    如果用数学语言,可以表示为:

    2/6

    似然函数表示的是当前参数条件下,出现待估计数据的可能性大小,可以看做是x确定条件下的参数θ的一个函数,通常用x的概率密度函数或者概率质量函数来表示,数学公式表示为:

    3/6

    似然比检验需要使用两个模型的似然函数来构建统计量,其统计量公式可表示如下,其中sup表示集合的上确界:

    4/6

    接下来,我们要确定的是该检验的拒绝域(reject region)或者关键值(critical value),常用的办法是,根据θ0和θ的参数个数差,来获得对应的卡方分布的分位数,卡方分布的自由度为参数个数差,分位数的比例为统计检验的置信度(通常为0.95),假设自由度为df,则对应的数值为:

    5/6

    接下来,我们就可以利用python来获得上述的关键值了,scipy这一包为我们提供了各类常用分布的各类函数,其中ppf函数为我们计算对应百分比的分位数,如果要是用卡方分布的功能,需要通过:

    from scipy.stats import chi2

    来引入chi2

    然后使用:

    chi2.ppf(q, df)

    来获取对应的卡方分布分位数。

    图示为自由度为21的卡方分布95%分位数:

    6/6

    接下来我们将-2log(X)来于对应的关键值c(刚才获得的分位数)进行比较,如果-2log(X)大于关键值c,则不拒绝原假设,反之拒绝原假设。

    注意事项

    如果大家觉得有帮助,随手点赞或投票哦~~

    本文关键词:

    版权声明:

    1、本文系转载,版权归原作者所有,旨在传递信息,不代表看本站的观点和立场。

    2、本站仅提供信息发布平台,不承担相关法律责任。

    3、若侵犯您的版权或隐私,请联系本站管理员删除。

    4、文章链接:http://www.ff371.cn/art_238210.html

    相关资讯

    ©2019-2020 http://www.ff371.cn/ 国ICP备20009186号05-13 11:15:30  耗时:0.022
    0.0225s