機率與統計·p值

p值

p值是在假定虛無假設H0成立時,得到至少與觀測結果同樣極端之結果的機率。這裡機器應平均每袋裝500克;曲線顯示假如它確實如此,n袋的平均重量會落在哪裡。

這幅圖中的p值

假設檢定問的是:假如預設的說法,也就是虛無假設H0H_0成立,眼前的資料是否令人意外。這裡H0H_0表示機器平均每袋裝μ=500\mu = 500克;各袋重量彼此獨立,服從σ=10\sigma = 10克的常態分布。每次研究稱nn袋,算出平均重量xˉ\bar x。p值是在假定H0H_0成立時,得到至少與觀測平均數同樣極端之平均數的機率;「極端」的方向由對立假設H1H_1指定。顯著水準α=0.05\alpha = 0.05的檢定會在p<0.05p < 0.05時發出警報。

p究竟是哪件事的機率?H0H_0成立時,平均數Xˉ\bar X服從平均數500、標準誤10/n10/\sqrt{n}的常態分布,因此觀測平均數的位置是z=(xˉ−500)/(10/n)z = (\bar x - 500)/(10/\sqrt{n})。對H1 ⁣:μ<500H_1\!: \mu < 500,p值是zz左側的面積Φ(z)\Phi(z);對μ>500\mu > 500,是右側的面積;對μ≠500\mu \neq 500,則是兩端尾部的面積2Φ(−∣z∣)2\Phi(-|z|)。重做研究用計數呈現同一件事:機器確實平均每袋裝500克時,平均重量至少這麼極端的研究比例會穩定在p附近。方向必須在看資料前決定。看完後才挑兩個單側p值中較小的那個,H0H_0為真時會有10%的機率發出警報,是宣稱比率的兩倍。方向事先固定時,H0H_0成立下的p在0到1之間均勻分布,因此檢定誤報的機率恰好是α\alpha。

為何p不是H0H_0為真的機率?它完全是在H0H_0成立的世界中計算的:它是P(平均重量至少同樣極端∣H0)P(\text{平均重量至少同樣極端} \mid H_0),不是P(H0∣資料)P(H_0 \mid \text{資料})。反過來推論需要兩項p從未使用的資訊:少裝機器原本有多常見,以及檢定力,也就是檢定對少裝機器發出警報的機率。稱25袋、少裝3克時,檢定力約為Φ(1.5−1.645)\Phi(1.5 - 1.645)。依貝氏定理的勝算形式,一次警報會把「少裝對正常」的勝算乘以檢定力÷α\alpha,這裡約為8.85,使先驗勝算1:4變成約2.2:1。同樣收到警報,少裝機器占80%時,警報中來自正常機器的比例約為2.7%;少裝機器占5%時,則約為68%。這些比例只用了「發出警報」這項資訊;xˉ\bar x的具體數值還能提供更多資訊。

為何很小的差異也可能顯著?觀測到的少裝量固定時,∣z∣=∣xˉ−500∣n/σ|z| = |\bar x - 500|\sqrt{n}/\sigma隨n\sqrt{n}增大,因此稱的袋子越多,左尾p值越小;袋數變為四倍,∣z∣|z|就變為兩倍。在左尾檢定中,比500克少半克,從1083袋起就在α=0.05\alpha = 0.05下顯著。p小只表示假如μ\mu恰好是500,這樣的資料並不常見;不表示差異在實務上重要。反過來也不成立:真的少裝也常常不發出警報,因此p>0.05p > 0.05不能證明機器正常。應搭配信賴區間報告差異的大小:α=0.05\alpha = 0.05的雙側檢定拒絕μ=500\mu = 500,恰好等同於95%區間xˉ±z∗σ/n\bar x \pm z^*\sigma/\sqrt{n}(z∗≈1.96z^* \approx 1.96)不包含500。

觀察重點

  • 上層袋子是示例資料,經安排使平均重量恰好為xˉ\bar x;檢定只使用xˉ\bar x、nn與σ\sigma。
  • 下層橫軸始終涵蓋500兩側各6個標準誤,所以曲線形狀不變。上層橫軸下方的框標示這個範圍,會隨nn增大而縮窄。

操作方式

  • 在下層拖曳可調整xˉ\bar x。移動xˉ\bar x或方向撥桿,會對同一批重做結果重新計數;改變nn則會讓兩種模式重新開始。

相關概念