機率與統計·信賴區間

信賴區間

95%信賴區間是由樣本透過一種方法算出的區間;重複抽樣時,這種方法有95%的機率涵蓋真值。這裡每個區間來自一份含n個值的樣本;藍線是真實平均數μ,黃色區間有的涵蓋它,有的沒有。

這幅圖中的信賴區間

信賴區間用一份樣本算出的範圍,估計未知的數值;這裡估計的是母體平均數μ\mu。從樣本產生區間的規則須事先固定;其信賴水準(例如95%)就是反覆抽樣時,區間涵蓋真值的比例。本頁每份樣本由常態母體中nn次獨立抽取的值組成,母體的μ=50\mu = 50、σ=10\sigma = 10;堆疊圖每一列是一份樣本的區間,藍線則是μ\mu。

為何不能說某個區間有95%的機率包含μ\mu?這裡μ\mu是固定數值,隨樣本改變的是區間。抽樣前,算出的區間涵蓋μ\mu的機率是95%;抽樣後,這個區間不是包含50,就是不包含。帶狀區域說明了σ\sigma已知時此規則為何成立:區間xˉ±ME\bar x \pm \mathrm{ME}涵蓋μ\mu,若且唯若∣xˉ−μ∣≤ME|\bar x - \mu| \le \mathrm{ME},也就是xˉ\bar x落在同寬的帶狀區域μ±ME\mu \pm \mathrm{ME}內。因為Xˉ\bar X服從標準差為σ/n\sigma/\sqrt{n}的常態分布,當ME=z∗σ/n\mathrm{ME} = z^*\sigma/\sqrt{n}時,它有95%的機率落在帶內。若要給μ\mu本身一個機率,就需要μ\mu的先驗分布;那是貝氏可信區間,是另一種建構方式。

寬度由什麼決定?σ\sigma已知時,半寬就是誤差範圍z∗σ/nz^*\sigma/\sqrt{n};z∗z^*使標準常態分布兩端的尾部各留下剩餘機率的一半。越接近100%,再往前一步所需增加的寬度越大:z∗z^*在90%時約為1.645,99%時約為2.576,99.9%時約為3.29,沒有有限區間能百分之百保證涵蓋。平均數的標準誤為σ/n\sigma/\sqrt{n}(見抽樣分布),所以水準固定時,寬度與1/n1/\sqrt{n}成正比縮小。涵蓋率不隨nn改變,因為誤差範圍始終以xˉ\bar x的散布程度為單位。

為何σ\sigma未知就需要tt?以樣本的ss取代σ\sigma,會讓中心和寬度都變成隨機的。ss偏小的樣本會得到較短的區間,因此未涵蓋真值的情況集中在短區間。若仍使用z∗z^*,涵蓋率為2Fn−1(z∗)−12F_{n-1}(z^*) - 1,低於設定水準;其中Fn−1F_{n-1}是自由度為n−1n - 1的學生氏tt分布的累積分布函數。對常態資料而言,這正是(Xˉ−μ)/(S/n)(\bar X - \mu)/(S/\sqrt{n})的分布。它的尾部比常態分布厚,因此分位點t∗t^*大於z∗z^*,而xˉ±t∗s/n\bar x \pm t^* s/\sqrt{n}的涵蓋率恰好等於設定水準。只修正ss的平均值還不夠:s/c4s/c_4雖然是σ\sigma的不偏估計量,但配合z∗z^*、設定95%時,在n=5n = 5的情況下涵蓋率仍只有約89.5%。隨著nn增大,ss趨近σ\sigma,t∗t^*也趨近z∗z^*。

觀察重點

  • 區間估計的是平均數,而非資料的散布範圍。σ\sigma已知、n=25n = 25、水準95%時,一份樣本平均只有約31%的值落在自己的區間內。
  • 當方法的真實涵蓋率為95%時,最近100個區間的涵蓋比例會波動,標準誤約為2.2個百分點;累計10000份樣本時,標準誤約為0.22個百分點。

操作方式

  • 水準旋鈕依距離100%的遠近配置刻度:90%、99%與99.9%等距,讓區間寬度變化最快的部分也能調整。
  • 點選堆疊圖的一列,或在圖取得焦點時使用↑↓,即可在上方顯示該區間的樣本;執行中點選會暫停。改變水準或公式,會用同一批樣本重算區間;改變nn則需要新樣本。

相關概念