線性回歸-置信區間為何在中部收窄?

前言

不知道你是否見過這樣的圖:

企鵝分組.png

這是對三組數據的線性回歸,來自經典數據集“帕爾默群島的企鵝”

  • 橫軸(bill_depth)表示嘴峰深度。
  • 縱軸(bill_length)表示嘴峰長度。
  • 不同顏色表示不同種類的企鵝。
    在看到這樣的圖時,你是否會有疑問:置信區間為什么會出現兩邊寬,中間窄的情況咧?
    我們從線性模型開始

線性模型

假設響應變量Y和自變量X_1, X_2, X_3...X_k存在簡單線性關系(兩者可以用最簡單的線性模型描述):
Y = \beta_0 + \beta_1X + \varepsilon
其中\beta_0 , \beta_1回歸參數(回歸系數)\varepsilon誤差項
根據這個式子我們可以得到預測值與自變量之間的關系:
\hat{y}_i = \hat{\beta_1}x_i +\hat{\beta_0}

最小二乘估參數

  • 定義殘差e_i = y_i - \hat{y}_i

如果我們找到一條直線,使得每個預測值和實際值之間的差的平方和(或者絕對值、正負項和)最小,也就是殘差平方和最小, 此時的直線最接近實際的數據,由此而來只需要計算殘差平方和(residual sum of squares ,RSS),并求其取最小值時的\beta_0 , \beta_1, 即可找到擬合直線,前人已經計算過了,使RSS最小的參數估計值是:
\hat{\beta}_1 = \frac{\sum_{i=1}^{n} (x_i - \overline{x})(y_i - \overline{y})}{\sum_{i=1}^{n} (x_i - \overline{x})^2}
\hat{\beta}_0 = \overline{y} - \hat{\beta}_1 \overline{x}

置信區間

某一個預測值\hat{y}的置信區間CI_\hat{y}(confidence interval)可以用以下的式子表達:
CI_{\hat{y}} = \hat{y} \pm t_{\alpha/2, n-2} \cdot SE_{\hat{y}}
其中某一個預測值\hat{y}的標準誤差SE_\hat{y}(standard error)可以用以下的式子表達:
SE_{\hat{y}} = S \sqrt{1 + \frac{1}{n} + \frac{(x - \overline{x})^2}{\sum (x_i - \overline{x})^2}}
其中:

  • S 是殘差的標準誤。
  • n 是樣本大小。
  • \overline{x} 是自變量的均值。
  • x 是特定的自變量。

合在一起就得到了預測值置信區間的表達式:
CI_{\hat{y}} = (\beta_0 + \beta_1x) \pm t_{\alpha/2, n-2} \cdot S \sqrt{1 + \frac{1}{n} + \frac{(x - \overline{x})^2}{\sum (x_i - \overline{x})^2}}

結論

由上式可知當x取到\overline{x}附近時,(x - \overline{x})^2逐漸變小,使得偏差部分整體變小,那么此時的置信區間就會變小,反映在圖上就是中間窄、兩邊寬了;

還有一種直觀的理解是:左側點對右側預測值的影響沒有那么大,右側點對左側預測值的影響也沒有那么大,但是兩邊的點都會為預測中間的值作出貢獻,因此中間的預測值實際上是在獲得了更多信息(兩側的點)后預測得到的,加上回歸直線一定會過(\overline{x}, \overline{y})點,所以預測中間值的信心就很足,置信區間就窄一點。

?著作權歸作者所有,轉載或內容合作請聯系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 227,250評論 6 530
  • 序言:濱河連續發生了三起死亡事件,死亡現場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機,發現死者居然都...
    沈念sama閱讀 97,923評論 3 413
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事。” “怎么了?”我有些...
    開封第一講書人閱讀 175,041評論 0 373
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經常有香客問我,道長,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 62,475評論 1 308
  • 正文 為了忘掉前任,我火速辦了婚禮,結果婚禮上,老公的妹妹穿的比我還像新娘。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 71,253評論 6 405
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上,一...
    開封第一講書人閱讀 54,801評論 1 321
  • 那天,我揣著相機與錄音,去河邊找鬼。 笑死,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 42,882評論 3 440
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 42,023評論 0 285
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后,有當地人在樹林里發現了一具尸體,經...
    沈念sama閱讀 48,530評論 1 331
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 40,494評論 3 354
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發現自己被綠了。 大學時的朋友給我發了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 42,639評論 1 366
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 38,177評論 5 355
  • 正文 年R本政府宣布,位于F島的核電站,受9級特大地震影響,放射性物質發生泄漏。R本人自食惡果不足惜,卻給世界環境...
    茶點故事閱讀 43,890評論 3 345
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 34,289評論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春,著一層夾襖步出監牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 35,552評論 1 281
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人。 一個月前我還...
    沈念sama閱讀 51,242評論 3 389
  • 正文 我出身青樓,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 47,626評論 2 370

推薦閱讀更多精彩內容