python概率分析_请问如何通过Python做R*C列表的Fisher确切概率分析?
因为你这里主要想用到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确切概率分析?相关推荐
- 有没有必要买python课_请问自学 Python 有必要买课程吗?
既然说到了找工作,那肯定是要去培训班的. 网上优质.实用的免费课程有很多,但都是基础教程,真正能找到Python正式程序员水平的课程,根本没有! 那些高赞回答中的教程中高阶部分,也是收费教程!或者看书 ...
- 3 x 10的python表达式_这道数学题用PYTHON编程语言怎么写? 编程语言python是用
我觉着,这个应该这样解决比较符合计算机解题思路. 下面的回答的,思考的东西太多. # -*- coding: utf-8 -*- __author__ = 'lpe234' __date__ = '2 ...
- 零基础学python 视频_全网最全Python视频教程真正零基础学习Python视频教程 490集...
Python Web开发-进阶提升 490集超强Python视频教程 真正零基础学习Python视频教程 [课程简介] 这是一门Python Web开发进阶课程,手把手教你用Python开发完整的商业 ...
- 为什么要学python语言_我们为什么要学习Python语言?
原标题:我们为什么要学习Python语言? 聊到我们为什么要学习Python语言?小编不禁又想起大佬潘石屹准备开启Python学习旅程时所发布的微博. 我们为什么要学习Python语言? 在农业社会时 ...
- anaconda3对应的python版本_利用SWIG实现Python调用C/C++!速度快就算了,代码量还特少!...
SWIG简介 SWIG是Simplified Wrapper and Interface Generator的简称,它是一个能将C和C++的程序与其他各种高级语言诸如Perl,Python,Ruby和 ...
- python图像分类_用于实现用python和django编写的图像分类的Keras UI
KerasUI是一种可视化工具,可以在图像分类中轻松训练模型,并允许将模型作为服务使用,只需调用API. https://github.com/zeppaman/KerasUI 主要特点: 用oaut ...
- 计算机导论python知识点_如何系统地自学 Python?
我是自学的Python.从对Python一无所知,到在博客上写Python相关的系列文章(Python快速教程),前后有将近三年的时间.期间有不少门槛,但也充满乐趣.乐趣是自学的最大动力.Python ...
- pythoncookbook和流畅的python对比_为什么你学Python效率比别人慢?因为你没有这套完整的学习资料...
以下资源免费获取方式! 关注!转发!私信"资料"即可免费领取! 入门书籍 1.<Python基础教程>(Beginning Python From Novice to ...
- 杭州哪里学python好_杭州哪里学python好
1. Python 易于学习 好吧,相较于其它许多你可以拿来用的编程语言而言,它"更容易一些". Python 的语言没有多少仪式化的东西,所以就算不是一个 Python 专家,你 ...
- 使用python数据分析_我如何使用Python分析《权力游戏》
使用python数据分析 By Rocky Kev 洛基·凯夫(Rocky Kev) I wanted to learn Python for a long time, but I could nev ...
最新文章
- .NET Micro Framework动态调用C/C++底层代码(原理篇)
- linux 生成字母序列,讲解在Linux平台中的Bash序列表达式(Sequence Expression)
- Golang map输出排序
- java sqlmap_sqlmap 学习指南
- PHP怎样表示几时几分,PHP将时间戳转换为刚刚、N分钟前、今天几点几分、昨天几点几分......
- java ajax传输图片_Java使用Ajax实现跨域上传图片功能
- 【李宏毅2020 ML/DL】补充:Structured Learning: Introduction Structured Linear Model
- Javascript的一个生产PDF的库: unicode和中文问题的解决
- Java基础 IO流——第三部分
- disc性格测试cs适合职业_DISC职业性格测试
- matlab 结果导入excel,Matlab导入Excel文件中的数据的详细教程分享
- ableton 中文_Ableton live 中文 PDF.pdf
- 设置 Google Analytics(分析)全局网站统计代码
- 反射Modifier
- Notepad++ 提示:An exception occurred due to plugin: HexEditor.dll 出错
- 使用ffmpeg读取本地文件,进行推流
- 创业失败后的他们,是怎么找工作的?
- 最全Java各修饰符详解(接口、类、方法、变量)
- css 动画 呼吸 呼吸灯 效果
- AUTOSAR 通信栈配置
热门文章
- Mapgis67软件实用操作
- aspectjweaver和aspectjrt的作用?(作用、说明、案例)
- 教你如何修改树莓派的时区和网络对时
- 面对不稳定的百度网盘离线下载,这些备用方案可以帮你应急
- 我为什么要学习 Python(2020版)
- 饥荒一直服务器没有响应,饥荒总是启动服务器进不去 | 手游网游页游攻略大全...
- php 中文 验证码,PHP中文验证码函数
- edptrayicon怎么卸载_谁知道怎样删除北信源内网管理软件客户端?
- 内网穿透软件NPS--客户端NPC SDK交叉编译ARM64位库
- C#最简单的NetMeeting