“获胜概率”的实时计算(或估计)很困难。
我们经常在足球比赛中,在选举中看到这种情况。
考虑经典的多项选择考试。在每个问题之后,想象您尝试计算学生通过考试的概率。
-
加法定律
几个互斥事件中任何一个发生的概率是各个事件发生的概率之和。
-
相乘定律
一系列独立事件都发生的概率是每个独立事件概率的乘积。
有了这两个工具,就足以回答日常生活中大部分关于“几率是多少”的问题。
当一个统计学家说出”1/10的可能性”这类话时,他就是对未来进行了一次预测。
这或许是对一系列永远都无法检验的事件所做的假设性陈述,或许是对即将发生的事件不掺半点水分的如实解说。
不管是哪种,他都是在对可能的结果进行统计学解释,所有的统计学家所说的话都无非如此。
如果你能够理解以下表述,那么你就具有了像统计黑客一样行动和思考的必备能力:
“如果有10件事情可能要发生,并且这10件事情发生的可能性相等,那么这10件事中任何一件发生的几率是1/10。”
科学研究中充满了可用统计来回答的问题,当然,还有概率定律的运用,但在实验室之外还有很多难题,比愚笨陈旧的科学问题更加重要的问题,比如骰子游戏。
假设你是一名业余赌徒,家里的小孩想要双新鞋子。你下次掷出一对骰子的值会决定你的未来。那你也许想知道骰子扔出各种结果的可能性,而且是非常准确地知道这种可能性!
只凭这两件概率工具,就能回答你可能问到的三类最重要的概率问题。
你提出的问题很可能是以下三种类型之一。
-
下一步出现某个特定结果的可能性是多少?比如,下面会掷出一个7吗?
-
下一步出现某组结果的可能性是多少?比如,下面会出现7或11吗?
-
下一步出现一系列结果的可能性是多少?比如,一对没被动过手脚的骰子真的能够整晚都不出现7吗(我说的是永远都不出现)?我的意思是,那真的可能吗?可能吗?!
在这里考虑我们有 50 个问题的情况。学生在答对 25 个以上时通过。为了模拟,我假设学生在每个问题上只掷硬币,我有 n 个学生,50 个问题
M=matrix
令 Xi,j 表示学生 i在问题 j 的分数。让 Si,j 表示累积分数,即 . 在第 j 步,我可以使用 T^i,j =50×Si,j /j 对最终得分进行某种预测。这是代码
S=apply B
我们可以绘制它
plot(B) abline for(i in 2:n) lines lines
但这 只是 对每一步的最终分数的预测。这不是通过概率的计算!
如果在 j 个问题之后,学生有 25 个正确的答案,那么概率应该是 1——即如果 Si,j ≥25。另一个简单的例子是:如果在j题之后,他直到最后都答对了,他能得到的分数不够,他就会失败。这意味着如果 Si,j +(50−i+1)<25,概率应该是 0。否则,要计算成功的概率,就很简单了。它是当成功的概率实际上是 Si,j /j 时,在 50-j 个问题中获得至少 25-Si,j 正确答案的概率。我们认识到二项式分布的生存概率。然后代码很简单
for(i in 1:50){ for(j in 1:n){ if() P\[i,j\]=1 if() P\[i,j\]=0 if()B\[i,j\]=1-pbinom
所以如果我们绘制它,我们得到
plot(P abline for(i in 2:n) lines lines
这比我们之前获得的曲线更不稳定!所以,计算“获胜概率”是一项复杂的工作!
随时关注您喜欢的主题
当然,如果我的学生不抛硬币,情况就略有不同……这是我们得到的结果,如果一半的学生是好的(有2/3的概率答对问题),一半是不好的(1/3的概率)。
如果我们看通过的概率,我们通常不必等到最后(50道题)就知道谁通过了,谁没通过
可下载资源
关于作者
Kaizong Ye是拓端研究室(TRL)的研究员。在此对他对本文所作的贡献表示诚挚感谢,他在上海财经大学完成了统计学专业的硕士学位,专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。
本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。
非常感谢您阅读本文,如需帮助请联系我们!