知識グラフ全体このページを中心とした周辺知識グラフ
応用分野: 変動の分解

線形単回帰における「変動の分解」

データNo. データ X
(説明変数)
データ Y
(目的変数)
1 x 1 y 1
2 x 2 y 2
3 x 3 y 3
n x n y n

X Y のデータの組が n 個あるとする. このデータから求められる線形単回帰式

y ^ = a x + b ・・・・・・(1)

ただし, b0

で表されるとする.

備考 x i 説明変数の実際の値, y i 目的変数の実際の値, y ^ i を予測値(回帰式より求めた値), ε i 残差とし.これらの平均値 y ¯ y ^ ¯ ε ¯ で表わす.

目的変数 Y分散 σ y 2 = 1 n i = 1 n y i y ¯ 2 ・・・・・・(2)

予測値 Y ^ 分散 σ y ^ 2 = 1 n i=1 n y ^ i y ^ ¯ 2 ・・・・・・(3)

残差 ε 分散 σ ε 2 = 1 n i=1 n ε i ε ¯ 2 ・・・・・・(4)

ただし

y ¯ = 1 n i = 1 n y i  ・・・・・・(5)

y ^ ¯ = 1 n i=1 n y ^ i  ・・・・・・(6)

ε ¯ = 1 n i=1 n ε i  ・・・・・・(7)

y ^ i = a x i + b  ・・・・・・(8)

ε i = y i y ^ i  ・・・・・・(9)

とおくと

σ y 2 = σ y ^ 2 + σ ε 2  ・・・・・・(10)

が成り立つ.(10)の両辺をn倍すると,(2) ,(3),(4)より

i=1 n y i y ¯ 2 = i=1 n y ^ i y ^ ¯ 2 + i=1 n ε i ε ¯ 2  ・・・・・・(11)

備考: y ¯ = y ^ ¯ ε ¯ =0 (証明の(19),(25)を参照のこと)

の関係が得られる.

i=1 n y i y ¯ 2 全変動(Total Variation),あるいは,全平方和(Total Sum of Squares)(以下,TSSと表記)

i=1 n y ^ i y ^ ¯ 2 回帰変動(Regression Variation),あるいは,回帰平方和(Regression Sum of Squares,あるいは,Sum of Squares due to Regression)(以下,SSRと表記)

i = 1 n ε i ε ¯ 2 残差変動(Residual Variation),あるいは,残差平方和(Residual Sum of Squares ( RSS),あるいは,Sum of Squared Errors)(以下,SSEと表記)

(11)を書き替えると

全変動(TSS)=回帰変動(SSR)+残差変動(SSE) ・・・・・・(12)

のような表現もある.(11),(12)の関係を変動の分解という.

■証明

(1)が線形単回帰式より

a= σ xy σ x 2  ・・・・・・(13)

b= y ¯ a x ¯  ・・・・・・(14)

ただし

σ x 2 = 1 n i=1 n x i x ¯ 2  ・・・・・・(15)

σ xy = 1 n i=1 n x i x ¯ y i y ¯  ・・・・・・(16)

x ¯ = 1 n i=1 n x i  ・・・・・・(17)

である.

とおくと,(6),(8),(17)より

y ^ ¯ = 1 n i=1 n y ^ i = 1 n i=1 n a x i +b =a 1 n i=1 n x i +b =a x ¯ +b  ・・・・・・(18)

となる.

(18)に(14)を代入すると

y ^ ¯ =a x ¯ + y ¯ a x ¯ = y ¯  ・・・・・・(19)

(7),(9),(19)より

ε ¯ = 1 n i=1 n ε i = 1 n i=1 n y i y ^ i = 1 n i=1 n y i 1 n i=1 n y ^ i = y ¯ y ^ ¯ =0  ・・・・・・(20)

となる.

(2)と(9)より

σ y 2 = 1 n i=1 n y i y ¯ 2

= 1 n i=1 n y ^ i + ε i y ¯ 2

= 1 n i=1 n y ^ i y ¯ + ε i 2

= 1 n i=1 n y ^ i y ¯ 2 +2 y ^ i y ¯ ε i + ε i 2

= 1 n i=1 n y ^ i y ¯ 2 + 2 n i=1 n y ^ i y ¯ ε i + 1 n i=1 n ε i 2  ・・・・・・(21)

第1項 1 n i=1 n y ^ i y ¯ 2 について

(19),(3)より

1 n i=1 n y ^ i y ¯ 2 = 1 n i=1 n y ^ i y ^ ¯ 2 = σ y ^ 2  ・・・・・・(22)

第2項 2 n i=1 n y ^ i y ¯ ε i について

(8),(14)より

y ^ i y ¯ = a x i +b a x ¯ +b =a x i x ¯  ・・・・・・(23)

(9),(8),(14)より

ε i = y i a x i +b = y i a x i + y ¯ a x ¯ = y i y ¯ a x i x ¯  ・・・・・・(24)

(23),(24)を第2項に代入すると

2 n i=1 n y ^ i y ¯ ε i = 2 n i=1 n a x i x ¯ y i y ¯ a x i x ¯

=a 2 n i=1 n x i x ¯ y i y ¯ a 2 2 n i=1 n x i x ¯ 2

(16),(15)より

=2a σ xy a σ x 2

(13)を代入すると

=2 σ xy σ x 2 σ xy σ xy σ x 2 σ x 2

=0  ・・・・・・(25)

となる.

第3項 1 n i=1 n ε i 2 について

(20),(4)より

1 n i=1 n ε i 2 = 1 n i=1 n ε i ε ¯ 2 = σ ε 2  ・・・・・・(26)

(22),(25),(26)より

σ y ^ 2 + σ ε 2

となり,(10)が成り立つことが導かれた.

(10)の両辺をn倍すると,(11)が得られる.

 

ホーム>>カテゴリー分類>>確率>>統計>>線形単回帰における「変動の分解」

 最終更新日: 2026年5月17日

[ページトップ]

金沢工業大学

利用規約,誤植や誤りなどがございましたら、こちらからご連絡ください。

google translate (English version)