因为你这里主要想用到Fisher精确检验,我就只以Fisher使用来说明。

说到Fisher精确检验,不得不说“女士品茶”的经典案例。

奶和茶混合的杯子里,到底是先倒的茶还是先倒的奶呢?尝起来有啥不同呢?

假设现在让一位女士来品茶,我们记录下实际的倒茶顺序和这位女士所说的顺序。

其中n是制作的茶的总杯数。先倒茶的杯数是a+c,这位女士将其中的a+b分类为先倒茶。理论情况下,如果她能够尝到区别,b和c的计数应该很小。相反,如果她实际上不能分辨,那么a和c大致上会相等。

如果现在为了检验这位女士的能力,我们准备8杯茶,4杯先倒茶,4杯先倒奶,并告知这位女士其中4杯是先倒茶的,4杯是先倒奶的。杯子以随机的顺序递给她。她的目标就是分辨出哪4杯先倒茶,哪4杯先倒奶。

这个设计中,行和列的总和固定为4,即

a+b=a+c=c+d=b+d=4

在这些限制下,当a、b、c、d中任何一个确定下来,剩下的3个也被确定了:

b=4-a,c=4-a,d=a

如果我们认为她没有辨别能力,那么这种将顺序随机化的方法将让她选择为先倒茶的4杯的可能性等于从8杯中任意选4杯的可能性。

从8杯中选择4杯为“先倒茶”一共有(8选4)=70种可能的方式。(Python可以用comb实现)

结果:

如果假定她没有辨别能力,正确分类为先倒茶(表格中的a)的数量服从概率分布(Python可以用hypergeom实现)。对a来说,可以分为0、1、2、3、4五种可能性。

用Python来实现如图:

结果:

有了这些概率,我们就能计算检验该女士无法分辨这两种准备的假设的p值了。p值是在无效假设下,观察到一个比目前观察到的结果同样极端或更极端的概率。

如果她每次都说对了,p值就是1/70,如果她至少说对3次,p值就是1/70+16/70~0.24 。

上述检验成为“Fisher 精确检验”,Python用stats.fisher_exact()实现:

oddsratio,p = stats.fisher_exact(obs, alternative='greater')

其中obs是包含观测值的矩阵。

这里假设女士只猜对1次,则:

那么可以这样写:

将得到OR值和p值:

实际上我对Fisher的理解也很不到位,只能把我知道的告诉你,希望对你有所帮助!

python概率分析_请问如何通过Python做R*C列表的Fisher确切概率分析?相关推荐

  1. 有没有必要买python课_请问自学 Python 有必要买课程吗?

    既然说到了找工作,那肯定是要去培训班的. 网上优质.实用的免费课程有很多,但都是基础教程,真正能找到Python正式程序员水平的课程,根本没有! 那些高赞回答中的教程中高阶部分,也是收费教程!或者看书 ...

  2. 3 x 10的python表达式_这道数学题用PYTHON编程语言怎么写? 编程语言python是用

    我觉着,这个应该这样解决比较符合计算机解题思路. 下面的回答的,思考的东西太多. # -*- coding: utf-8 -*- __author__ = 'lpe234' __date__ = '2 ...

  3. 零基础学python 视频_全网最全Python视频教程真正零基础学习Python视频教程 490集...

    Python Web开发-进阶提升 490集超强Python视频教程 真正零基础学习Python视频教程 [课程简介] 这是一门Python Web开发进阶课程,手把手教你用Python开发完整的商业 ...

  4. 为什么要学python语言_我们为什么要学习Python语言?

    原标题:我们为什么要学习Python语言? 聊到我们为什么要学习Python语言?小编不禁又想起大佬潘石屹准备开启Python学习旅程时所发布的微博. 我们为什么要学习Python语言? 在农业社会时 ...

  5. anaconda3对应的python版本_利用SWIG实现Python调用C/C++!速度快就算了,代码量还特少!...

    SWIG简介 SWIG是Simplified Wrapper and Interface Generator的简称,它是一个能将C和C++的程序与其他各种高级语言诸如Perl,Python,Ruby和 ...

  6. python图像分类_用于实现用python和django编写的图像分类的Keras UI

    KerasUI是一种可视化工具,可以在图像分类中轻松训练模型,并允许将模型作为服务使用,只需调用API. https://github.com/zeppaman/KerasUI 主要特点: 用oaut ...

  7. 计算机导论python知识点_如何系统地自学 Python?

    我是自学的Python.从对Python一无所知,到在博客上写Python相关的系列文章(Python快速教程),前后有将近三年的时间.期间有不少门槛,但也充满乐趣.乐趣是自学的最大动力.Python ...

  8. pythoncookbook和流畅的python对比_为什么你学Python效率比别人慢?因为你没有这套完整的学习资料...

    以下资源免费获取方式! 关注!转发!私信"资料"即可免费领取! 入门书籍 1.<Python基础教程>(Beginning Python From Novice to ...

  9. 杭州哪里学python好_杭州哪里学python好

    1. Python 易于学习 好吧,相较于其它许多你可以拿来用的编程语言而言,它"更容易一些". Python 的语言没有多少仪式化的东西,所以就算不是一个 Python 专家,你 ...

  10. 使用python数据分析_我如何使用Python分析《权力游戏》

    使用python数据分析 By Rocky Kev 洛基·凯夫(Rocky Kev) I wanted to learn Python for a long time, but I could nev ...

最新文章

  1. .NET Micro Framework动态调用C/C++底层代码(原理篇)
  2. linux 生成字母序列,讲解在Linux平台中的Bash序列表达式(Sequence Expression)
  3. Golang map输出排序
  4. java sqlmap_sqlmap 学习指南
  5. PHP怎样表示几时几分,PHP将时间戳转换为刚刚、N分钟前、今天几点几分、昨天几点几分......
  6. java ajax传输图片_Java使用Ajax实现跨域上传图片功能
  7. 【李宏毅2020 ML/DL】补充:Structured Learning: Introduction Structured Linear Model
  8. Javascript的一个生产PDF的库: unicode和中文问题的解决
  9. Java基础 IO流——第三部分
  10. disc性格测试cs适合职业_DISC职业性格测试
  11. matlab 结果导入excel,Matlab导入Excel文件中的数据的详细教程分享
  12. ableton 中文_Ableton live 中文 PDF.pdf
  13. 设置 Google Analytics(分析)全局网站统计代码
  14. 反射Modifier
  15. Notepad++ 提示:An exception occurred due to plugin: HexEditor.dll 出错
  16. 使用ffmpeg读取本地文件,进行推流
  17. 创业失败后的他们,是怎么找工作的?
  18. 最全Java各修饰符详解(接口、类、方法、变量)
  19. css 动画 呼吸 呼吸灯 效果
  20. AUTOSAR 通信栈配置

热门文章

  1. Mapgis67软件实用操作
  2. aspectjweaver和aspectjrt的作用?(作用、说明、案例)
  3. 教你如何修改树莓派的时区和网络对时
  4. 面对不稳定的百度网盘离线下载,这些备用方案可以帮你应急
  5. 我为什么要学习 Python(2020版)
  6. 饥荒一直服务器没有响应,饥荒总是启动服务器进不去 | 手游网游页游攻略大全...
  7. php 中文 验证码,PHP中文验证码函数
  8. edptrayicon怎么卸载_谁知道怎样删除北信源内网管理软件客户端?
  9. 内网穿透软件NPS--客户端NPC SDK交叉编译ARM64位库
  10. C#最简单的NetMeeting