傾向スコア・マッチングのソースを表示

'''傾向スコア・マッチング'''（プロペンシティ・スコア・マッチング、propensity score matching、PSM）は、観察データの[[統計学|統計]]分析の分野において、治療を受けることを予測する共変量を考慮して、処置（treatment）、方針、その他介入の効果を推定しようとする[[マッチング (統計学)|マッチング]]手法。処置を受けた人々と受けなかった人々の結果を単純に比較して治療効果を推定すると[[交絡]]変数によるバイアス（[[偏り]]）が発生する。このバイアスを軽減するための手法が傾向スコア・マッチングであり、[[1983年]]、[[ポール・ローゼンバウム]]と[[ドナルド・ルービン]]が発表した <ref name="Rosenbaum 1983 41–55">{{Cite journal|last=Rosenbaum|first=Paul R.|last2=Rubin|first2=Donald B.|year=1983|title=The Central Role of the Propensity Score in Observational Studies for Causal Effects|journal=[[Biometrika]]|volume=70|issue=1|pages=41–55|DOI=10.1093/biomet/70.1.41}}</ref>。

処置群とコントロール群（非処置群）の処置結果（[[平均処置効果]]など）の違いは、処置そのものではなく処置を予測する要因によって引き起こされる可能性があり、その場合はバイアスが発生する。[[ランダム化比較試験]]では、[[無作為]]割り付けによってバイアスなく処置効果を推定することができる。無作為割り付けによって、各共変量のバランスが取れることを[[大数の法則]]が保証する。残念ながら、観察研究の多くで、処置の無作為割り付けはなされていない。マッチングでは、観察された共変量が同じくらいの標本を処置群とコントロール群のそれぞれから抽出することにより、割り付けバイアスを減らして、無作為割り付けに近いものにする。

たとえば、[[たばこ病|喫煙]]の影響を知りたい場合を考える。人々を喫煙群に無作為に割り付けることは非倫理的であるため、観察研究が必要である。喫煙群と非喫煙群とを単純に比較することによって処置効果を推定すると、喫煙率に影響する要因（性別や年齢など）によるバイアスが生じる。傾向スコア・マッチングでは、処置群とコントロール群の制御変数（この例では性別や年齢など）を同じくらいにすることによって、これらのバイアスを制御することを目指す。

== 概要 ==
傾向スコア・マッチングは、以下の様な非実験的設定における[[帰納|因果推論]]および選択バイアスに対して使用される。

# コントロール群には、処置群の被験者と同等といえる被験者がほとんどいない
# 高次元の特徴量を元に比較する必要があるため、処置群の被験者とよく似たコントロール群の被験者を選択することが難しい

通常のマッチングでは、治療群と対照群を区別する単一の特性が照合される。ただし、2つのグループに実質的な重複がない場合は、かなりの[[誤差]]が発生する可能性がある。たとえば、コントロール群の最悪のケースと処置群の最良のケースとの比較になってしまうと、[[平均への回帰]]によって、実際によりも良く見えたり悪く見えたりする可能性がある。

傾向スコア・マッチングでは割り付けの予測確率（たとえば、処置群とコントロール群のどちらに割り付けられるかの確率）を採用する。この予測確率は、観測された予測因子に基づいた[[ロジスティック回帰]]によって取得されることが多く、反実仮想群を作成することができる。傾向スコアは、単独または他のマッチング変数または共変量と一緒に、マッチングに使ったり共変量として扱ったりすることができる

== 一般的な手順 ==
1.  [[ロジスティック回帰]]を実行する

* 従属変数：治療群では <math>Z = 1</math>、コントロール群（非治療群）では <math>Z = 0</math> となる。
* 適切な交絡因子（治療と結果の両方に関連すると仮定された変数）を選択する。
* 傾向スコアの[[推定量]]（予測確率 <math>p </math>）を得る。

2. 傾向スコアで層別化した上で、層内の治療群とコントロール群との間で共変量のバランスが取れていることを確認する

* 標準化群間差（standardized difference）またはグラフで分布を調べる

3. 次のいずれかの方法を用いて、傾向スコアに基づいて治療群の各被験者をコントロール群の被験者 1 人（または複数）とマッチングする

* [[最近傍探索]]
* キャリパーマッチング
* [[マハラノビス距離]]マッチング
* [[層化抽出法]]
* 差分マッチング
* 完全一致

4. マッチングした被験者において、処置群とコントロール群との間で共変量のバランスが取れていることを確認する

5. 新しいサンプルに基づいて多変量解析を行う

* 処置群の被験者のマッチング対象をコントロール群の被験者から複数選ぶ場合は、独立性を仮定しない手法で解析する。すなわち、通常の最小二乗法ではなく、加重最小二乗法を使用する。

== 正式な定義 ==

=== 基本設定 ===
基本的なケースとして、2 つの処置（処置群が 1、コントロール群が 0 ）を <math> N  </math>人の被験者に割り付け、[[独立同分布]]（i.i.d.）を仮定する <ref name="Rosenbaum 1983 41–55">{{Cite journal|last=Rosenbaum|first=Paul R.|last2=Rubin|first2=Donald B.|year=1983|title=The Central Role of the Propensity Score in Observational Studies for Causal Effects|journal=[[Biometrika]]|volume=70|issue=1|pages=41–55|DOI=10.1093/biomet/70.1.41}}</ref>。

<math> i  </math>  番目の被験者（被験者 <math> i  </math> とする）の処置への反応 <math>r_{1i}</math> とコントロールへの反応 <math>r_{0i}</math> とを考えて、[[平均処置効果|平均治療効果]] <math>E[r_1]-E[r_0]</math> を推定する。

変数 <math>Z_i</math> を用いて、被験者 <math> i  </math> が処置群に割り付けられた場合は <math> Z_i = 1 </math>、コントロール群に割り付けられた場合は <math> Z_i = 0 </math> と表現する。<math>X_i</math> を被験者 <math>i </math> の処置の割り当て前に観測された値を表すベクトルとする。処置の割り当てを決定する場合に参考にする変数の中には、<math>X_i</math> に含まれないもののあるかもしれない。番号付け（<math> i = 1, 2, 3, \dots, N </math>）の値そのものには、<math>X_i</math> の内容に関する情報は含まれていないものとする。以降も個々の被検者について議論するが、<math>i </math>  の記載は省略する。

=== 「強く無視可能」な治療の割り当て (SITA) ===
ある被験者が共変量 <math>X </math>を持つ（すなわち、条件付き交絡がない）ものとする。
治療群およびコントロール群の'''潜在アウトカム'''を、それぞれ、<math>r_0 </math>、<math>r_1</math> とする。
潜在アウトカムが背景変数 <math>X</math> を条件とする治療 <math>Z</math> とは[[独立 (確率論)|独立]]である場合、治療の割り当ては「'''強く無視可能'''である」とされる。
Strongly ignorable treatment assignment (SITA) は、以下のように簡潔に記載することができる。

: <math> r_0, r_1 \perp Z \,|\, X </math>

ここで、<math>\perp</math> は[[独立 (確率論)|統計的独立性]]を示す <ref name="Rosenbaum 1983 41–55">{{Cite journal|last=Rosenbaum|first=Paul R.|last2=Rubin|first2=Donald B.|year=1983|title=The Central Role of the Propensity Score in Observational Studies for Causal Effects|journal=[[Biometrika]]|volume=70|issue=1|pages=41–55|DOI=10.1093/biomet/70.1.41}}</ref>。

=== バランシングスコア ===
'''バランシングスコア'''（balancing score）<math>b(X)</math> は観測された共変量 <math>X </math> の関数であり、<math>b(X) </math> に対する <math>X </math> の[[条件付き確率|条件付き分布]]は、処置群の被験者（<math>Z = 1</math>）とコントロール群の被験者（<math>Z = 0</math>）とで等しい。

: <math> Z \perp X \,|\, b(X)</math>

最もシンプルな関数は <math>b(X) = X</math> である。

=== 傾向スコア ===
'''傾向スコア'''は、所与の観察値に応じて、被験者（人、教室、学校など）が特定の処置に割り当てられる[[確率]]を示す。傾向スコアは、これらの共変量に基づき処置群とコントロール群とを均質化して、[[観測選択効果|選択バイアス]]を減らすために用いられる。

2 値変数の処置[[指示関数|指標]] <math> Z </math>、応答変数 <math> r </math>、およびバックグラウンドで観測された共変量 ''<math>X </math>'' があるとする。傾向スコアは、バックグラウンド変数に対する処置の[[条件付き確率]]として定義される。

: <math>e(x) \ \stackrel{\mathrm{def}}{=}\ \Pr(Z=1 \mid X=x)</math>

=== 主な定理 ===
[[1983年]]、[[ポール・ローゼンバウム]]と[[ドナルド・ルービン]]が下記内容を示した <ref name="Rosenbaum 1983 41–55">{{Cite journal|last=Rosenbaum|first=Paul R.|last2=Rubin|first2=Donald B.|year=1983|title=The Central Role of the Propensity Score in Observational Studies for Causal Effects|journal=[[Biometrika]]|volume=70|issue=1|pages=41–55|DOI=10.1093/biomet/70.1.41}}</ref>

* 傾向スコア <math>e(x)</math> はバランシングスコアである。
* 関数 <math>f</math> を用いて <math>e(X) = f(b(X))</math> と表されるような、傾向スコア <math>e(X)</math> よりも細かい（finer）スコア <math>b(X)</math> は、バランシングスコアである。
** 最も粗い（coarsest）バランシングスコア関数は傾向スコアである。（多次元オブジェクト <math>X_i</math> を 1 次元に変換する）
** 最も細かい（finest）バランシングスコア関数は <math>b(X) = X</math> である。
* 任意の <math>X</math> に対し、処置の割り付けが'''強く無視可能'''な場合は、次のようになる。
** 任意のバランシングスコア関数に対して、'''強く無視可能'''である。具体的には、任意の傾向スコアに対して <math> (r_0, r_1) \perp Z \mid e(X)</math>
** バランシングスコアの任意の値について、バランシングスコアの値が同じである被験者に基づく、標本中の処置群とコントロール群の平均の差 <math>\bar{r}_1-\bar{r}_0</math> は、[[平均処置効果]]の不偏推定量 <math>E[r_1]-E[r_0]</math> として機能する。
* バランシングスコアの標本推定量を使用すると、<math>X</math> に関する標本の均衡が得られる。

=== 十分統計量との関係 ===
<math>Z</math> の値を <math>X</math> の分布に影響を与える母集団のパラメータと考えると、バランススコアは <math>Z</math> の[[十分統計量]]として機能する。
さらに、上記の定理は、<math>Z</math> を <math>X</math> のパラメーターとして考える場合、傾向スコアは最小[[十分統計量|十分統計量で]]あることを示している。
最後に、<math>X</math> に対して処置割り付け <math>Z</math> が強く無視可能な場合、傾向スコアは、同時分布 <math>(r_0, r_1)</math> の最小[[十分統計量]]となる。

=== 交絡変数の存在を検出するためのグラフィカルテスト ===
[[ジューディア・パール]]は、交絡変数の存在を検出する、バックドア基準と呼ばれる単純なグラフィカルテストが存在することを示した。治療の効果を推定するには、回帰する際に交絡変数として補正するか、交絡変数をマッチングすることによって、グラフ内の全てのバックドア経路をブロックすることが必要である <ref name="pearl">{{Cite book|last=Pearl|first=J.|year=2000|title=Causality: Models, Reasoning, and Inference|url=https://archive.org/details/causalitymodelsr0000pear|location=New York|publisher=Cambridge University Press|isbn=978-0-521-77362-1}}</ref>。

== 短所 ==
傾向スコア・マッチングはモデルの「不均衡、非効率、モデル依存、バイアス」を増加させる場合があることが示され、Gary King によって他のマッチング方法に比して推奨されない、とされた <ref>{{Cite journal|last=King|first=Gary|last2=Nielsen|first2=Richard|date=2019-05-07|title=Why Propensity Scores Should Not Be Used for Matching|journal=Political Analysis|volume=27|issue=4|pages=435–454|DOI=10.1017/pan.2019.11|ISSN=1047-1987}} | [https://gking.harvard.edu/files/gking/files/psnot.pdf link to the full article] (from the author's homepage)</ref>。マッチングの背後にある洞察は有効であっても、他のマッチング法が必要になる。また、傾向スコアには重み付けや二重ロバスト推定などマッチング以外の用途もある。

他のマッチング法と同様、傾向スコア・マッチングは観測データから[[平均処置効果]]を推定する。傾向スコア・マッチングが紹介された当初は、単一のスコアに対して共変量の線形結合を使用することにより、多数の観測値を失うことなく、多数の共変量について処置群とコントロール群を均衡させることができることが利点として挙げられた。処置群とコントロール群とが多くのの共変量において均衡している場合、「[[次元の呪い]]」を克服するために膨大な観測データが必要になる。すなわち、均衡させる共変量の項目を増やせば増やすほど、必要となる観測データが[[等比数列|幾何級数的]]に増大する。

傾向スコア・マッチングの欠点の1つは、観測された（および観測可能な）共変量のみを考慮し、潜在的な変数は考慮しないことである。処置への割り付けとアウトカムに影響を与えるが観察できない要因は、マッチングでは考慮できない <ref>{{Cite journal|year=2014|title=Methods for Constructing and Assessing Propensity Scores|journal=Health Services Research|volume=49|issue=5|pages=1701–20|DOI=10.1111/1475-6773.12182|PMID=24779867|PMC=4213057}}</ref>。観測データのみをコントロールするため、マッチングしても潜在変数による隠れたバイアスが残る可能性がある <ref>{{Cite book|last=Shadish|first=W. R.|last2=Cook|first2=T. D.|last3=Campbell|first3=D. T.|year=2002|title=Experimental and Quasi-experimental Designs for Generalized Causal Inference|location=Boston|publisher=Houghton Mifflin|isbn=978-0-395-61556-0}}</ref>。

さらに、傾向スコア・マッチングでは、処置群とコントロール群とが十分に重なり合った、膨大な数の標本が必要である。

[[ジューディア・パール]]は、マッチングに関する一般的な懸念事項として、観測された変数に基づいてマッチングすることで、観測されていない交絡因子によるバイアスが表面化する可能性があると主張した。さらに、パールは、バイアスの低減は、処置、アウトカム、観察された共変量と観察されなかった共変量の間の定性的な因果関係をモデル化することによってのみ（漸近的に）保証できると主張した <ref name="pearl:ch11-3-5">{{Cite book|last=Pearl|first=J.|chapter=Understanding propensity scores|title=Causality: Models, Reasoning, and Inference|location=New York|publisher=Cambridge University Press|edition=Second|year=2009|isbn=978-0-521-89560-6}}</ref>。交絡は、実験者が独立変数と従属変数の間に観察された関係について、本来の因果関係以外の原因をコントロールできない場合に発生する。これをコントロールするためには、パールの「バックドア基準」を満たす必要がある <ref name="pearl">{{Cite book|last=Pearl|first=J.|year=2000|title=Causality: Models, Reasoning, and Inference|url=https://archive.org/details/causalitymodelsr0000pear|location=New York|publisher=Cambridge University Press|isbn=978-0-521-77362-1}}</ref>。

== 統計パッケージでの実装 ==
* [[R言語|R]] – 傾向スコアマッチングは <code>MatchIt</code> パッケージで利用することができるほか <ref>{{Cite journal|last=Ho|first=Daniel|last2=Imai|first2=Kosuke|last3=King|first3=Gary|last4=Stuart|first4=Elizabeth|year=2007|title=Matching as Nonparametric Preprocessing for Reducing Model Dependence in Parametric Causal Inference|journal=[[Political Analysis (journal)|Political Analysis]]|volume=15|issue=3|pages=199–236|DOI=10.1093/pan/mpl013}}</ref> <ref>{{Cite web|title=MatchIt: Nonparametric Preprocessing for Parametric Causal Inference|website=R Project|date=|url=https://cran.r-project.org/package=MatchIt|accessdate=2020-12-30}}</ref>、手動で実装することもできる <ref>{{Cite book|first=Andrew|last=Gelman|first2=Jennifer|last2=Hill|title=Data Analysis Using Regression and Multilevel/Hierarchical Models|location=New York|publisher=Cambridge University Press|year=2007|isbn=978-0-521-68689-1|pages=206–212|url=https://books.google.com/books?id=lV3DIdV0F9AC&pg=PA206}}</ref>。
* SAS – PSMatchプロシージャおよび  <code>OneToManyMTCH</code> マクロは、傾向スコアに基づいてマッチングする <ref>{{Cite web|first=Lori|author=Parsons|title=Performing a 1:N Case-Control Match on Propensity Score|publisher=SAS Institute|url=http://www2.sas.com/proceedings/sugi29/165-29.pdf|accessdate=June 10, 2016}}</ref>。
* [[Stata]] – ユーザー作成の <code>psmatch2</code> コマンドなどが傾向スコア・マッチングを実装している <ref>[http://fmwww.bc.edu/RePEc/usug2001/psmatch.pdf Implementing Propensity Score Matching Estimators with STATA]. Lecture notes 2001</ref> <ref>{{Cite journal|last=Leuven|first=E.|last2=Sianesi|first2=B.|date=2003|title=PSMATCH2: Stata module to perform full Mahalanobis and propensity score matching, common support graphing, and covariate imbalance testing|url=http://ideas.repec.org/c/boc/bocode/s432001.html}}</ref> 。Stataバージョン13以降では、組み込みの <code>teffects psmatch</code> コマンドも提供されている <ref>{{Cite web|title=teffects psmatch — Propensity-score matching|website=Stata Manual|url=https://www.stata.com/manuals15/teteffectspsmatch.pdf|accessdate=2020-12-30}}</ref>。
* [[SPSS]] – IBM SPSS Statistics の「傾向スコアによる一致」で傾向スコア・マッチングのダイアログボックスが開き、「適合の許容」、「完全一致を優先」、「実行パフォーマンスを最大化」、「一致の抽出時にケース順序をランダム化」「乱数のシード」などの設定が可能である。FUZZY Python プロシージャは、拡張機能ダイアログボックスからソフトウェアの拡張機能として追加することもできる。このプロシージャは、指定されたキー変数のセットに基づいて、コントロール群からの無作為抽出を利用して、処置群とコントロール群とをマッチングする。 FUZZYコマンドは、完全一致とあいまい一致をサポートする。

== 関連項目 ==
* [[ルービン因果モデル]]
* 無視可能性
* ヘックマン修正
* [[マッチング (統計学)]]
*[[統計的因果推論]]

== 書籍 ==

* 星野 崇宏(2009)『調査観察データの統計科学―因果推論・選択バイアス・データ融合 (シリーズ確率と情報の科学)』岩波書店

== 脚注 ==
{{Reflist|30em}}

== 参考文献 ==
* {{Cite journal|last=Abadie|first=Alberto|author-link=Alberto Abadie|last2=Imbens|first2=Guido W.|year=2006|title=Large Sample Properties of Matching Estimators for Average Treatment Effects|journal=[[Econometrica]]|volume=74|issue=1|pages=235–267|DOI=10.1111/j.1468-0262.2006.00655.x}}
* {{Cite book|last=Leite|first=Walter L.|year=2017|title=Practical Propensity Score Methods using R|location=Washington, DC|publisher=Sage Publications|isbn=978-1-4522-8888-8}}

== 外部リンク ==
* [https://cran.r-project.org/web/packages/MatchIt/ MatchIt: Nonparametric Preprocessing for Parametric Causal Inference] - [[CRAN]]
* [https://cran.r-project.org/web/packages/Matching/ Matching: Multivariate and Propensity Score Matching with Balance Optimization] - [[CRAN]]


{{統計学}}

{{デフォルトソート:けいこうすこあまつちんく}}
[[Category:因果推論]]
[[Category:統計学]]
[[Category:疫学]]
[[Category:回帰分析]]