傾向スコア・マッチング

傾向スコア・マッチング（プロペンシティ・スコア・マッチング、propensity score matching、PSM）は、観察データの統計分析の分野において、治療を受けることを予測する共変量を考慮して、処置（treatment）、方針、その他介入の効果を推定しようとするマッチング手法。処置を受けた人々と受けなかった人々の結果を単純に比較して治療効果を推定すると交絡変数によるバイアス（偏り）が発生する。このバイアスを軽減するための手法が傾向スコア・マッチングであり、1983年、ポール・ローゼンバウムとドナルド・ルービンが発表した 。
処置群とコントロール群（非処置群）の処置結果（平均処置効果など）の違いは、処置そのものではなく処置を予測する要因によって引き起こされる可能性があり、その場合はバイアスが発生する。ランダム化比較試験では、無作為割り付けによってバイアスなく処置効果を推定することができる。無作為割り付けによって、各共変量のバランスが取れることを大数の法則が保証する。残念ながら、観察研究の多くで、処置の無作為割り付けはなされていない。マッチングでは、観察された共変量が同じくらいの標本を処置群とコントロール群のそれぞれから抽出することにより、割り付けバイアスを減らして、無作為割り付けに近いものにする。
たとえば、喫煙の影響を知りたい場合を考える。人々を喫煙群に無作為に割り付けることは非倫理的であるため、観察研究が必要である。喫煙群と非喫煙群とを単純に比較することによって処置効果を推定すると、喫煙率に影響する要因（性別や年齢など）によるバイアスが生じる。傾向スコア・マッチングでは、処置群とコントロール群の制御変数（この例では性別や年齢など）を同じくらいにすることによって、これらのバイアスを制御することを目指す。
概要.
傾向スコア・マッチングは、以下の様な非実験的設定における因果推論および選択バイアスに対して使用される。
通常のマッチングでは、治療群と対照群を区別する単一の特性が照合される。ただし、2つのグループに実質的な重複がない場合は、かなりの誤差が発生する可能性がある。たとえば、コントロール群の最悪のケースと処置群の最良のケースとの比較になってしまうと、平均への回帰によって、実際によりも良く見えたり悪く見えたりする可能性がある。
傾向スコア・マッチングでは割り付けの予測確率（たとえば、処置群とコントロール群のどちらに割り付けられるかの確率）を採用する。この予測確率は、観測された予測因子に基づいたロジスティック回帰によって取得されることが多く、反実仮想群を作成することができる。傾向スコアは、単独または他のマッチング変数または共変量と一緒に、マッチングに使ったり共変量として扱ったりすることができる
一般的な手順.
1. ロジスティック回帰を実行する
2. 傾向スコアで層別化した上で、層内の治療群とコントロール群との間で共変量のバランスが取れていることを確認する
3. 次のいずれかの方法を用いて、傾向スコアに基づいて治療群の各被験者をコントロール群の被験者 1 人（または複数）とマッチングする
4. マッチングした被験者において、処置群とコントロール群との間で共変量のバランスが取れていることを確認する
5. 新しいサンプルに基づいて多変量解析を行う
正式な定義.
基本設定.
基本的なケースとして、2 つの処置（処置群が 1、コントロール群が 0 ）を formula_4人の被験者に割り付け、独立同分布（i.i.d.）を仮定する 。
formula_5 番目の被験者（被験者 formula_5 とする）の処置への反応 formula_7 とコントロールへの反応 formula_8 とを考えて、平均治療効果 formula_9 を推定する。
変数 formula_10 を用いて、被験者 formula_5 が処置群に割り付けられた場合は formula_12、コントロール群に割り付けられた場合は formula_13 と表現する。formula_14 を被験者 formula_15 の処置の割り当て前に観測された値を表すベクトルとする。処置の割り当てを決定する場合に参考にする変数の中には、formula_14 に含まれないもののあるかもしれない。番号付け（formula_17）の値そのものには、formula_14 の内容に関する情報は含まれていないものとする。以降も個々の被検者について議論するが、formula_15 の記載は省略する。
「強く無視可能」な治療の割り当て (SITA).
ある被験者が共変量 formula_20を持つ（すなわち、条件付き交絡がない）ものとする。
治療群およびコントロール群の潜在アウトカムを、それぞれ、formula_21、formula_22 とする。
潜在アウトカムが背景変数 formula_23 を条件とする治療 formula_24 とは独立である場合、治療の割り当ては「強く無視可能である」とされる。
Strongly ignorable treatment assignment (SITA) は、以下のように簡潔に記載することができる。
ここで、formula_26 は統計的独立性を示す 。
バランシングスコア.
バランシングスコア（balancing score）formula_27 は観測された共変量 formula_20 の関数であり、formula_29 に対する formula_20 の条件付き分布は、処置群の被験者（formula_1）とコントロール群の被験者（formula_2）とで等しい。
最もシンプルな関数は formula_34 である。
傾向スコア.
傾向スコアは、所与の観察値に応じて、被験者（人、教室、学校など）が特定の処置に割り当てられる確率を示す。傾向スコアは、これらの共変量に基づき処置群とコントロール群とを均質化して、選択バイアスを減らすために用いられる。
2 値変数の処置指標 formula_35、応答変数 formula_36、およびバックグラウンドで観測された共変量 "formula_20" があるとする。傾向スコアは、バックグラウンド変数に対する処置の条件付き確率として定義される。
主な定理.
1983年、ポール・ローゼンバウムとドナルド・ルービンが下記内容を示した 
十分統計量との関係.
formula_24 の値を formula_23 の分布に影響を与える母集団のパラメータと考えると、バランススコアは formula_24 の十分統計量として機能する。
さらに、上記の定理は、formula_24 を formula_23 のパラメーターとして考える場合、傾向スコアは最小十分統計量であることを示している。
最後に、formula_23 に対して処置割り付け formula_24 が強く無視可能な場合、傾向スコアは、同時分布 formula_58 の最小十分統計量となる。
交絡変数の存在を検出するためのグラフィカルテスト.
ジューディア・パールは、交絡変数の存在を検出する、バックドア基準と呼ばれる単純なグラフィカルテストが存在することを示した。治療の効果を推定するには、回帰する際に交絡変数として補正するか、交絡変数をマッチングすることによって、グラフ内の全てのバックドア経路をブロックすることが必要である 。
短所.
傾向スコア・マッチングはモデルの「不均衡、非効率、モデル依存、バイアス」を増加させる場合があることが示され、Gary King によって他のマッチング方法に比して推奨されない、とされた 。マッチングの背後にある洞察は有効であっても、他のマッチング法が必要になる。また、傾向スコアには重み付けや二重ロバスト推定などマッチング以外の用途もある。
他のマッチング法と同様、傾向スコア・マッチングは観測データから平均処置効果を推定する。傾向スコア・マッチングが紹介された当初は、単一のスコアに対して共変量の線形結合を使用することにより、多数の観測値を失うことなく、多数の共変量について処置群とコントロール群を均衡させることができることが利点として挙げられた。処置群とコントロール群とが多くのの共変量において均衡している場合、「次元の呪い」を克服するために膨大な観測データが必要になる。すなわち、均衡させる共変量の項目を増やせば増やすほど、必要となる観測データが幾何級数的に増大する。
傾向スコア・マッチングの欠点の1つは、観測された（および観測可能な）共変量のみを考慮し、潜在的な変数は考慮しないことである。処置への割り付けとアウトカムに影響を与えるが観察できない要因は、マッチングでは考慮できない 。観測データのみをコントロールするため、マッチングしても潜在変数による隠れたバイアスが残る可能性がある 。
さらに、傾向スコア・マッチングでは、処置群とコントロール群とが十分に重なり合った、膨大な数の標本が必要である。
ジューディア・パールは、マッチングに関する一般的な懸念事項として、観測された変数に基づいてマッチングすることで、観測されていない交絡因子によるバイアスが表面化する可能性があると主張した。さらに、パールは、バイアスの低減は、処置、アウトカム、観察された共変量と観察されなかった共変量の間の定性的な因果関係をモデル化することによってのみ（漸近的に）保証できると主張した 。交絡は、実験者が独立変数と従属変数の間に観察された関係について、本来の因果関係以外の原因をコントロールできない場合に発生する。これをコントロールするためには、パールの「バックドア基準」を満たす必要がある 。