線形回帰

線形回帰（せんけいかいき、）とは、説明変数（独立変数ともいう）に対して目的変数（従属変数、あるいは反応変数ともいう）が線形またはそれから近い値で表される状態。
線形回帰は統計学における回帰分析の一種であり、非線形回帰と対比される。
また線形回帰のうち、説明変数が1つの場合を単純線形回帰、2つ以上の場合を重回帰と呼ばれる。
概要.
線形回帰では，データから推定される線形予測関数を用いて関係性がモデル化される。このようなモデルは線形モデルと呼ばれる。
説明変数（または予測変数）に対して目的変数の条件付き期待値は、アフィン写像で与えられる。（通常は条件付き期待値だが、条件付メジアンまたは他の分位数を用いることもある。）
線形回帰が非線形回帰に比べて用いられる頻度が高いのは、未知のパラメータに線形に依存するモデルの方が、パラメータに非線形に依存するモデルよりもフィッティングが容易で、推定値の統計的性質を決定しやすいためである。
線形回帰が取り扱う範囲は、予測変数の値を与えられた応答の条件付き確率分布に限る。
全ての変数の同時確率分布は多変量解析の領域として、ここでは扱わない。
線形回帰の用途.
線形回帰は多くの実用的な用途があり、大まかには以下の二種類の用途に分類される。
○予測、予想、またはエラーの削減を目的とする。
→線形回帰は、応答変数と説明変数の値の観測されたデータセットに予測モデルを適合させるために使用できる。
説明変数の追加値が収集された場合、このモデルから応答変数を予測できる。
○説明変数の変動に起因する応答変数の変動を説明することを目的とする。
→線形回帰分析を適用して、応答と説明変数の関係の強さを定量化できる。
これにより各説明変数が応答と全く線形関係を持たないかどうかを判断したり、説明変数のどのサブセットに応答に関する冗長な情報が含まれているかを特定できる。
線形モデルのフィッティング方法.
線形回帰モデルは多くの場合、最小二乗法を用いてフィッティングされる。
それ以外のフィッティング方法としては、最小絶対値法や、リッジ回帰（L2ノルムペナルティ）やラッソ回帰（L1ノルムペナルティ）のように、最小二乗コスト関数のペナルティ付きバージョンを最小化する方法などがある。
逆に最小二乗法は、線形モデルではないモデルのフィットにも使用できる。
このように、「最小二乗法」と「線形モデル」という言葉は密接に関連しているが、同義ではない。
基本モデル.
線形回帰モデルは、目的変数 と説明変数 および擾乱項 の関係を以下のようにモデル化したものである。
ここで は切片（「定数」項）、 は各々の説明変数の係数であり、 は説明変数の個数である。線形回帰においては、説明変数の係数および切片の組 をパラメタとするモデルを与える。また、擾乱項 は説明変数 とは独立である。
ベクトル・行列記法を用いれば、線形回帰モデルは以下のように表せる。
線形回帰が「線形」であるのは、目的変数 が説明変数 の係数 に対して線形であるためである。たとえば
という回帰は に対して明らかに線形ではないが、係数 に対して線形であるから、線形回帰の問題に分類される。
の場合、説明変数は1つだけであり回帰パラメタは2つである。上式は以下のようになる。
同等な定式化に、線形回帰を条件付き期待値のモデルとして陽に表すものがある。
ここで、所与の に対する の条件付き確率分布は擾乱項の確率分布に一致する。
線形回帰の種類.
最小二乗モデル.
最小二乗法はカール・フリードリッヒ・ガウスが1820年代に発展させた。本方法は、擾乱項 の振る舞いに次のような仮定をする（ガウス＝マルコフ仮定）。
以上の仮定は、最小二乗法がある意味で最適なパラメタの推定量を与えることを保証する。
説明変数の個数が 個のモデルを考えると、線形回帰によって決定すべきパラメタは係数 と切片 の 個である。目的変数と説明変数の測定結果の組 を1つのデータとし、 個のデータを用いた線形回帰は以下のように表すことができる。
上記の連立方程式は、目的変数の観測値を 成分の列ベクトル 、説明変数の観測値および切片 の係数 () を 行列 、回帰パラメタを 成分の列ベクトル、観測ごとの擾乱を 成分の列ベクトル とすれば、行列の記法を用いて以下のように表せる。
 の場合、回帰パラメタの標準誤差は算出できない。 が より小さい場合、パラメタは算出できない。 
回帰パラメタの推定量は、
で与えられ、ガウス＝マルコフの定理より推定量formula_9は最良線形不偏推定量になる。つまり、任意の線形不偏推定量formula_10に対して
が成立する。
回帰の二乗和 "SSR" は下式で与えられる。
ここで formula_11 であり formula_12 は "n" × 1 の1ベクトル（各要素が1）である。
項 formula_13 は formula_14　と等価である。
誤差の二乗和 "ESS" は下式で与えられる。
二乗和の全和 "TSS' は下式で与えられる。
決定係数, "R"² は下式で与えられる。
擾乱項が正規分布に従うモデル.
以下では擾乱項が互いに独立な平均formula_15, 分散formula_16の正規分布に従うと仮定する。
残差は、観測値とモデルによる予測値の差を表し、以下のように決定される。
この時、統計量formula_17は分散formula_16の不偏推定量（formula_19）になる。また、最小二乗推定量formula_9と統計量formula_21について以下が成立することが知られている。証明は久保川(2017)や解説記事が詳しい。
上記の事実をもとに回帰係数の有意性検定、信頼区間や予測区間を構成できる。
回帰係数の有意性検定.
回帰係数の推定量formula_29は正規分布formula_30に従うことから
は自由度formula_25のformula_33分布に従う。ここでformula_34は行列formula_35の第formula_36成分である。（添え字formula_37は0から始まることに注意。）
これより適当な有意水準formula_38で
を検定することできる。
信頼区間と予測区間.
値 formula_41におけるformula_42 の信頼区間は下式で表される。
同様に値 formula_41におけるformula_42 の予測区間は下式で表される。