逆確率重み付け （ぎゃくかくりつおもみづけ、Inverse probability weighting、IPW）は、サンプリング母集団（データ収集に用いた母集団）とは異なる、擬似的な母集団に標準化された統計量を算出する統計手法。サンプリング母集団とターゲット母集団（推論の対象となる母集団）が異なることは多い。コストや時間、倫理的な問題などのため、ターゲット母集団から直接サンプリングすることが難しい場合もある。層化抽出法 などの代替デザイン戦略は一つの解決策である。重み付けを正しく適用することで、推定量のバイアスを低減することができる。
最初期の重み付け推定量の1つに、平均値のHorvitz–Thompson推定量 がある[1] 尤度 、重み付き推定方程式、重み付き確率密度 であり、これらから多くの統計量が派生する。これらの応用により、限界構造モデル 、標準化死亡比 、粗データや集合データに対するEMアルゴリズム など、他の統計および推定量の理論が体系化された。
逆確率重み付けは、欠測データのある被験者を一次分析に含めることができない場合に、欠測データを考慮するためにも用いられる。サンプリング確率（その因子が別の測定で測定される確率）の推定値があれば、逆確率重み付けを用いて、欠測データ が多いために過小評価されている被験者の重みを増加させることができる。
逆確率重み付け推定量（IPWE）
編集
対象実験は実施できないがモデル化可能な観測データはある場合、逆確率重み付け推定量（Inverse Probability Weighted Estimator, IPWE）を用いることで因果関係を示すことができる。治療の無作為割付は仮定しておらず、母集団のすべての被験者に対して特定の治療を割り当てた場合の反実仮想アウトカム（潜在的アウトカム）を推定することが目標となる。
観測データ
{
(
X
i
,
A
i
,
Y
i
)
}
i
=
1
n
{\displaystyle \{\left(X_{i},A_{i},Y_{i}\right)\}_{i=1}^{n}}
独立同分布 で未知の分布
P
{\displaystyle P}
X
∈
R
p
{\displaystyle X\in \mathbb {R} ^{p}}
A
∈
{
0
,
1
}
{\displaystyle A\in \{0,1\}}
Y
∈
R
{\displaystyle Y\in \mathbb {R} }
治療の無作為割付は想定していない。 まずは、潜在アウトカム
Y
∗
(
a
)
{\displaystyle Y^{*}(a)}
a
{\displaystyle a}
a
{\displaystyle a}
μ
a
=
E
(
Y
∗
(
a
)
)
{\displaystyle \mu _{a}=\mathbb {E} \left(Y^{*}(a)\right)}
{
(
X
i
,
A
i
,
Y
i
)
}
i
=
1
n
{\displaystyle \{\left(X_{i},A_{i},Y_{i}\right)\}_{i=1}^{n}}
μ
a
{\displaystyle \mu _{a}}
推定量の式
編集
μ
^
a
,
n
I
P
W
E
=
1
n
∑
i
=
1
n
Y
i
1
A
i
=
a
p
^
n
(
A
i
∣
X
i
)
{\displaystyle {\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})}}}
IPWEの構築
編集
μ
a
=
E
(
Y
1
A
=
a
p
(
a
∣
x
)
)
{\displaystyle \mu _{a}=\mathbb {E} \left({\frac {Y\,\mathbf {1} _{A=a}}{p(a\mid x)}}\right)}
p
(
a
∣
x
)
=
P
(
A
=
a
,
X
=
x
)
P
(
X
=
x
)
{\displaystyle p(a\mid x)={\frac {P(A=a,X=x)}{P(X=x)}}}
任意のプロペンシティモデル（多くはロジスティック回帰モデル）を用いて
p
^
n
(
a
∣
x
)
{\displaystyle {\hat {p}}_{n}(a\mid x)}
p
(
a
∣
x
)
{\displaystyle p(a\mid x)}
μ
^
a
,
n
I
P
W
E
=
Σ
i
=
1
n
Y
i
1
A
i
=
a
n
p
^
n
(
A
i
∣
X
i
)
{\displaystyle {\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})}}}
各治療群の平均値を算出した後、t検定 またはANOVA 検定を用いて群平均間の差を判定し、治療効果の統計的有意性を判定することができる。
一貫性：
Y
=
Y
∗
(
A
)
{\displaystyle Y=Y^{*}(A)}
未測定交絡因子がない：
{
Y
∗
(
0
)
,
Y
∗
(
1
)
}
⊥
A
∣
X
{\displaystyle \{Y^{*}(0),Y^{*}(1)\}\perp A\mid X}
治療の割り当ては、共変量データのみに基づいており、潜在的アウトカムとは無関係である。
正値性：すべての
a
{\displaystyle a}
x
{\displaystyle x}
P
(
A
=
a
∣
X
=
x
)
>
0
{\displaystyle P(A=a\mid X=x)>0}
制限事項
編集
推定される傾向が小さい場合、逆確率均等推定量（IPWE）は不安定になる可能性がある。
いずれかの治療割り当ての確率が小さい場合、ロジスティック回帰モデルはテール周辺で不安定になり、IPWEが安定しない。
拡張逆確率重み付け推定量（AIPWE）
編集
その他の推定量として、拡張逆確率重み付き推定量（Augumented Inverse Probability Weighted Estimator AIPWE）がある。これは、回帰ベースの推定量と逆確率重み付き推定量の両方のプロパティを組み合わせたもので、傾向モデルまたは結果モデルのいずれかを正しく指定するだけで、両方を指定する必要がないという点で、「二重に堅牢な」方法である。IPWEを拡張して、変動を減らし、推定効率を向上させる。このモデルは、逆確率均等推定量（IPWE）と同じ仮定を保持している。 [2]
推定量の式
編集
μ
^
a
,
n
A
I
P
W
E
=
1
n
∑
i
=
1
n
(
Y
i
1
A
i
=
a
p
^
n
(
A
i
∣
X
i
)
−
1
A
i
=
a
−
p
^
n
(
A
i
∣
X
i
)
p
^
n
(
A
i
∣
X
i
)
Q
^
n
(
X
i
,
a
)
)
=
1
n
∑
i
=
1
n
(
Q
^
n
(
X
i
,
a
)
)
+
1
n
∑
i
=
1
n
1
A
i
=
a
p
^
n
(
A
i
∣
X
i
)
(
Y
i
−
Q
^
n
(
X
i
,
a
)
)
{\displaystyle {\begin{aligned}{\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{aligned}}}
ただし、
1
A
i
=
a
{\displaystyle 1_{A_{i}=a}}
指示関数 である。回帰推定量
Q
^
n
(
x
,
a
)
{\displaystyle {\hat {Q}}_{n}(x,a)}
X
{\displaystyle X}
A
{\displaystyle A}
i
{\displaystyle i}
Y
{\displaystyle Y}
プロペンシティ推定値
p
^
n
(
A
i
∣
X
i
)
{\displaystyle {\hat {p}}_{n}(A_{i}\mid X_{i})}
ロジスティック回帰 を使用する。
AIPWEとして組み合わせて
μ
^
a
,
n
A
I
P
W
E
{\displaystyle {\hat {\mu }}_{a,n}^{\mathrm {AIPWE} }}
解釈と「二重ロバスト性」
編集
式を並べ替えると、根本的なアイデアが明らかになる。推定量は、モデルを用いて予測されたアウトカムの平均値（下記）に基づいている。
1
n
∑
i
=
1
n
(
Q
^
n
(
X
i
,
a
)
)
{\displaystyle {\frac {1}{n}}\sum _{i=1}^{n}{\Biggl (}{\hat {Q}}_{n}(X_{i},a){\Biggr )}}
しかし、モデルにバイアスがある場合、モデルの残差は0付近にはならない。モデル Q の平均残差の項を追加することにより、この潜在的なバイアスを修正できる。
1
n
∑
i
=
1
n
1
A
i
=
a
p
^
n
(
A
i
∣
X
i
)
(
Y
i
−
Q
^
n
(
X
i
,
a
)
)
{\displaystyle {\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)}
Yの値が欠落しているため、各残差の相対的な重要性を膨らませるために重みを付ける（これらの重みは、各被験者の観測値を見る逆の傾向、別名確率に基づく）（ [3]
「二重ロバスト性」は、推定量が不偏であるためには、
Q
^
n
(
X
i
,
a
)
{\displaystyle {\hat {Q}}_{n}(X_{i},a)}
p
^
n
(
A
i
∣
X
i
)
{\displaystyle {\hat {p}}_{n}(A_{i}\mid X_{i})}
[3] [4] [5]
二重ロバスト推定量のバイアスは2次バイアス と呼ばれ、
1
p
^
n
(
A
i
∣
X
i
)
−
1
p
n
(
A
i
∣
X
i
)
{\displaystyle {\frac {1}{{\hat {p}}_{n}(A_{i}\mid X_{i})}}-{\frac {1}{p_{n}(A_{i}\mid X_{i})}}}
Q
^
n
(
X
i
,
a
)
−
Q
n
(
X
i
,
a
)
{\displaystyle {\hat {Q}}_{n}(X_{i},a)-Q_{n}(X_{i},a)}
機械学習 推定器を使用して、二重ロバスト推定量の全体的なバイアスを下げることができる[6]
関連項目
編集
^ Horvitz, D. G.; Thompson, D. J. (1952). “A generalization of sampling without replacement from a finite universe”. Journal of the American Statistical Association 47 (260): 663–685. doi :10.1080/01621459.1952.10483446 . ^ Cao, Weihua; Tsiatis, Anastasios A.; Davidian, Marie (2009). “Improving efficiency and robustness of the doubly robust estimator for a population mean with incomplete data” . Biometrika 96 (3): 723–734. doi :10.1093/biomet/asp033 . ISSN 0006-3444 . PMC 2798744 . PMID 20161511 . https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2798744/ . ^ a b Kang, Joseph DY, and Joseph L. Schafer.
^ Kim, Jae Kwang, and David Haziza.
^ Seaman, Shaun R., and Stijn Vansteelandt.
^ Hernán, Miguel A., and James M. Robins.
外部リンク
編集
