
作者 | Amond Lee
即使是沒有任何統(tǒng)計學基礎的讀者朋友可能也聽說過「p 值」,但是鮮有文章能夠清楚解釋 p 值是什么,以及 p 值在統(tǒng)計學中的作用。本文是 TowardDataScience 的一篇博文,作者條理清楚地解釋了 p 值的相關內容,并給出了一個簡單的例子,適合讀者參考。
還記得我作為暑期實習生第一次在 CERN 海外實習時,大多數(shù)人都在討論,要超過「5-sigma」閾值(這意味著 p 值為 0.0000003)才能確認發(fā)現(xiàn)了希格斯玻色子。
那時我對 p 值、假設檢驗甚至統(tǒng)計顯著一無所知。
直到進入數(shù)據(jù)科學領域后,我終于意識到了 p 值的含義,以及在某些實驗中,p 值是如何成為決策工具的一部分的。
因此,我決定在這篇文章中解釋什么是 p 值以及如何在假設檢驗中使用 p 值。希望能幫你更好、更直觀地理解 p 值。
本文共分四個部分,從假設檢驗到理解 p 值,以及根據(jù) p 值指導我們的決策過程。我強烈建議你仔細閱讀全文,以便詳細地了解 p 值:
在討論 p 值的意義之前,我們先理解一下假設檢驗。在假設檢驗中,常用 p 值確定結果的統(tǒng)計顯著性。
我們的最終目標是確定結果的統(tǒng)計顯著性。而統(tǒng)計顯著性建立在這 3 個簡單概念之上:
假設檢驗是用來通過一組數(shù)據(jù)檢驗針對總體的聲明(零假設)有效性的。如果零假設不成立,我們就會相信備擇假設。
換句話說,我們需要提出聲明(零假設),并用樣本數(shù)據(jù)來檢驗聲明是否有效。如果聲明是無效的,就選擇備擇假設。就這么簡單。
而要知道聲明是否有效,就要用 p 值來衡量證據(jù)的強度,從而了解到它是否有統(tǒng)計顯著性。如果證據(jù)支持備擇假設,那就拒絕零假設并接受備擇假設。后面的章節(jié)中會解釋這些內容。
我們舉個例子來更清晰地說明這一概念,這個例子會貫穿全文同時說明其他概念。
假設某個披薩店聲稱,他們的平均配送時間小于等于 30 分鐘,但你認為他們的配送時間不止 30 分鐘。所以你做了假設檢驗,對配送時間隨機采樣來檢驗這一說法:
這里的目標是確定樣本數(shù)據(jù)中的證據(jù)能更好地支持哪種假設(零假設或備擇假設)。
本例中用的是單尾檢驗,因為我們只想知道平均配送時間是否大于 30 分鐘。
因為配送時間小于等于 30 分鐘都是可以接受的,因此我們忽略另一個方向的可能性。這里想要檢驗的是平均配送時間是否會大于 30 分鐘。換句話說,我們想知道披薩店是否在某種角度上騙了我們。
假設檢驗的常用方法之一是 Z 檢驗。這里我們不討論細節(jié),因為我們想要先理解表面的內容,然后再深入。
平均值為 μ 標準差為 σ 的正態(tài)分布
正態(tài)分布是用來觀察數(shù)據(jù)分布的概率密度函數(shù)。
正態(tài)分布有兩個參數(shù)——平均值(μ)和標準差(σ)。
均值是分布的集中趨勢。它決定了正態(tài)分布峰值的位置。標準差是衡量可變性的標準,它決定了均值到值的下降幅度。
正態(tài)分布通常和 68-95-99.7 規(guī)則(上圖所示)相關:
還記得文章開頭說的發(fā)現(xiàn)希格斯玻色子的「5-sigma」閾值嗎?在科學家證實發(fā)現(xiàn)希格斯玻色子之前,5-sigma 約為數(shù)據(jù)的「99.9999426696856%」。設置這么嚴格的閾值是為了避免潛在的錯誤信號。
好了?,F(xiàn)在你可能想知道「正態(tài)分布是如何應用在假設檢驗中的」。
因為是用 Z 檢驗進行假設檢驗的,因此要計算 Z 分數(shù)(用于檢驗統(tǒng)計量),這是數(shù)據(jù)點到平均值的標準偏差數(shù)。在本文的例子中,每個數(shù)據(jù)點都是收集到的披薩配送時間。
計算每個數(shù)據(jù)點的 Z 分數(shù)的公式。
對每個披薩配送時間點計算 Z 分數(shù),并繪制出標準正態(tài)分布曲線時,x 軸上的單位從分鐘變成了標準差單位,因為已經(jīng)通過計算(變量減去平均值再除以標準差,見上述公式)將變量標準化了。
標準正態(tài)分布曲線是很有用的,因為我們可以比較測試結果和在標準差中有標準單位的「正態(tài)」總體,特別是在變量的單位不同的情況下。
Z 分數(shù)的標準正態(tài)分布
Z 分數(shù)可以告訴我們整個數(shù)據(jù)相對于總體平均值的位置。
我喜歡 Will Koehrsen 的說法——Z 分數(shù)越高或越低,結果就越不可能偶然發(fā)生,結果就越有可能有意義。
但多高(低)才足以說明結果是有意義的呢?
這就是解決這個難題的最后一片拼圖——p 值。根據(jù)實驗開始前設定的顯著水平(alpha)檢驗結果是否具有統(tǒng)計學意義。
什么是 P 值
與其用維基百科給出的定義來解釋 p 值,不如用文中的披薩配送時間為例來解釋它。
對披薩配送時間隨機采樣,目的是檢查平均配送時間是否大于 30 分鐘。如果最終的結果支持披薩店的說法(平均配送時間小于等于 30 分鐘),那就接受零假設。否則,就拒絕零假設。
因此,p 值的工作就是回答這個問題:
如果我生活在披薩配送時間小于等于 30 分鐘(零假設成立)的世界中,那我在真實世界中得到的證據(jù)有多令人驚訝?
p 值用數(shù)字(概率)回答了這一問題。
p 值越低,證據(jù)越令人驚訝,零假設越荒謬。
當零假設很荒謬的時候還能做什么?可以拒絕零假設并轉而選擇備擇假設。
如果 p 值低于之前定義的顯著水平(人們一般將它稱為 alpha,但我將它稱之為荒謬閾值——別問為什么,我只是覺得這樣更容易理解),那么就可以拒絕零假設。
現(xiàn)在我們理解了 p 值是什么意思。接下來把 p 值用到文中的例子中。
現(xiàn)在已經(jīng)抽樣得到了一些配送時間,計算后發(fā)現(xiàn)平均配送時間要長 10 分鐘,p 值為 0.03。
這意味著在披薩配送時間小于等于 30 分鐘(零假設成立)的世界中,由于隨機噪聲的影響,我們有 3% 的概率會看到披薩配送時間延長了至少 10 分鐘。
p 值越低,結果越有意義,因為它不太可能是由噪聲引起的。
大多數(shù)人對于 p 值都有一個常見的誤解:
p 值為 0.03 意味著有 3%(概率百分比)的結果是偶然決定的——這是錯誤的。
人們都想得到確切的答案(包括我),而這也是我在很長時間內都對 p 值的解釋感到困惑的原因。
p 值不能證明任何事。這只是一種根據(jù)驚訝程度做出合理決策的基礎方法。Cassie Kozyrkov
我們是如何用 0.03 的 p 值來做出合理決策的(重點):
到現(xiàn)在為止,你可能已經(jīng)注意到了,在上面的例子中,p 值不能證明或決定任何事。
在我看來,當結果有統(tǒng)計學意義時,p 值可以作為挑戰(zhàn)初始信念(零假設)的工具。在我們認為自己的信念荒謬(假設 p 值表明結果具有統(tǒng)計顯著性)的那一刻,就放棄了自己的初始信念(拒絕零假設)并做出了更合理的決定。
統(tǒng)計顯著性
這是最后一步,將所有內容放在一起,并檢驗結果是否有統(tǒng)計學意義。
只有 p 值是不夠的,還要設定閾值(即顯著水平——alpha)。為了避免偏差,實驗開始之前就應該設定 alpha。如果觀測的 p 值小于 alpha,那就可以得出結論——結果具有統(tǒng)計顯著性。
經(jīng)驗法則一般將 alpha 設定為 0.05 或 0.01(同樣,值取決于你的問題)。
如上文所述,假設在實驗開始前將 alpha 設置為 0.05,得到的結果具有統(tǒng)計顯著性,因為 p 值(0.03)小于 alpha。
為便于參考,整個實驗的基本步驟如下:
步驟 5 計算檢驗統(tǒng)計量的公式。
數(shù)據(jù)分析咨詢請掃描二維碼
若不方便掃碼,搜微信號:CDAshujufenxi
SQL Server 中 CONVERT 函數(shù)的日期轉換:從基礎用法到實戰(zhàn)優(yōu)化 在 SQL Server 的數(shù)據(jù)處理中,日期格式轉換是高頻需求 —— 無論 ...
2025-09-18MySQL 大表拆分與關聯(lián)查詢效率:打破 “拆分必慢” 的認知誤區(qū) 在 MySQL 數(shù)據(jù)庫管理中,“大表” 始終是性能優(yōu)化繞不開的話題。 ...
2025-09-18CDA 數(shù)據(jù)分析師:表結構數(shù)據(jù) “獲取 - 加工 - 使用” 全流程的賦能者 表結構數(shù)據(jù)(如數(shù)據(jù)庫表、Excel 表、CSV 文件)是企業(yè)數(shù)字 ...
2025-09-18DSGE 模型中的 Et:理性預期算子的內涵、作用與應用解析 動態(tài)隨機一般均衡(Dynamic Stochastic General Equilibrium, DSGE)模 ...
2025-09-17Python 提取 TIF 中地名的完整指南 一、先明確:TIF 中的地名有哪兩種存在形式? 在開始提取前,需先判斷 TIF 文件的類型 —— ...
2025-09-17CDA 數(shù)據(jù)分析師:解鎖表結構數(shù)據(jù)特征價值的專業(yè)核心 表結構數(shù)據(jù)(以 “行 - 列” 規(guī)范存儲的結構化數(shù)據(jù),如數(shù)據(jù)庫表、Excel 表、 ...
2025-09-17Excel 導入數(shù)據(jù)含缺失值?詳解 dropna 函數(shù)的功能與實戰(zhàn)應用 在用 Python(如 pandas 庫)處理 Excel 數(shù)據(jù)時,“缺失值” 是高頻 ...
2025-09-16深入解析卡方檢驗與 t 檢驗:差異、適用場景與實踐應用 在數(shù)據(jù)分析與統(tǒng)計學領域,假設檢驗是驗證研究假設、判斷數(shù)據(jù)差異是否 “ ...
2025-09-16CDA 數(shù)據(jù)分析師:掌控表格結構數(shù)據(jù)全功能周期的專業(yè)操盤手 表格結構數(shù)據(jù)(以 “行 - 列” 存儲的結構化數(shù)據(jù),如 Excel 表、數(shù)據(jù) ...
2025-09-16MySQL 執(zhí)行計劃中 rows 數(shù)量的準確性解析:原理、影響因素與優(yōu)化 在 MySQL SQL 調優(yōu)中,EXPLAIN執(zhí)行計劃是核心工具,而其中的row ...
2025-09-15解析 Python 中 Response 對象的 text 與 content:區(qū)別、場景與實踐指南 在 Python 進行 HTTP 網(wǎng)絡請求開發(fā)時(如使用requests ...
2025-09-15CDA 數(shù)據(jù)分析師:激活表格結構數(shù)據(jù)價值的核心操盤手 表格結構數(shù)據(jù)(如 Excel 表格、數(shù)據(jù)庫表)是企業(yè)最基礎、最核心的數(shù)據(jù)形態(tài) ...
2025-09-15Python HTTP 請求工具對比:urllib.request 與 requests 的核心差異與選擇指南 在 Python 處理 HTTP 請求(如接口調用、數(shù)據(jù)爬取 ...
2025-09-12解決 pd.read_csv 讀取長浮點數(shù)據(jù)的科學計數(shù)法問題 為幫助 Python 數(shù)據(jù)從業(yè)者解決pd.read_csv讀取長浮點數(shù)據(jù)時的科學計數(shù)法問題 ...
2025-09-12CDA 數(shù)據(jù)分析師:業(yè)務數(shù)據(jù)分析步驟的落地者與價值優(yōu)化者 業(yè)務數(shù)據(jù)分析是企業(yè)解決日常運營問題、提升執(zhí)行效率的核心手段,其價值 ...
2025-09-12用 SQL 驗證業(yè)務邏輯:從規(guī)則拆解到數(shù)據(jù)把關的實戰(zhàn)指南 在業(yè)務系統(tǒng)落地過程中,“業(yè)務邏輯” 是連接 “需求設計” 與 “用戶體驗 ...
2025-09-11塔吉特百貨孕婦營銷案例:數(shù)據(jù)驅動下的精準零售革命與啟示 在零售行業(yè) “流量紅利見頂” 的當下,精準營銷成為企業(yè)突圍的核心方 ...
2025-09-11CDA 數(shù)據(jù)分析師與戰(zhàn)略 / 業(yè)務數(shù)據(jù)分析:概念辨析與協(xié)同價值 在數(shù)據(jù)驅動決策的體系中,“戰(zhàn)略數(shù)據(jù)分析”“業(yè)務數(shù)據(jù)分析” 是企業(yè) ...
2025-09-11Excel 數(shù)據(jù)聚類分析:從操作實踐到業(yè)務價值挖掘 在數(shù)據(jù)分析場景中,聚類分析作為 “無監(jiān)督分組” 的核心工具,能從雜亂數(shù)據(jù)中挖 ...
2025-09-10統(tǒng)計模型的核心目的:從數(shù)據(jù)解讀到?jīng)Q策支撐的價值導向 統(tǒng)計模型作為數(shù)據(jù)分析的核心工具,并非簡單的 “公式堆砌”,而是圍繞特定 ...
2025-09-10