如何理解贝叶斯
“还不睡?”,室友小憨憨在昨晚十二点多的时候问我。
“你先睡吧,我写个bug”,小憨憨一顿操作,关了灯放下了手机,我知道我明天起来又见不到他了。
和小憨憨成为室友以来四百多天里,我起得比他早的次数一只手数得过来。好吧,我承认没有写bug,写了个理论上快点的程序,一跑竟然比对比程序慢了60多倍。。。这怎么能忍?于是在小憨憨关了灯后我也轻轻地将台灯开到了最暗,调起了我的bug;时间,一点点流逝,伴随着小憨憨舒缓的呼吸声,困意也一点点袭来,思维的效率也降到了极点。此时,我觉得幸福就是关了台灯趴在床上。于是我默默地打开了百度:“幸福是指一个人自我的价值得到满足而产生的喜悦”,显然我对幸福有误解。伴随着隔壁室友的呼噜和梦话,困意也消失了一丝。虽然很困,但我没有拉开窗帘,因为我怕我跑到小憨憨的床上。于是,终于,在两点多,我的程序在不到1s里就给出了正确结果,比对比程序快了两千六百多倍。幸福不幸福我不知道,那一刻喜悦还是有的。啊,终于可以睡觉了。
意外总是很意外,今天醒来的时候我看到了小憨憨,哈哈哈哈。咦?哦,今天原来是周六。
我打破了沉默,“还没走啊?”
“一会儿,昨晚几点睡的?”
“两点多,调出来就睡了”,
“怎么样?”
“快了两千多倍”
“牛逼,是啥呀?”
“一个贝叶斯的程序”
“什么是贝叶斯?”
“今天晚上回来你就知道了”
于是就有了这篇博客。
以下,我们正式来看贝叶斯:
我们就拿笔者工作日到岗的时间来说,正常时间笔者的到岗时间为9点,这也就意味着笔者要在8点多起床,通常来说,起床的时候室友小憨憨已经不在寝室了。当然,笔者也有在7点多早起的时候,这种情况很少发生,要么就是有事,要么就是心血来潮,间歇性踌躇满志,想努力。拿2019年来说,笔者总共早起了10次,其中7次是因为有事,3次是因为间歇性努力症发作。Ok,你也看出来了,笔者早起时有事的概率为
7
/
10
=
70
%
7/10=70\%
7/10=70%,早起时心血来潮的概率为
3
/
10
=
30
%
3/10=30\%
3/10=30%。当然,这是正向概率,什么意思?因为有事,所以我早起了。然而小憨憨有点不一样,怎么说呢?比如在2019年12月31日笔者就早起了,小憨憨微睁着眼问我:“有事?”,看,小憨憨的贝叶斯学得多好,虽然他自己也不知道,甚至还在问我啥事啥是贝叶斯。
我微微一笑,“都2019年最后一天了,还不努力一下?”
我们再来看看贝叶斯公式:
P ( B i ∣ A ) = P ( B i ) P ( A ∣ B i ) ∑ j = 1 n P ( B j ) P ( A ∣ B j ) P({B_i}|A) = \frac{{P({B_i})P(A|{B_i})}}{{\sum\limits_{j = 1}^n {P({B_j})P(A|{B_j})} }} P(Bi∣A)=j=1∑nP(Bj)P(A∣Bj)P(Bi)P(A∣Bi)
好了,让我们来建个模吧:
我们这样表示:
B
1
B_1
B1:有事
B
2
B_2
B2:心血来潮
A
A
A:早起
让我们按照贝叶斯公式来算算笔者2019年的早起情况:
2019年早起的概率:
P
(
A
)
=
10
/
365
P(A)=10/365
P(A)=10/365
2019年有事早起的概率:
P
(
A
∣
B
1
)
=
7
/
365
P(A|B_1)=7/365
P(A∣B1)=7/365
2019年心血来潮早起的概率:
P
(
A
∣
B
2
)
=
3
/
365
P(A|B_2)=3/365
P(A∣B2)=3/365
对,2019年影响早起的只有两种情况,有事和心血来潮,令有事和心血来潮的概率一样大,都为
1
/
2
1/2
1/2,即:
P
(
B
1
)
=
P
(
B
2
)
=
1
/
2
P({B_1})=P(B_2)=1/2
P(B1)=P(B2)=1/2。
纵观整个2019年,假如我某一天早起了,是因为有事的概率为:
P
(
B
1
∣
A
)
=
P
(
A
∣
B
1
)
P
(
B
1
)
P
(
A
)
=
7
/
10
P(B_1|A)= \frac{{P(A|B_1)P(B_1)}}{{P(A)}}=7/10
P(B1∣A)=P(A)P(A∣B1)P(B1)=7/10;(
P
(
B
1
)
P
(
A
∣
B
1
)
+
P
(
B
2
)
P
(
A
∣
B
2
)
=
10
/
365
=
P
(
A
)
P({B_1})P(A|{B_1}) + P({B_2})P(A|{B_2}) = 10/365=P(A)
P(B1)P(A∣B1)+P(B2)P(A∣B2)=10/365=P(A)),是因为心血来潮的概率为:
P
(
B
2
∣
A
)
=
P
(
A
∣
B
2
)
P
(
B
2
)
P
(
A
)
=
3
/
10
P(B_2|A)= \frac{{P(A|B_2)P(B_2)}}{{P(A)}}=3/10
P(B2∣A)=P(A)P(A∣B2)P(B2)=3/10
对,这就是贝叶斯,逆概率,在结果发生的时候发推原因是什么------“执果寻因”。
所以,去年我哪天早起了,70%的概率是因为我有事,30%是因为我间歇性踌躇满志了。
好了,今天的分享就到这里啦~
等等,小憨憨问我“有事?”,他竟然有70%的概率,我的天呐。可以看的出来小憨憨对我有多了解。
我去问问小憨憨!
才疏学浅,难免有错误和不当之处,欢迎交流批评指正!
同时有问题的话欢迎留言或邮箱联系(ljt_IT@163.com)。