統計シェアするTwitterGoogle+Pocket媒介分析シリーズ、第二段です。前回は、よく使われる媒介分析の手法の問題点についてまとめました。Analysis)を紹介するイントロとして、そもそも「媒介効果」なるものをどうやって定義するのかについてまとめます。話の性質上、少しだけ数式が登場しますが、できるだけわかりやすく書いてみようと思います。Model)の考え方および数式での表現の仕方(notation)を知っている必要があります。詳しくは前回の記事をご参照ください。前回までと同じ表記方法を使います。「介入AがアウトカムYに与える影響」に関心があるとしましょう。このような因果効果は反事実モデルの下で、と書くことができます。右下の小さい文字がポイントで、Ya=1は「もしもa=1だったときにとるであろうYの値」という意味です。Ya=0は片方しかわかりません。前回と同じ例、ジムへの入会(A=1で会員、A=0で非会員)と体脂肪率(Y%)の関係について考えます。合計200人の人がいて、A=1とA=0の人数が半々であるとすると、現実のデータは次のように見えると思います。*1個人ひとりひとりを見ると、事実と反するシナリオ下でのアウトカム(A=1のひとにとってのYa=0、A=0の人にとってのYa=1)は絶対に知ることができません。平均同士なら計算できます。ID=1~100の人たちのYa=1の平均がE[Ya=1|A=1]=E[Y|A=1]で、ID101~200の人たちのYa=0の平均がE[Ya=0|A=0]=E[Y|A=0]となります。Exchangeabilityを整えてあげれば、平均因果効果の推定完了です。繰り返しですが、忘れた人は前回記事を参照ください。媒介効果の定義ここからついに媒介効果の定義に入ります。Mediationと呼ばれることもあります。大きくわけて、二通りの定義の仕方があります。どちらかが正しいというのではなく、媒介分析の目的によって使い分けるべきものです。Assumption)も異なるので状況に応じて適切な方を選ぶといいと思います。反事実モデルに基づく媒介効果の定義は以下の論文で最初に提唱されました。疫学の世界では有名なR&Gコンビです。Scienceで有名なImai先生の以下の論文でも分かりやすくまとめてあります。765-789.AがYに影響するとき、そのような効果を引き起こす中間要因を媒介因子Mと呼ぶことにします。ジム入会(A)は運動量の増加(M)を通して体脂肪率(Y)に影響するかもしれません。DAGで書くと次のような感じです。交絡がない状況を考えています。AがMを通してYに与える影響(A→M→Y)を間接効果、AがMを介さずYに与える影響(A→Y)を直接効果と呼びます。Ya=0と表記しました。媒介効果を考えるために少し拡張して、Yamという表記を導入します。右下に小さい文字が入っているので反事実モデルですが、今度はA=aだけでなくM=mも書いてあります。Intervention下のアウトカムとして理解することができます。説明を簡単にするため、以下ではMが二値(一日30分以上の運動あり=1、30分未満または運動なし=0)であると仮定しましょう。*2Y11は「もしA=1かつM=1であったときに観測されていたであろうYの値」Y10は「もしA=1かつM=0であったときに観測されていたであろうYの値」といった具合です。Y10は「ジムに入会させておきながら、一日30分以上の運動を禁止する」という不思議な介入のシナリオです。(CDE)と呼ばれるものです。次のように定義されます。mに同じ値が入るというのがポイントです。CDE(0)はY10-Y00となり、「M=0に固定したとき(30分以上の運動の運動を禁止)のA=1とA=0(ジム入会の有無)の比較」と解釈することができます。mの値が一定に固定されているので、Y1mとY0mの差にMは全く無関係で、純粋にAの違いから生じるものだと理解することができます。運動をしていないのにジムに入会した人とそうでない人で体脂肪率に差が出たとしたら、それはジムに入会することが運動以外のメカニズムを通して体脂肪率に影響していることを示唆します。CDE(m)は「M=mに固定したときのA→Yで示された直接効果」となります。Interaction)が存在するシチュエーションを考えました。仮にそのような交互作用が存在した場合、CDE(m)はmの値によって異なる可能性があります。交互作用については以下の記事でまとめてあります。定義に入る前に、上で新たに定義したYamに加えてもうひとつ反事実モデルに基づく新たな表記方法、Maを導入します。Maとは「もしA=aだったときにとっていたであろうMの値」という意味です。Counterfactualが入り組んできたので、下の表を見ながら考えてみましょう。Consistencyを仮定すると、先ほどと同様に実際にA=1だった人たちのYaは現実のYと一致します。同じ考え方で、実際にA=1だった人たちのMa=1は現実のMと一致します。表中の上半分の人たちを見ると、実際に同じ値が入っていることがわかります。Ya=0およびMa=0には?が入っています。?はその値が観測できないことを意味します。現実にA=1だったのだから、もしA=0だったときに何が起きていたかは個人レベルではわからないのです。通常の統計的因果推論ではExchangeabilityという条件を整えて、集団レベルでE[Ya=1]-E[Ya=0]の計算を目指していくという話を前回しました。因果媒介分析では集団から集めたデータを駆使して、YaだけでなくMaも考えていく必要があるので、通常よりもきつい仮定を置く必要があります。Analysis)と必ずセットで行われるべきなのです。一般的な感度分析については次の記事を参照ください。Analysisについては次回説明します。NDEとNIEの定義さて、さきほど導入したMaを使ってNDEとNIEは以下のように定義されます。言葉でいうと、NDEは「Mをその人がA=0であったときにとっていたであろう値(M0)に固定して、A=1とA=0のときのアウトカムYを比較したもの」NIEは「A=1に固定して、MがA=1でとっていたであろう値(M1)のときとA=0でとっていたであろう値(M0)のときのアウトカムを比較したもの」と解釈できます。NDEとNIEはAがYに与える効果全体を分解したものととらえることができます。Decompositionと呼び、NDE/NIEを用いた媒介分析の特徴です。先ほど説明したCDE(m)=Y1m-Y0mとはどう違うのでしょうか。CDE(m)は全員に対してMをm同じ値に固定している一方、NDEとNIEの定義にあるM1やM0の値は一人ひとりで異なる、というのがポイントです。入れ子構造になっていて分かりにくいので、例を使って説明します。次の三人を考えてみましょう。ここでは仮想的に、(タイムマシンが存在するなど)通常観測されない反事実のアウトカムおよび曝露因子を観測可能な世界にいるとします。ID=1の人に注目しましょう。この人にとっての各値は以下のように解釈できます。M1=1であるこの人にとってのNDE/NIEはという風に計算できます。ちなみにこの人のCDEはとなります。3の人ではとなります。M1やM0の値が人によって違う点に注意です。三人の平均を考えることで、平均因果効果の媒介分析を行うことができます。Assumptionのもとで集団のデータからCDEやNDE/NIEを計算することができます。引き続き上記の例を使いましょう。EffectもNIEも0なのですがCDE(0)=1となっています。効果がないのか、あるのか困惑してしまいそうな結果です。InteractionがID=2の人に生じているからです。Effectがすべて媒介効果で説明されているように見えるのに、CDE(0)=1で直接効果があるような結果となっています。これも交互作用によるものです。最後に三人の平均をとった結果をみます。Effectとなっており、これは全て直接効果で間接効果がまったくないことを指しているようにも見えます。このようなCDEとNDE/NIEの違いはとても紛らわしいのですが、そもそも媒介分析の目的がCDEをつかうときとNDE/NIEをつかうときで全く異なることを理解すると矛盾していないことがわかります。CDEで答えることができる問いとは「Aの効果のうち、どの程度がM=mに固定する介入をしたあとに残るのか?」というものです。それに対して、NDE・NIEが答えている問いは「Aの効果のうち、どの程度がMによって説明されるか?」になります。前者は媒介因子Mに対する積極的な「介入」を想定している問いなのに対して、後者は現象としてMによるAの効果の媒介を問うています。Decompositionに関する問いということになります。したがって、なんらかの媒介因子Mに対する介入の効果を推定するなど、ポリシーメイキングに貢献するための媒介分析が目的ならばCDEのほうが相性がいいですし、逆に現象としての媒介の理解、メカニズム・Etiologyの解明に興味がある場合はNDE/NIEのフレームで媒介分析を行った方が相性がいいと思います。EffectからCDEを引いたものでもって媒介効果の推定値とするのは誤りです。Effect-CDEで媒介効果を定義してしまうと、0でない値が得られる場合があるからです(ID=2の人のパターン)。今回はA・Mがどちらも二値の場合を考えていましたが、当然それ以外のパターンにも応用できます。Regression以外をつかうときでもスケールを変えるだけでCDE・NDE・NIEは定義可能です。詳細は下記の参考資料を参照ください。Analysisを発展させたVanderWeeleが媒介分析についてまとめた一冊。反事実モデルに基づいた媒介分析に関する唯一の教科書だと思います。SASによる解析コード付きです。Kindle版この商品を含むブログを見る同じくVanderWeeleによる、比較的わかりやすいまとめ論文。https://www.annualreviews.org/doi/pdf/10.1146/annurev-publhealth-032315-021402Imaiらによる、素晴らしいまとめ論文。https://imai.princeton.edu/research/files/mediationP.pdfBaron-KennyのKennyのブログ。Analysisについて触れています。http://davidakenny.net/cm/mediate.htmKennyのブログに対する、Pearlのコメント(ブログ)。PearlはDAGの生みの親です。Mediation*1:前回紹介したConsistencyを仮定しています。*2:二値じゃない場合でも因果媒介分析はもちろん可能ですが数学的に若干ややこしくなるので。Analysis)入門2~反事実モデルに基づく媒介効果の定義~関連記事2018-12-28回帰分析を使った因果推論の仮定:パラメトリックモデルを使うということお久しぶりです。Analysis)入門1...お久しぶりです。anal...さて、今回からイデオロギー色の少ないブログ活動に戻ります。(3)スターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。