你好，我是林驥。

最近，我正在讀《統(tǒng)計學(xué)關(guān)我什么事》這本書，副標題是：生活中的極簡統(tǒng)計學(xué)，作者小島寬之，是日本的經(jīng)濟學(xué)博士，知名的數(shù)學(xué)隨筆作家。

這是一本關(guān)于「貝葉斯統(tǒng)計學(xué)」的超級入門書，其中沒有復(fù)雜的公式，作者運用簡單的「面積圖」和四則運算，就把復(fù)雜的統(tǒng)計學(xué)理論，講得簡單、生動、有趣，我認為值得學(xué)習(xí)。

1. 一個例子

書中舉了一個很有趣的例子，假如有位女同事給你送了一盒巧克力，你推測她把你當(dāng)成「真命天子」的概率有多大？

根據(jù)貝葉斯推理的方法，只需要簡單的 3 步，就能算出這個概率。

第 1 步

先假設(shè)女同事把你當(dāng)成「真命天子」的概率為10%，這個概率可以根據(jù)自己的主觀感受進行修改；

第 2 步

通過調(diào)查統(tǒng)計知道，職場女性對「真命天子」送出巧克力的平均概率為42.5%，對「無關(guān)路人」送出巧克力的平均概率為22%；

第 3 步

根據(jù)女同事送巧克力的行為，計算她把你當(dāng)成「真命天子」的概率 = 10%*42.5% / (10%*42.5% + 90%*22%) =?17.67%。

也就是說，當(dāng)你認為女同事對你有意思的可能性比較低時，即使她送了巧克力給你，也只是增加了她對你有意思的概率，但并不能因此說明她對你真的有意思。

2. 矩陣樹圖

為了直觀地展現(xiàn)上述推理過程的概率分布情況，我用 pyecharts 制作了一張矩陣樹圖：

# 導(dǎo)入庫from pyecharts import options as optsfrom pyecharts.charts import TreeMapfrom pyecharts.globals import ThemeType 
# 真命天子的概率a = 0.1# 無關(guān)路人的概率b = 1 - a# 送給真命天子的概率c = 0.425# 不送真命天子的概率d = 1 - c# 送給無關(guān)路人的概率e = 0.22# 不送無關(guān)路人的概率f = 1 - e
# 定義數(shù)據(jù)data = [    {        "value": 100,         "name": "后驗概率",        "children": [            {                "value": 100,                "name": "送巧克力",                "children": [                    {"value": round(a*c/(a*c+b*e)*100,2), "name": "送給真命天子"},                    {"value": round(b*e/(a*c+b*e)*100,2), "name": "送給無關(guān)路人"},                ],            }        ],    },    {        "value": 100,        "name": "先驗概率",        "children": [            {                "value": 100,                "name": "條件概率",                "children": [                    {"value": round(a*c*100,2), "name": "送給真命天子"},                    {"value": round(a*d*100,2), "name": "不送真命天子"},                    {"value": round(b*e*100,2), "name": "送給無關(guān)路人"},                    {"value": round(b*f*100,2), "name": "不送無關(guān)路人"},                ],            }        ],    },]
# 繪制矩形樹圖c = (    TreeMap(init_opts=opts.InitOpts(width='800px', height='800px', theme=ThemeType.LIGHT))    .add("推測送巧克力的女同事的心意",          data,         levels=[            opts.TreeMapLevelsOpts(                treemap_itemstyle_opts=opts.TreeMapItemStyleOpts(                    border_color="#00589F", border_width=4, gap_width=4                )            ),            opts.TreeMapLevelsOpts(                color_saturation=[0.3, 0.5],                treemap_itemstyle_opts=opts.TreeMapItemStyleOpts(                    border_color_saturation=0.7, gap_width=2, border_width=2                ),            ),            opts.TreeMapLevelsOpts(                color_saturation=[0.3, 0.5],                treemap_itemstyle_opts=opts.TreeMapItemStyleOpts(                    border_color_saturation=0.6, gap_width=1                ),            ),            opts.TreeMapLevelsOpts(color_saturation=[0.3, 0.5]),        ],    )    .set_global_opts(title_opts=opts.TitleOpts(        title="推測送巧克力的女同事的心意，她把你當(dāng)成真命天子的概率是"+'%.2f%%' % (a*c/(a*c+b*e)*100),        pos_left="center", pos_top="30",    ))    .render("推測送巧克力的女同事的心意.html"))