R语言估计获胜概率:模拟分析学生多项选择考试通过概率可视化

“获胜概率”的实时计算(或估计)很困难。

由Kaizong Ye,Liao Bao撰写

我们经常在足球比赛中,在选举中看到这种情况。

考虑经典的多项选择考试。在每个问题之后,想象您尝试计算学生通过考试的概率。

×
两个基本定律
加法定律和相乘定律



  • 加法定律



几个互斥事件中任何一个发生的概率是各个事件发生的概率之和。



  • 相乘定律



一系列独立事件都发生的概率是每个独立事件概率的乘积。



有了这两个工具,就足以回答日常生活中大部分关于“几率是多少”的问题。



关于未来的问题
中心极限定理非常简单,但非常强大。



当一个统计学家说出”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

但这 只是 对每一步的最终分数的预测。这不是通过概率的计算!


课程

R语言数据分析挖掘必知必会

从数据获取和清理开始,有目的的进行探索性分析与可视化。让数据从生涩的资料,摇身成为有温度的故事。

立即参加

如果在 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


R语言对布丰投针(蒲丰投针)实验进行模拟和动态可视化生成GIF动画

阅读文章


这比我们之前获得的曲线更不稳定!所以,计算“获胜概率”是一项复杂的工作!


随时关注您喜欢的主题


当然,如果我的学生不抛硬币,情况就略有不同……这是我们得到的结果,如果一半的学生是好的(有2/3的概率答对问题),一半是不好的(1/3的概率)。

如果我们看通过的概率,我们通常不必等到最后(50道题)就知道谁通过了,谁没通过


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498

关注有关新文章的微信公众号


永远不要错过任何见解。当新文章发表时,我们会通过微信公众号向您推送。

技术干货

最新洞察

This will close in 0 seconds