重回帰分析のソースを表示

'''重回帰分析'''（じゅうかいきぶんせき）は、[[多変量解析]]の一つ。[[回帰分析]]の[[独立変数]]が複数になったもの。適切な変数を複数選択することで、計算しやすく[[誤差]]の少ない予測式を作ることができる。

一般的によく使われている[[最小二乗法]]、[[線形モデル]]の重回帰は、数学的には[[線形分析]]の一種であり、[[分散分析]]などと数学的に類似している。

==結果の見方==
[[中学生]]を対象に調査を行いその結果を重回帰分析したところ下の式が得られたとする。

: 中学で勉強した時間数 * 3 + 小学生の時代の塾の学習時間数 * 5 + 20 = 知っている英単語の数

この場合、Aさんが中学で100時間、小学生時代20時間勉強していたら 100 * 3 + 20 * 5 + 20 = 英単語420語 を知っているという計算になる。中学で1時間勉強すると平均的には3単語を覚えているという解釈ができる。

ここでは*3や*5という数値を重回帰分析で計算・算出するが、あくまで平均的な値であり個々のサンプルにおいてその通りに一致するとは限らない。例えば、Aさんの英単語数は420語ではなく、実際には450語かもしれない。全体の平均を取ると、3や5という値を取ると全体が最もうまく説明できデータによく適合するということから出てきた数値になる。

また英単語数を決めるのが勉強時間だという関係は、分析者が自分で決めるため絶対的なものではない。あくまで勉強時間が独立変数（説明変数）だと仮定した上で分析している。そのため、予測を行うことはできてもその方向に[[因果関係]]があることは保証されない。

==多重共線性==
独立変数（説明変数）を選択する際、[[マーケティング]]や[[アンケート]]でよく使う一般的な重回帰の場合、複数の説明変数同士は無[[相関係数|相関]]という[[仮定]]が入っている。そのため、説明変数同士が関連性の高い場合、多重共線性と呼ばれる状態になるため、係数が直感に反する値になることがある。

例えば、小学校での定期テスト得点から重回帰で分析する場合、理科の点数を[[従属変数]]に、数学と国語とを説明変数にした場合、数学が増えると理科の点数が増え、国語の点数が高ければ理科の点数が減るといった意味の係数が出ることがある。これは数学と国語との点数の間に強い相関がある（一般に、どちらの成績も学習習慣や[[知能]]の影響を強く受ける）ことで起こりうる。この場合のように説明変数間の相関が高いと係数が不安定になりやすい。

実務的対応としては、一方を除いて分析するのが最も手軽である。また、数学と国語の平均点と、数学と国語の得点の差というように和と差に数字を加工すると、この2つは相関がたいてい低く、かつ解釈しやすい。数学と国語の得点の差は、数学の方が高い生徒の方が理科の点数が高い傾向があるというように理解できるためである。ただし、このような正の相関を持つ変数同士の差得点は元の変数よりも信頼性が落ちるので、サンプル数を増やすなどの対応が求められる。

また、適切な予測力を実質的には持たない変数であっても、説明変数に加えると予測式自体の説明力 (R<sup><small>2</small></sup>) は上がることが多い。そのため、単なるR<sup><small>2</small></sup>ではなく、その分を調整した修正R<sup><small>2</small></sup>を参照する、ステップワイズ法等で投入する説明変数を取捨選択する、[[赤池情報量規準|AIC]]を見るなどの対応が求められる。

==ダミー変数を用いた重回帰分析==
性別や民族といった[[名義尺度]]データを説明変数に用いたい場合は、ダミー変数を導入して重回帰分析を行うことになる。日本国内で見かける[[数量化理論|数量化I類]]は、実質的にはこれと同じ分析である。

ダミー変数は通常、条件ごとに説明変数を作り1と0を持つデータを設定して、分析する。その際、重回帰では一つは回帰式に含めない（多重共線性を避けるためである）。

信号の色による車の通過速度の違いを例に取ると、

: [[信号機]]のある交差点を通過する車の平均時速 (km/h) = 50*青なら1 + 15*黄色なら1 + 0

という式が得られうる。

この場合、青でも黄色でもない条件のデータは、赤になり、計算値は定数の0になる。つまり、赤なら平均時速0km/hとなる。解釈としては交差点を抜ける車はいないという意味になる。また、青の場合なら平均時速は50km/hと推定されたことになる。

男女の体重の予測では、

: 体重 (kg) = 12*男なら1 + 50

という予測式が得られたら、男でない「女」の平均体重は50kgと計算される。男性なら12kg多く、62kgが平均になるという意味になる。ちなみにこれは同じデータを男女別に単純平均したものと一致する。

性別、学年など複数の変数を組み合わせて、分析することもできる。

: 例: 体重 = 性別の違い + 学年

また、「男性で1年生なら1」というように細かく分けてダミー変数を増やして重回帰を行うことも可能ではある。ただし、説明変数の数が大きく増すので、連関の強いダミー変数同士で多重共線性の問題が生じやすいこと、十分な信頼性を確保するためにはサンプル数がかなり求められることなどを考えると、実用性に乏しい。

==ソフトウェア==
* [[SAS Institute|SAS]]、[[Stata]]、[[SPSS]]、[[Microsoft Office]] [[Excel]]といった主要な統計パッケージで、重回帰分析も実行が可能である。

===フリーソフト===
* [http://www.heisei-u.ac.jp/ba/fukui/analysis.html College Analysis]
* [http://www.vector.co.jp/soft/win95/edu/se203904.html 多変量解析入門(95)]
* [[R言語]] - 統計解析言語。重回帰分析だけでなく[[多変量解析]]ほか多くの統計関数を標準装備した[[フリーウェア]]。『モデル式』でモデル記述や当てはめが容易。他アプリケーションのファイル取込やODBC接続対応。[[FDA]]公認。CRANなる仕組で世界の膨大なソフトを無償利用可能。[[可視化]]機能に優れ、日本語対応。[[マルチプラットフォーム]]。[http://www.okada.jp.org/RWiki/index.php?R%A4%CE%B4%F0%CB%DC%A5%D1%A5%C3%A5%B1%A1%BC%A5%B8%C3%E6%A4%CE%B2%F3%B5%A2%A1%A2%CA%AC%BB%B6%CA%AC%C0%CF%B4%D8%BF%F4%B0%EC%CD%F7 Rの基本パッケージ中の回帰、分散分析関数一覧]。重回帰分析はlm関数で行えるほか、独自に書かれた関数もある: [http://aoki2.si.gunma-u.ac.jp/R/mreg2.html][http://aoki2.si.gunma-u.ac.jp/R/All_possible_subset_selection.html]。

==関連する分析手法==
* [[分散分析]] - 重回帰モデルの有意性の検定に用いられる
* [[回帰分析]] - 重回帰分析の説明変数が一つだけの場合
* [[数量化理論|数量化I類]] - ダミー変数を用いた重回帰分析と同等
* [[ロジスティック回帰]] - 目的変数は名義尺度、ロジスティック関数に基づく予測
* [[プロビット分析]] - 目的変数は名義尺度、正規累積関数に基づく予測
* [[正準相関分析]] - 変数群の一方が一変数のみであれば、重回帰分析と同様
* [[パス解析]] - 古典的なパス解析は、複雑なモデルについて重回帰分析の繰り返しから個々の係数を推定
* [[共分散構造分析]]（構造方程式モデリング） - 重回帰モデルに対し、適合度統計量等の算出が可能

{{DEFAULTSORT:しゆうかいきふんせき}}
[[Category:統計学|しゆうかいきふんせき]]
[[Category:数学に関する記事|しゆうかいきふんせき]]