はじめに
この記事では、共変量を所与とし、介入効果と結果変数とが無関係であるという条件付き独立の仮定(Conditional Independence Assumption;CIA)を扱います。
内容はJoshua D. Angrist & Jorn-steffen Pischke (2008)『Mostly harmless econometrics』Ch.3.2.1"The Conditional Independence Assumption"を参考にしています。
条件付き独立の仮定とは?
観察されたデータを元に因果推論を行うとき、潜在結果(potential outcome)という概念を用います。我々は大学進学が将来の年収に与える影響について関心があるとしましょう。
大学に進学するなら、進学しないならです。は大学に行った場合の年収、は行かなかった場合の年収です。ここでとの差は「同じさんが大学に行った場合と行かなかった場合の年収の差」なので因果効果(causal effect)を表します。実際に観察したデータである年収を、この因果効果を使って表してみましょう。
ここで重要なのは、我々はとのどちらか一方しか観察できないということです。同じさんは大学に行く・行かないの両方を選択できません*1。潜在結果とは、潜在的にが存在しますが、の取りうる値によって実現する値(つまり我々が観察できる値)が異なる、という意味です。
因果関係を検証したい、けれどとの片方しか観察できない。ならば同じ個人さんではなく、すべての人を対象にして「すべての人が大学に行った場合の平均年収とすべての人が大学に行かなかった場合の平均年収の差」もしくは「大学に行った人グループにおける平均的な因果効果」つまりを代わりに見ることで因果関係を検証しよう、という発想に至る訳です。
一般的には、大学に行ったグループと行かなかったグループの平均的な年収を比較しても、因果関係は分かりませんし、直観的には「当然大卒者の方が年収は高いだろう」とも考えれます。この場合、正のセレクションバイアスが発生していることになります。つまり「元々高い年収を得る素質や特徴を持つ人が大学に進学する」というバイアスが存在するため、大学に行った人と大学に行かなかった人の平均年収を比べても、因果効果を観察することにはならないのです。
前置きが長くなってしまいましたが、このセレクションバイアスに対処するための仮定が、条件付き独立の仮定(Conditional Independence Assupmtion; CIA)です。
CIAは、年収に影響を与える素質や特徴などの共変量を所与(条件)として、結果変数と介入効果が独立であることをサポートします。
(3)の両辺を共変量で条件付けてみましょう。
(3)の第3項のセレクションバイアスがなくなり、因果効果を求める式となりました。ただし(5)はによって因果効果の値が異なります。共変量をさんの家族の年収だとすると、元々実家が裕福な人が大学に行った場合の年収に与える効果は○○、元々実家が貧乏な人の効果は××と、共変量次第で因果効果の大きさが変わってしまうため、一般的にどれくらい効果があるのかが分かりません。そこで因果効果を表す包括的な1つの指標を考えます。(5)に繰り返し期待値の法則をあてはめてみましょう。
これを母集団の条件なしの平均因果効果(unconditional average causal effect)と呼びます。同様に、大学に行ったグループの平均的な因果効果は、
となり、これは大学に行った人において、大卒が年収に与えた平均的な効果を表します。(6)はの周辺分布を用いて、すべてのによる効果の平均をとることで求められます。同様に(7)は大学に行ったグループにおけるの周辺分布を使い、同じくによる効果の平均をとることで求まります。
これらの変数はどちらも実証分析上、マッチング推計量(matching estimator)と呼びます。マッチングの話題は別の節で扱います。簡単に言うと、大学進学という選択以外は同じようなを持つさんたちのグループを2つ作り、それら別々のさん同士をマッチングさせ、平均的な年収の違いを比較する、という考え方です。
ここまで介入効果を、大学進学の有無の例を使って説明しましたが、二項変数である必要はありません。例えば教育年数をとし、教育年数の違いが年収に与える影響についても同様に考えることができます。このとき潜在結果変数をと定義し、(6)、(7)式と同じように因果効果を表すことができます。例えば高卒と大卒の違いが収入に与える影響を考えましょう。小学校から数えると高卒までの教育年数は、大卒まではですね。条件なしの平均因果効果、大卒グループの平均因果効果はそれぞれ次のように表すことができます。
回帰分析で因果推論ができる条件
条件付き独立の仮定(CIA)が成立することにより、回帰分析を用いて因果推論を行うことが可能になります。この「CIA→回帰分析→因果推論」のつながりについて、以下2通りの考え方を紹介します。ここからは潜在結果を使って考えましょう。
まずは が線形かつすべてのさんにとって関数形が同じ、というケースです。このとき線形回帰はを推定するモデルそのものになりますね。ですが現実に即すとこの仮定は厳しいものです。個々人にとっての形は異なるでしょうし、ましてや一律線形とも考えにくいです。なのでもう1つは、さんによっての形は異なり、非線形なも認めるケースです。このとき回帰分析は個々人の教育年数による年収の違い、の重みづけ平均の推定値を導出することになります。
まずは1ケース目を考えましょう。
(7)は線形な因果モデルです。右辺のにはが含まれていません。何を言いたいかと言うと、(8)では実際にさんが受けた教育年数についての情報は含んでおらず、は仮に年の教育を受けた場合はどれぐらい稼ぐのかを表しています。観察されたデータを(7)に代入してみましょう。
左辺も観察された年収(実現値)となっております。ここで重要なのがは潜在結果変数や誤差項のとも相関している可能性があるということです。観察された共変量についてCIAが成立するとしましょう。(9)の誤差項をの線形結合と平均独立(mean independence)な(つまりを満たす)に区別します。
はを満たす回帰係数ベクトルです。
CIAの仮定を用いることで、潜在結果変数の条件付き期待値は以下のように表せます。
これで誤差項が説明変数と独立しており、が検証したい因果効果を表す線形モデルを導出することができました。
最後に重要な点ですが、(11)を導出するまでの一連の流れにおいては、こそがと(同様にと)の相関を生んでいる、という仮定があることです。このを明示的に回帰モデルに登場させることで、CIAの下と独立な介入効果が潜在変数に与える影響を「因果効果」として捉えることができるのです。セレクションバイアスもそもそも共変量ありきで発生する、と考えられます。
終わりに・感想
次回は適切な共変量をモデルに含めないことで生じる欠落変数バイアス(omitted variable bias)を扱います。読んでいただきありがとうございました。
*1:本文ではロバート・フロストの詩「選ばれざる道(Road Not Taken)」が引用されています。めちゃくちゃ良かったです。