逆確率重み付けのソースを表示

'''逆確率重み付け'''（ぎゃくかくりつおもみづけ、Inverse probability weighting、IPW）は、サンプリング母集団（データ収集に用いた母集団）とは異なる、擬似的な母集団に標準化された統計量を算出する統計手法。サンプリング母集団とターゲット母集団（推論の対象となる母集団）が異なることは多い。コストや時間、倫理的な問題などのため、ターゲット母集団から直接サンプリングすることが難しい場合もある。[[層化抽出法]]などの代替デザイン戦略は一つの解決策である。重み付けを正しく適用することで、推定量のバイアスを低減することができる。

最初期の重み付け推定量の1つに、平均値の[[Horvitz–Thompson推定量]]がある<ref>{{Cite journal|last=Horvitz|first=D. G.|last2=Thompson|first2=D. J.|year=1952|title=A generalization of sampling without replacement from a finite universe|journal=[[Journal of the American Statistical Association]]|volume=47|issue=260|pages=663–685|DOI=10.1080/01621459.1952.10483446}}</ref>。ターゲット母集団からサンプリング母集団として抽出されるサンプリング確率が既知の場合、この確率の逆数を観測値の重み付けに使用する。このアプローチは、さまざまなフレームワークのもと、統計学の多くの分野で広く用いられている。重み付き[[尤度]]、重み付き推定方程式、重み付き[[確率密度]]であり、これらから多くの統計量が派生する。これらの応用により、[[限界構造モデル]]、[[標準化死亡比]]、粗データや集合データに対する[[EMアルゴリズム]]など、他の統計および推定量の理論が体系化された。

逆確率重み付けは、欠測データのある被験者を一次分析に含めることができない場合に、欠測データを考慮するためにも用いられる。サンプリング確率（その因子が別の測定で測定される確率）の推定値があれば、逆確率重み付けを用いて、[[欠測データ]]が多いために過小評価されている被験者の重みを増加させることができる。

== 逆確率重み付け推定量（IPWE） ==
対象実験は実施できないがモデル化可能な観測データはある場合、逆確率重み付け推定量（Inverse Probability Weighted Estimator, IPWE）を用いることで因果関係を示すことができる。治療の無作為割付は仮定しておらず、母集団のすべての被験者に対して特定の治療を割り当てた場合の反実仮想アウトカム（潜在的アウトカム）を推定することが目標となる。

観測データ <math>\{\left( X_i, A_i, Y_i \right) \}_{i=1}^{n}</math> が[[独立同分布]]で未知の分布 <math>P</math> に従うと仮定する。ここで、

* <math>X \in \mathbb{R}^{p}</math>：共変量
* <math>A \in \{0, 1\}</math>：2つの可能な治療法
* <math>Y \in \mathbb{R}</math>：反応
* 治療の無作為割付は想定していない。

まずは、潜在アウトカム <math>Y^{*}(a)</math> 、すなわち1人の被験者に治療 <math>a</math> が割り当てられた場合のアウトカムを推定する。次に、母集団のすべての患者に治療 <math>a</math> が割り当てられた場合の平均アウトカム <math>\mu_a = \mathbb{E} \left( Y^*(a) \right)</math> を比較する。観測データ<math>\{\left( X_i,A_i,Y_i \right) \}_{i=1}^n</math> から <math>\mu_a</math> を推定したい。

=== 推定量の式 ===
:<math>\hat{\mu}_{a,n}^\mathrm{IPWE} = \frac{1}{n} \sum_{i=1}^n \, Y_{i} \frac{\mathbf{ 1}_{A_i=a}}{\hat{p}_n (A_i \mid X_{i})}</math>

==== IPWEの構築 ====

# <math>\mu_{a} = \mathbb{E} \left( \frac{Y \, \mathbf{1}_{A=a}}{p(a \mid x)} \right)</math> ここで、<math>p(a \mid x) = \frac{P(A=a, X=x)}{P(X=x)}</math>
# 任意のプロペンシティモデル（多くはロジスティック回帰モデル）を用いて <math>\hat{p}_{n}(a \mid x)</math> ないし <math>p(a \mid x)</math> を構築する
# <math>\hat{\mu}_{a,n}^\mathrm{IPWE} = \frac{\Sigma_{i=1}^{n} Y_i \, 1_{A_i=a}}{n\hat{p}_{n}(A_i \mid X_i)}</math>

各治療群の平均値を算出した後、[[t検定]]または[[ANOVA]]検定を用いて群平均間の差を判定し、治療効果の統計的有意性を判定することができる。

==== 仮定 ====
# 一貫性：<math>Y = Y^{*}(A)</math>
# 未測定交絡因子がない：<math>\{Y^{*}(0), Y^{*}(1)\} \perp A \mid X</math>
#* 治療の割り当ては、共変量データのみに基づいており、潜在的アウトカムとは無関係である。
# 正値性：すべての <math>a</math> および <math>x</math> に対して <math>P(A=a \mid X=x)>0 </math>

==== 制限事項 ====
推定される傾向が小さい場合、逆確率均等推定量（IPWE）は不安定になる可能性がある。
いずれかの治療割り当ての確率が小さい場合、ロジスティック回帰モデルはテール周辺で不安定になり、IPWEが安定しない。

== 拡張逆確率重み付け推定量（AIPWE） ==
その他の推定量として、拡張逆確率重み付き推定量（Augumented Inverse Probability Weighted Estimator AIPWE）がある。これは、回帰ベースの推定量と逆確率重み付き推定量の両方のプロパティを組み合わせたもので、傾向モデルまたは結果モデルのいずれかを正しく指定するだけで、両方を指定する必要がないという点で、「二重に堅牢な」方法である。IPWEを拡張して、変動を減らし、推定効率を向上させる。このモデルは、逆確率均等推定量（IPWE）と同じ仮定を保持している。 <ref>{{Cite journal|last=Cao|first=Weihua|last2=Tsiatis|first2=Anastasios A.|last3=Davidian|first3=Marie|year=2009|title=Improving efficiency and robustness of the doubly robust estimator for a population mean with incomplete data|journal=Biometrika|volume=96|issue=3|pages=723–734|DOI=10.1093/biomet/asp033|ISSN=0006-3444|PMID=20161511|PMC=2798744}}</ref>

=== 推定量の式 ===
:<math>\begin{align}
\hat{\mu}_{a,n}^\mathrm{AIPWE}
&= \frac{1}{n}\sum_{i=1}^n \left( \frac{Y_{i} \, 1_{A_i=a}}{\hat{p}_{n}(A_i \mid X_i)} -
\frac{1_{A_i=a}-\hat{p}_n(A_i \mid X_i)}{\hat{p}_n(A_i \mid X_i)}\hat{Q}_n(X_i,a) \right) \\
&= \frac{1}{n}\sum_{i=1}^n \left( \hat{Q}_n(X_i,a) \right) + \frac{1}{n} \sum_{i=1}^n\frac{1_{A_i=a}}{\hat{p}_n(A_i \mid X_i)} \left( Y_i - \hat{Q}_n(X_i,a) \right)
\end{align}
</math>

ただし、

# <math>1_{A_{i}=a}</math> は、被験者 i が治療群 a に属するか否かを示す[[指示関数]]である。
# 回帰推定量 <math>\hat{Q}_n(x,a)</math> を構築し、共変量<math>X</math> と治療 <math>A</math> に基づいて被験者 <math>i</math> におけるアウトカム <math>Y</math> を予測する。たとえば、通常の最小二乗回帰を使用する。
# プロペンシティ推定値 <math>\hat{p}_n(A_i \mid X_i)</math> を求める。たとえば、[[ロジスティック回帰]]を使用する。
# AIPWEとして組み合わせて <math>\hat{\mu}_{a,n}^\mathrm{AIPWE}</math> を得る。

=== 解釈と「二重ロバスト性」 ===
式を並べ替えると、根本的なアイデアが明らかになる。推定量は、モデルを用いて予測されたアウトカムの平均値（下記）に基づいている。

: <math>\frac{1}{n}\sum_{i=1}^n\Biggl(\hat{Q}_n(X_i,a)\Biggr)</math>

しかし、モデルにバイアスがある場合、モデルの残差は0付近にはならない。モデル Q の平均残差の項を追加することにより、この潜在的なバイアスを修正できる。

: <math>\frac{1}{n} \sum_{i=1}^n \frac{1_{A_i=a}}{\hat{p}_n(A_i \mid X_i)} \left( Y_i - \hat{Q}_n(X_i,a) \right)</math>

Yの値が欠落しているため、各残差の相対的な重要性を膨らませるために重みを付ける（これらの重みは、各被験者の観測値を見る逆の傾向、別名確率に基づく）（ <ref name="kang2007">Kang, Joseph DY, and Joseph L. Schafer. </ref>10ページを参照）。

「二重ロバスト性」は、推定量が不偏であるためには、 <math>\hat{Q}_n(X_i,a)</math> および <math>\hat{p}_{n}(A_i \mid X_i)</math> という 2つのモデルのうち、いずれかが正しく規定されていれば十分であるという事実に由来する。これは、アウトカムモデルが適切に規定されていれば、その残差は（重み付けに関係なく）0付近になるためである一方、モデルが不偏でなくても、重み付けモデルが適切に規定されている場合、そのバイアスは重み付け平均残差によって適切に推定（および補正）される<ref name="kang2007">Kang, Joseph DY, and Joseph L. Schafer. </ref> <ref>Kim, Jae Kwang, and David Haziza. </ref> <ref>Seaman, Shaun R., and Stijn Vansteelandt. </ref>。

二重ロバスト推定量のバイアスは'''2次バイアス'''と呼ばれ、<math>\frac{1}{\hat{p}_n(A_i \mid X_i)} - \frac{1}{p_n(A_i \mid X_i)}</math> と<math>\hat{Q}_n(X_i,a) - Q_n(X_i,a)</math> の積に依存する。この性質により、「十分に大きい」サンプルサイズがある場合、（パラメトリックモデルの代わりに）[[機械学習]]推定器を使用して、二重ロバスト推定量の全体的なバイアスを下げることができる<ref>Hernán, Miguel A., and James M. Robins. </ref>。

== 関連項目 ==
* [[傾向スコア・マッチング]]
* [[限界構造モデル]]

== 脚注 ==
{{Reflist}}

== 外部リンク ==
* [https://cran.r-project.org/web/packages/ipw/index.html ipw: Estimate Inverse Probability Weights] - CRAN

{{統計学}}
{{デフォルトソート:きやくかくりつおもみつけ}}
[[Category:疫学]]
[[Category:調査]]
[[Category:統計学]]
[[Category:因果推論]]
[[Category:数学に関する記事]]