「学術的に理解する」ってどういうこと? ：統計的因果推論のフレームワークから見た 「量的×質的」に関する眺望の (いささかとっちらかった)スケッチを共有する

統計的因果推論のフレームワークとはいかなるものかを理解することで、統計的因果推論がいかに質的な知見を（しばしば暗黙に）その基盤としているか、統計的因果推論のスコープがしばしば いかに”狭い&薄い”のかを議論しました。前半はDAG-潜在結果モデル-構造的因果モデルの統一的な解説です。後半はSUTVAと質的研究の話、固有性-法則性と斎藤清二先生のナラエビの話、筒井淳也(2019)と"キャンセルアウト"の是非の話、です。

Transcript

1. 先端課題研究会＠ 2021/7/21 @国⽴環境研究所 社会システム領域 林岳彦 「学術的に理解する」ってどういうこと？ 統計的因果推論のフレームワークから⾒た 「量的×質的」に関する眺望の （いささかとっちらかった）スケッチを共有する
2. 統計的因果推論のフレームワークとは いかなるものかを理解することで 統計的因果推論がいかに質的な知⾒を （しばしば暗黙に）その基盤としているか 統計的因果推論のスコープがしばしば いかに”狭い＆薄い”のか を議論したい 本⽇のねらい 質的×量的に関して内容的には全く新規性のある話ではないが、「統計的因果推論の枠組みからは こう⾒える」という幾つかの眺望のスケッチを共有でき、それが何らかのヒントとなれば幸いです
3. ＊私は質的研究には詳しいわけではありません にわか勉強としてひとまず以下の本は読みました （その結果、本⽇の内容は何らかの形で既に⾔及されていると思った次第）
4. ⾃⼰紹介｜どんな研究をしているのか もともとの出⾝は⽣態学（性淘汰・種分化理論） 『性選択理論を整理する』 （性淘汰理論を統⼀的に解説） 種分化過程の個体ベース シミュレーション研究 Figure 2. The distribution of locus effects in the genetic distance model with two loci. (A) Differentiation in the female locus only (th “Buridan’s Ass” regime). The population is from the simulation shown in Figure 1B at generation 7500. (B) Differentiation in both femal and male loci. The population is from the simulation shown in Figure 1B at generation 25,000. (C) Sympatric speciation with strong sexua conflict (Popt = 0.2). Upper white bars: female locus. Lower black bars: male locus. Parameters are the same as in Figure 1B unless specified case) of a haploid network when genetic diversification in both sexes occurred. It shows a population consisting of small, loosely connected clusters, rather than large and distinct clusters. The web structure among clusters shows that recombinant genotypes were common in the population. With the same parameter configuration, genetic differentia- tion is more likely in the eight-locus model (Table 2a) than in the two-locus model (Table 1a). This partly reflects the difference in the strength of selection pressure on each locus. When the number of loci is small, the selection pressure on each locus is strong and tends to suppress the genetic variance required to initiate genetic differentiation. Moreover, differentiation rarely occurs once co- evolutionary chase begins. In contrast, genetic differentiation is still possible even when coevolutionary chase occurs in one of the Figure 3. Examples of haplotype networks. The ovals represent groups of individuals with the same haplotype. The length of a branc (more precisely, the number of nodes between clusters) represents the distance between haplotypes. These haplotype networks wer computed based on 50 individuals randomly sampled at generation 25,000. Parameters: Popt = 0.2, ␣ = 0.01, sc = 1.02, ␮ = 0.5 × 10−5 N = 10,000, unless specified. (A) Sympatric speciation in the two-locus genetic distance model with codominance. The data are from th simulation run shown in Figure 2C, ␣ = 0.05. (B) Genetic diversification without speciation in the eight-locus genetic distance mode with directional dominance. The data are from the simulation run shown in Figure 4A. (C) Sympatric speciation in the eight-locus geneti distance model with codominance, sc = 4 × 1.02. (D) Genetic diversification without speciation in the 32-locus genetic distance mode pairs of loci, provided that the number of loci is moderate (e.g L = 8). The codominance case.—Genetic differentiation in th loci of both sexes was often observed (Fig. 4B and Table 2b), al though the frequency of differentiation in male loci was smalle than in the directional dominance case. When selection in female was not sufficiently strong (i.e., sc = 1.02), no sympatric speci ation was observed (because of recombination among divergin loci). Strong female preference (␣ = 0.05) enhanced genetic dif ferentiations in the loci of both sexes but did not cause sympatri speciation (Table 2c). Decreased mutation rate (␮ = 10−5) sup pressed genetic differentiation, especially in male loci (Table 2d) No diversification was observed when population size was smal (N = 1000). 配偶者選択の基準が多様化することにより、 集団内に遺伝的⽂化が⽣じる過程を研究 ORIGINAL ARTICLE doi:10.1111/j.1558-5646.2007.00059.x GENETIC DIFFERENTIATION BY SEXUAL CONFLICT Takehiko I. Hayashi,1,2 Michael Vose,3,4 and Sergey Gavrilets5,6 1Research Center for Chemical Risk Management, National Institute of Advanced Industrial Science and Technology, 16-1 Onogawa, Tsukuba, Ibaraki 305-8589, Japan 2E-mail: ti-hayashi@aist.go.jp 3Department of Computer Sciences, University of Tennessee, Knoxville, Tennessee 37996 4E-mail: vose@cs.utk.edu 5Department of Ecology and Evolutionary Biology, Department of Mathematics, University of Tennessee, Knoxville, Tennessee 37996 6E-mail: gavrila@tiem.utk.edu Received June 6, 2006 Accepted October 31, 2006 Sexual conflict has been suggested as a general cause of genetic diversification in reproductive characters, and as a possible cause of speciation. We use individual-based simulations to study the dynamics of sexual conflict in an isolated diploid population with no spatial structure. To explore the effects of genetic details, we consider two different types of interlocus interaction between female and male traits, and three different types of intra-locus interaction. In the simulations, sexual conflict resulted in at least the following five regimes: (1) continuous coevolutionary chase, (2) evolution toward an equilibrium, (3) cyclic coevolution, (4) extensive genetic differentiation in female traits/genes only, and (5) extensive genetic differentiation in both male and female traits/genes. Genetic differentiation was hardly observed when the traits involved in reproduction were determined additively and interacted in a trait-by-trait way. When the traits interacted in a component-by-component way, genetic differentiation was frequently observed under relatively broad conditions. The likelihood of genetic differentiation largely depended on the number of loci and the type of within-locus dominance. With multiple loci per trait, genetic differentiation was often observed but sympatric speciation was typically hindered by recombination. Sympatric speciation was possible but only under restrictive conditions. Our simulations also highlight the importance of stochastic effects in the dynamics of sexual conflict. Hayashi et al. (2006) 研究の原点はフィールド観察にある、という感覚のある分野
5. ⾃⼰紹介｜どんな研究をしているのか ⽣態リスク評価、確率的リスク分析、因果推論 Yij θAlgae σj Means SD θFish Sensitivity differences among taxonomic groups Normal distributions Log(NOEC) Parameters were estimated by MCMC simulations θInvertebrate 階層ベイズモデルとモンテカルロシミュレーションを ⽤いた化学物質の⽣態リスクの定量化とリスク⽐較 Hayashi & Kashiwagi (2009) Hayashi & Kashiwagi (2010) Monte Carlo Analysis EPAF = F µ ECD - µ SSD s ECD 2 + s SSD 2 æ è ç ç ö ø ÷ ÷ µ ECD µ SSD s ECD s SSD Calculation of predictive distribution of EPAF Posterior distributions of ECD parameters Posterior distributions of SSD parameters Results: Quantitative Risk Comparison Median and 90% range of EPAF log10(EPAF) Large Risk→ ←Small Risk Chemicals Ammonia Copper Nickel Zinc Hayashi and Kashiwagi (2011) リスク／影響の確率論的解釈（個別-法則）
6. ⾃⼰紹介｜どんな研究をしているのか 『はじめてのバックドア基準』 （Pearlの関数的因果モデルの解説） ⽣態リスク評価、確率的リスク分析、因果推論 ネオニコチノイド系農薬の ⾚トンボへの影響の因果推論 !",$ = &",$ − &"(),$ &"(),$ = +) ,-./)," − ,-./),"() + ⋯ + +2 ,-./2," − ,-./2,"() ++" 3-45" + 6) 7)," + ⋯ + 68 78,"
   
    
                 背景知識からの因果モデルの構築 バックドア基準に基づく 統計モデルの構築 因果効果の推定 Takeshita et al. (2020)など 介⼊効果とドメイン知識
7. ⾃⼰紹介｜どんな研究をしているのか 社会対話、リスクコミュニケーション、EBPM @taiwa_kankyo の twitterアカウント 社会対話・協働推進オフィス＠国環研のメンバーとして コミュニケーション活動に参画 “エビデンス”とナラティブ
8. ⾃⼰紹介｜どんな研究をしているのか 社会対話、リスクコミュニケーション、EBPM “エビデンス”を政策利⽤する上で考慮されるべき論点を 科学哲学者との共同研究により5x3表の形で整理 Environmental Science and Policy 116 (2021) 86–95 Fig. 1. Cross table of managerial points for EBPM in environmental policymaking. Note. Each perspective comprises of two levels: the issues involved in each phase, and a checklist for managing each stage in EBPM. H. Kano and T.I. Hayashi Environmental Science and Policy 116 (2021) 86–95 A framework for implementing evidence in policymaking: Perspectives and phases of evidence evaluation in the science-policy interaction Hiroyuki Kano a,*, Takehiko I. Hayashi b a Graduate School of Human Sciences, Osaka University, 1-2 Yamadaoka, Suita, Osaka, 565-0871, Japan b Center for Health and Environmental Risk Research, National Institute for Environmental Studies, 16-2 Onogawa, Tsukuba, Ibaraki, 305-8506, Japan A R T I C L E I N F O Keywords: Evidence-based policy Science-policy interface Reference guide Institutionalization Mercury A B S T R A C T The use of scientific knowledge in policymaking has been a subject of debate in the environmental sector. An essential task for the effective use of evidence in policymaking is for scientists and policymakers to share a common understanding of how evidence should be produced and used. The purpose of this study is to establish a reference framework that enables scientists and policymakers to align their sights to deal with evidence in policymaking. To develop the framework, we introduced five perspectives that cover the domains of science, policy, and the science-policy interface as: (1) methodological rigorousness; (2) consistency; (3) proximity; (4) social appropriateness; and (5) legitimacy. We then examined how the issues from these perspectives will transit through three phases of interaction between scientific investigation and political institutionalization via: (1) a pre-institutionalization phase, in which the academic and social framing of an issue was unclear; (2) a mid- institutionalization phase, in which the academic and social framing of an issue was established and institu- tionalization was advanced; and (3) a post-institutionalization phase, in which issues were recursively defined within the evidence evaluation system itself. The framework encourages scientists and policymakers to shift their perspectives to each phase of institutionalization. A case study on mercury pollution shows how the framework serves as a checklist for the comprehensive evaluation of evidence, which provides specific guidance that appropriately promotes evidence-based policymaking and its implementation. 1. Introduction The use of scientific knowledge in policymaking has been a subject of debate in the environmental sector. Although evidence-based policy- making has been advocated over the last twenty years, the idea has not transparent process for incorporating research evidence into policy de- cisions (Dicks et al., 2014); a systematic climate policy evaluation (Haug et al., 2010); a participatory and inclusive structures for policymaking (Lalor and Hickey, 2014); social interactions, socio-political environ- ments, and power relations in the knowledge production and its use Contents lists available at ScienceDirect Environmental Science and Policy journal homepage: www.elsevier.com/locate/envsci Kano and Hayashi (2021) “エビデンス”と公共政策
9. ⾃⼰紹介｜どんな研究をしているのか 学際的に複数分野を漂流してきた 統計のエンドユーザーとして （しばしば明⽰的／⾮明⽰的に存在する） ドメインの知識／ニーズ／規範と 統計解析を適切に接続するための仕事をしてきた 「学術的に理解する」とは どういうことか？ （ドメイン知識と統計解析の狭間で） Issue-drivenの分析者はこういうことを考えることがままあると思う
10. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（IV部構成） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在結果モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
11. 本⽇のテーマの プレビュー
12. 本⽇のテーマのプレビュー (1)可能世界と反事実条件、および”世界間距離” (2)概念の同⼀性／異質性／妥当性 (3)法則性と固有性 (4)われわれは何が知りたいのか？ ちょっとこの件だけ少し前置きとして説明させてください
13. 前置き：反事実条件⽂と因果概念の定義 • 可能世界論からの「必然性」「可能性」の考え⽅ 「可能世界」の枠組みで考えれば⾊々と 捗るよ！ S.クリプキ (1940-) 「Xが必然である」とは？ 全ての可能世界においてXが成り⽴つ 「Xが可能である」とは？ 少なくとも１つの可能世界においてXが成り⽴つ ・世界について考えうる異なる「あり⽅」ごとに異なる 「可能世界」がある ・その中で我々が実際に暮らしているのこの世界が 「現実世界」 可能世界を考えることで、「必然性」や「可能性」という概念を 論理的に定式化できる image
14. 前置き：反事実条件⽂と因果概念の定義 • 可能世界論からの「因果」の考え⽅ D.ルイス (1941-) 「XがYの原因である」 「同じ状況（到達可能な近傍の可能世界） において、もしもXが起こらなければ、Yは 起こらなかっただろう」 反事実依存性に基づく因果解釈 「covid-19の流⾏がオンラインでの授業の実施の原因である」 「 covid-19の流⾏が起こらなければ、 オンラインでの授業は実施されなかっただろう」 image
15. 前置き：反事実条件⽂と因果概念の定義 • 反事実依存性と「因果推論の根本問題」 「同じ状況（到達可能な近傍の可能世界）において、もしも Xが起こらなければ、Yは起こらなかっただろう」 因果効果を次のように定義できそう： X→Yへの因果効果 ＝出来事Xが起きた世界におけるY - 出来事Xが起きなかった世界におけるY しかし、「起きた世界」と「起きなかった世界」の両⽅ を同時に観測することは原理的に不可能 “因果推論の根本問題” どちらかは必ず 反事実的な世界となる 例. コロナ禍
16. 前置き：反事実条件⽂と因果概念の定義 • 反事実依存性と「因果推論の根本問題」 「同じ状況（到達可能な近傍の可能世界）において、もしも Xが起こらなければ、Yは起こらなかっただろう」 因果効果を次のように定義できそう： X→Yへの因果効果 ＝出来事Xが起きた世界におけるY - 出来事Xが起きなかった世界におけるY しかし、「起きた世界」と「起きなかった世界」の両⽅ を同時に観測することは原理的に不可能 “因果推論の根本問題” どちらかは必ず 反事実的な世界となる 例. コロナ禍 “反事実的な世界”を 統計的に構成することにより この「根本問題」を解決するのが 「統計的因果推論」である （起こらなかった⽅の可能世界）
17. 反事実条件⽂による因果定義の特徴 “⽔がからからになり、林の居室の花は枯れた” もし林が⽔をあげていたら、居室の花は枯れ なかっただろう もしドナルド・トランプが⽔をあげていたら、 居室の花は枯れなかっただろう 林が⽔をあげなかったことが原因 トランプが⽔をあげなかったことが原因？ (A) (B) 反事実条件⽂の内部には可能世界間の”距離”に ついての情報が皆無 ←この⽂も真！ ←この⽂は真 c.f., p. 147 in Sloman (2005) Causal models.
18. 本⽇のテーマのプレビュー (1)可能世界と反事実条件、および”世界間距離” (2)概念の同⼀性／異質性／妥当性 (3)法則性と固有性 (4)われわれは何が知りたいのか？
19. 本⽇のテーマのプレビュー v IL Morgan and Winship (2014) 『Counterfactuals and Causal Inference: Methods and Principles for Social Research』 (2nd ed.) 本⽇の”こころの元ネタ” 今回具体的にこの本の内容の話をする わけではないが、話のスコープの 根っこにはこの本のモチーフがあります
20. 本⽇のテーマのプレビュー 因果推論本おすすめフローチャートの例 https://www.bradyneal.com/which-causal-inference-book より引⽤
21. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
22. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり 筑波⼤学での分担講義 「環境社会学：環境データ分析」の ⼀部をベースにしたもの 本⽇の構成（パートI-IV） 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る なので質的×量的への⾔及がところどころ挟まりますが、 基本的には授業っぽい内容です
23. 前置き｜統計的推測の主要な3⽬的 ⽬的の 分類 やりたいことの内容 事例のイメージ （⼤学の留年の例で） 理念的な理想をなす 解析法のイメージ しばしば 現実的な 解析法 予測 現在得られている情報 から、将来／未知の値 ／状況を予測したい （要因への介⼊は伴わない場合） 多数の変数の情報から、留年 リスクの⾼い学⽣を⾒つけた い（→その後の⽀援につなげ るなど） 機械学習 回帰分析！ 説明・記述 要素間の関連性を説明 したい なぜ学⽣が留年するかを説明 したい（要因を探索し理解／ 対策につなげたい） フルスタックの 構造⽅程式 （含む因⼦分析等） 回帰分析！ 因果／介⼊ 効果の推定 ある要因Xに介⼊したと きの、要因Yに及ぼす効 果を知りたい 特定のプログラムの導⼊が留 年を防⽌する効果を知りたい ランダム化⽐較試験 回帰分析！ ⼀⾒同じ「回帰分析」でも、⽬的に応じて係数の 解釈のあり⽅も、重要となる仮定群も異なる （異論は認める） ＊ ＊互いに排反な分類という意味ではないので注意 ↑ パートI,IIでは「因果効果の推定」が主テーマ ▪ （余談）
24. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
25. 前置き パートIのテーマ：「構造」と「差分」の視点の統合 「対象を理解すること」にとって 統計的因果推論とは何か？ ⾮巡回有向グラフ(DAG) 潜在結果モデル 構造的因果モデル 𝒀 = 𝑿𝒀𝟏 + (𝟏 − 𝑿)𝒀𝟎 X Y Z1 Z2 Z3 Average Causal Effect = 𝑬[𝒀𝟏] − 𝑬[𝒀𝟎] 𝑋: = 𝑔: 𝑝𝑎 𝑋: , 𝜀: 𝑖 = 1, … , 𝑝 「差異からみた因果」 「構造からみた因果」
26. バックドア基準は、背景にある因果構造を基に「バ イアスのない因果効果推定」のためにはどの変数を 考慮するべきかの理論的条件を与える (Preview) パートIの話の全体まとめ ▪ 潜在結果モデルは「潜在結果の差分」として因果効 果を定式化し、因果効果推定の解析プロトコルへの 明瞭な⾒通しを与える ▪ 構造的因果モデルは⾮常に広いクラスの因果効果の 推定アプローチを包含しうるものであり、issueの総 合的理解へ向けた統計的／理論的／理念的な枠組み の基盤となりうる ▪
27. 1. 構造から⾒た統計的因果推論 ：DAGとバックドア基準の⼊⾨ 2. 差分から⾒た統計的因果推論 ：潜在結果モデル⼊⾨ 3. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在結果モデルのSCMによる表現 (3) 潜在結果モデルとDAGとの繋がり 4. 差分と構造の統合 --- Issueの総合的理解へ パートIの⽬次 0. 序：因果と相関と「対象の多様なありよう」
28. 1. 構造から⾒た統計的因果推論 ：DAGとバックドア基準の⼊⾨ 2. 差分から⾒た統計的因果推論 ：潜在結果モデル⼊⾨ 3. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在結果モデルのSCMによる表現 (3) 潜在結果モデルとDAGとの繋がり 4. 差分と構造の統合 --- Issueの総合的理解へ パートIの⽬次 0. 序：因果と相関と「対象の多様なありよう」
29. 因果推論のはじまりの場所｜ 初⼼に戻り「対象のありよう」を丁寧に考える 10個のリンゴがあります 統計の教科書が想定してる「ありよう」は Exchangeable!
30. Exchangeable? 因果推論のはじまりの場所｜ 初⼼に戻り「対象のありよう」を丁寧に考える 10個のリンゴがあります 現実はしばしば「ふぞろいのリンゴ」である
31. 因果推論のはじまりの場所｜ 因果推論でまず重要なのは「対象集団のありよ う」について丁寧に省察することである ある意味、統計的因果推論とは集団内の異質 性／⾮均⼀性を巡る体系である（後述） 「ありよう」の質的な省察は推論の⼟台をなす Exchangeable?
32. 相関と因果がズレるとき｜ 仮想例として「肥料X→リンゴの糖度Y」の因果 効果を考えてみる（まずは均⼀品種バージョン） 肥料Xを与えると糖度Yは単純に+2される もともとのリンゴの糖度Yの平均は16、 分散は1.0 「肥料X=あり」で育てたリンゴは50個 「肥料X=なし」で育てたリンゴは50個 肥料X=あり/なしの糖度Yをプロットすると つまり「真の因果効果」は +2.0
33. 3 8 12 16 20 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 相関と因果がズレるとき｜ なし あり 肥料X 糖 度 Y 「肥料X→リンゴの糖度Y」の散布図（1品種ver） +2.1 各処理グループ 平均の差(+2.1) 「真の因果効果(+2)」 ≒ 散布図上での差を 因果効果として そのまま解釈可能 ＝バイアスなし ここでのバイアスの定義＝ 「真の因果効果」と「観測された処理 グループ平均の差」の系統的なズレ
34. 相関と因果がズレるとき｜ 「ぺこ」と「すまいる」の２つのリンゴ 品種がサンプル内に混在している 肥料X=あり/なしの糖度Yをプロットすると 元々の「ぺこ」の糖度Yの平均は16 元々の「すまいる」の糖度Yの平均は12 その他の設定・仮定は先程の例と同一 ひきつづき「肥料X→リンゴの糖度Y」の因果効 果を考えてみる（２品種への拡張バージョン）
35. 3 8 −1 0 1 2 3 3 8 12 16 20 −1 0 1 2 3 相関と因果がズレるとき｜ なし あり 肥料X 糖 度 Y 「肥料X→リンゴの糖度Y」の散布図（２品種ver） +4.4 各処理グループ 平均の差(+4.4) 「真の因果効果(+2)」 ≠ （他の設定は同⼀でも） ２種類の品種が 混在するだけで バイアスが⽣じる ここでのバイアスの定義＝ 「真の因果効果」と「観測された処理 グループ平均の差」の系統的なズレ
36. 相関と因果がズレるとき｜ 品種が混在するときは必ずバイアスが⽣じる？ 8 12 −1 0 1 2 3 8 12 −1 0 1 2 3 8 12 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 20 20 20 +4.4 +2.0 -0.8 なし あり 肥料X 糖 度 Y なし あり 肥料X なし あり 肥料X 40/10 10/40 25/25 25/25 40/10 10/40 「品種の⽐率」が処理グループ間で同じとき にはバイアスが⽣じない ぺこ/すまいる比（総数ではぺこ50個/すまいる50個）
37. 8 12 −1 0 1 2 3 8 12 −1 0 1 2 3 8 12 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 8 12 16 20 −1 0 1 2 3 相関と因果がズレるとき｜ 品種が混在するときは必ずバイアスが⽣じる？ +3.3 +2.2 +0.8 なし あり 肥料X 糖 度 Y なし あり 肥料X なし あり 肥料X ぺこ/すまいる比 8/42 2/48 5/45 5/45 8/42 2/48 「品種の⽐率」が処理グループ間で同じとき にはバイアスが⽣じない （総数ではぺこ10個/すまいる90個） ここ⼤事
38. 相関と因果がズレるとき｜ 「品種の⽐率」が処理グループ間で同じとき にはバイアスが⽣じない 「特性（共変量）のありよう」が処理グループ 間で同じときにはバイアスが⽣じない 統計的因果推論の要点：処理と共変量の独⽴性 処理グループ間で「共変量がバランシングしている」みたいな ⾔い⽅もします この例を⼀般化すると 数式だと P(共変量Z|X=あり)=P(共変量Z|X=なし) 統計的因果推論の主要な企みとは、上記の独⽴性 を実験計画や解析の⼯夫により達成すること
39. 相関と因果がズレるとき｜ 対処例：層別化による処理と共変量の「独⽴化」 鱿鱘ݸ
    
    鱞鲃鱉鲐ݸ 鱿鱘ݸ
    
    鱞鲃鱉鲐ݸ ʮංྉ9鱇鲏ʯ鱳಺༁ ʮංྉ9鱯鱜ʯ鱳಺༁ ౶ ౓ ංྉ9 鱯鱜 鱇鲏  ౶ ౓ ංྉ9 鱯鱜 鱇鲏  ౶ ౓ ංྉ9 鱯鱜 鱇鲏  『 』 層別解析 『 』 層別解析 層別化で処理間の共変量を 揃える→バイアスが消える
40. 相関と因果がズレるとき｜ ここまでの⼩まとめ サンプルの特性がまったく均⼀のときにはバ イアスが⽣じない ＊ここでのバイアスの定義＝「真の因果効果」と「観測された処理グループ平均の差」 の系統的なズレ 処理グループ間で特性（共変量）の分布が同じ ときにはバイアスが⽣じない この意味で、因果推論とは「対象集団におけ る異質性への対処」を巡る体系である 因果推論の主要な企みとは、この状況を実験 計画や解析の⼯夫により達成すること ＊また、そもそも異質性（たとえば品種）を⾒分けることができなければバイアスの 存在にも気づけない（→質的な省察が推論の⼟台となる）
41. 相関と因果がズレるとき｜ ひとくちに「特性の分布が同じ」といっても 無数の「特性」がありうるわけで… その全ての特性（全ての共変量）を 処理間で揃えなければならないの？ ここで少し考えてみよう 「揃えるべき特性」と「揃えなくてもよい特性」を 判別する理論的基準はあるのか？ ＊たとえばリンゴの「特性」を考えても、品種、産地、農法、農家、流通形態、 収穫時期、価格、などなどいくらでもありうるわけで… 細かく考えていくときりがない！
42. 導⼊｜相関と因果は違う 散布図上の相関関係からの素朴な予測や期待と 因果効果（介⼊効果）がズレることがある ▪ ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : p < 0.001 「Xが⼩さいとき、Yは⼤きい」 （=相関関係がある）とほぼ確実に⾔える
43. 導⼊｜相関と因果は違う 散布図上の相関関係からの素朴な予測や期待と 因果効果（介⼊効果）がズレることがある ▪ ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : p < 0.001 「Xが⼩さいとき、Yは⼤きい」 （=相関関係がある）とほぼ確実に⾔える 「Xを⼩さくすると、Yは⼤きくなる」 （=因果関係がある）と⾔えるか？
44. 導⼊｜相関と因果は違う もしXとYの相関が「Xと関係ない要因」で⽣じ ているならば、Xを減少させてもYは回復しない ▪ ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : ԼྲྀҬ ্ྲྀҬ
45. 導⼊｜相関と因果は違う もしYの状態がXにより⽣じているならば、Xを 減少させればYは回復する ▪ ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : ԼྲྀҬ ্ྲྀҬ
46. 導⼊｜相関と因果は違う もしYの状態がXにより部分的に⽣じているなら ば、Xを減少させればYは部分的に回復する ▪ ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : ԼྲྀҬ ্ྲྀҬ 「Xを⼩さくすると、Yは⼤きくなる」かは Yの状態が何によって⽣じているかに依存する
47. 導⼊｜ 相関と因果で必要とされる⽅法論が違う 「Xが⼩さいとき、Yは⼤きい」か？ ∈「相関関係(association)」についての問い p(Y|X=⼩さい) やるべきこと：Xが⼩さいときのYの計算 ←所与の数値のみから計算できる ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : p < 0.001
48. 導⼊｜相関と因果で必要とされる⽅法論が違う 「Xを⼩さくすると、Yは⼤きくなる」か？ ∈「因果的影響(causal effect)」についての問い p(Y|X=do(⼩さい)) やるべきこと：Xを⼩さくしたときのYの計算 ←計算に因果的情報が必要 ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : ԼྲྀҬ ্ྲྀҬ
49. 導⼊｜相関と因果で必要とされる⽅法論が違う p(Y|X=⼩さい) Xが⼩さいときのYの計算 p(Y|X=do(⼩さい)) Xを⼩さくしたときのYの計算 因果の問題 相関の問題 p(Y|X=see(⼩さい)) あるいは、いうなれば なぜズレるの？ 交絡のせい！ 内⽣性のせい！ ⼀歩込み⼊った議論をする際に判断に迷うことがある 最初の ⽬標 「バックドアパスが開いているせい」 「バックドア基準が満たされていないせい」 とい説明の仕⽅・概念を理解する association
50. 前置き｜バックドア基準は「何について」の話？ ▪ 例えば、重回帰モデルの場合 モデルに追加された説明変数の組が「X→Y」についてバックドア基準を満 たすとき、重回帰分析から得られたXの偏回帰係数をそのまま「X→Yの介 ⼊効果」のバイアスのない推定量とみなせる 例えば、「シンプソンのパラドックス」(ここでは、どの変数で層別化す るかによって推定結果が変るケースの意味で⽤いる）が⽣じている場合 興味の対象となる「処理X→結果Y」についてバックドア基準を満たす変数 で層別化して解析すれば「X→Yの介⼊効果」をバイアスなく推定できる 具体的には、バックドア基準を満たすと： ＊データが適切に測定されており、かつ適切なモデルが適⽤されているという⼤前提での話です バックドア基準は介⼊効果推定での変数選択の基準 をもたらす
51. 今⽇の話｜ バックドア基準とは？（ざっくり） “バックドア基準が満たされている” ▪ “バックドアパス”って何？ “開く／閉じる”ってどういうこと？ + (2) 処理X→結果Yの道がブロックされていない (1) 開きっぱなしのバックドアパスがない ≒ “ブロック”って何？
52. 今⽇の話｜ バックドア基準とは？（ざっくり） “バックドア基準が満たされている” ▪ “バックドアパス”って何？ “開く／閉じる”ってどういうこと？ + (2) 処理X→結果Yの道がブロックされていない (1) 開きっぱなしのバックドアパスがない ≒ “ブロック”って何？
53. 超ざっくり説明｜”バックドアパス”とは？ n 因果構造を丘にある「⼈⼯池」でイメージする 迂遠かもしれませんが、因果の「流れ」のアナロジーと して有効だと思うので少々お付き合いいただければ・・ X Y ＊降⾬の影響などは考えない 本講演を通して X: 処理 Y: 結果 Z: 共変量 で表記します ＊ Y=βX+γ+ε ⼈⼯池 ⼈⼯池 ⽔路
54. 超ざっくり説明｜”バックドアパス”とは？ n ⽤語説明：⾮巡回有向グラフ DAG⾃体は⾮常に広いクラスのモデルを表現しうる。たとえば、原理的には、時系列 的なsequentialな数値シミュレーションを書ける系はDAGの形で描ける X Y “親” “⼦” “祖先” “⼦孫” (Directed Acyclic Graph; DAG) 𝑎;<= 𝑎; 𝑎;>= 𝑌;<= 𝑌; 𝑌;>= 状態空間モデルの例
55. 超ざっくり説明｜”バックドアパス”とは？ 「バックドアパス」=『処理Xと結果Yの上 流側にある両者に影響を与える流れ』 ▪ X Y Z1 Z2 （X→Yの） バックドアパス 上流側にある「Z1」に インクをぶちまけると XにもYにも到達する (＝両者に影響を与える)
56. 超ざっくり説明｜”バックドアパス”とは？ 「バックドアパス」=『処理Xと結果Yの上 流側にある両者に影響を与える流れ』 ▪ X Y Z1 Z2 （X→Yの） バックドアパス ではない！ 上流側にある「Z1」に インクをぶちまけても Yにしか到達しない (＝両者には影響を与えない)
57. 超ざっくり説明｜”バックドアパス”とは？ 「バックドアパス」=『処理Xと結果Yの上 流側にある両者に影響を与える流れ』 ▪ X Y Z1 Z2 （X→Yの） バックドアパス ではない！ 上流側にある「Z1」に インクをぶちまけても Xにしか到達しない (＝両者には影響を与えない)
58. 超ざっくり説明｜”バックドアパス”とは？ 「バックドアパス」=『処理Xと結果Yの上 流側にある両者に影響を与える流れ』 ▪ X Y Z1 Z2 （X→Yの） バックドアパス ではない！ 上流側にある「Z1」に インクをぶちまけても Xにしか到達しない (＝両者には影響を与えない) ＊Xそのものを通してYへ繋がる 流れはノーカウント
59. 超ざっくり説明｜バックドアパスの何が問題？ 上流側の変動に伴い”シンクロ”が⽣じる ▪ X Y Z X Y Zの変動によりX-Y間に "シンクロ"が⽣じる Zが⼤ Zが⼩ （X→Yの） バックドアパス このZの変動によるシンクロが”疑似相関”を⽣み、 X→Yの介⼊効果の推定にバイアスをもたらす ⾮因果的 連関
60. 超ざっくり説明｜バックドアパスの何が問題？ 参考：バックドアパスがない場合の例 ▪ X Y Z X Y Zの変動でX-Y間に “シンクロ”は⽣じない Zが⼤ Zが⼩ バックドアパスがないとき、Zの変動はバイアス・ 擬似相関の原因とならない
61. 今⽇の話｜ バックドア基準とは？（ざっくり） “バックドア基準が満たされている” ▪ “バックドアパス”って何？ “開く／閉じる”ってどういうこと？ + (2) 処理X→結果Yの道がブロックされていない (1) 開きっぱなしのバックドアパスがない ≒ “ブロック”って何？
62. 超ざっくり説明｜バックドアパスが”閉じる”とは？ バックドアパスが”閉じる”=バックドアパス 上の変数を”固定”する ▪ X Y Z Zを“固定”する バックドアパスが”閉じ” 交絡の影響が消える シンクロを⽣む流れが “ブロック”される Zで層別化して解析 重回帰分析の共変量としてZをモデルに追加 “固定”= 条件付けconditioning の⼿法の例 ・ ・
63. 超ざっくり説明｜バックドアパスが”閉じる”とは？ 参考：Zでの層別化（X：汚染濃度, Y：種数, Z：流域） ▪ 「上流」「下流」の層別に X→Yの効果を計算して集計する Z:流域 X:汚染 濃度 Y:種数 流域Zを層別化により”固定”し解析することで交絡 の影響が消え、介⼊効果がバイアスなく推定できる ؀ڥԚછ෺࣭ͷՏ઒தೱ౓ 9 ఈ ੜ ࠛ ஬ ͷ छ ਺ : ԼྲྀҬ ্ྲྀҬ
64. 超ざっくり説明｜パスが”開きっぱなし”とは？ バックドアパスが開きっぱなし＝変数が固定 されずパスの流れが”ブロック”されていない ▪ X Y Z X Y Zの変動によりX-Y間に "シンクロ"が⽣じる Zが⼤ Zが⼩ （X→Yの） バックドアパス ⾮因果的 連関 介⼊効果の推定のバイアスが調整されぬまま残る
65. 超ざっくり追記｜固定により”開く”こともある 合流点を固定すると、“親”の間に双⽅向 パスが開く ▪ バックドアパスが “開いた”状態 X Y Z （X→Yの） バックドアパス ⼦ 親 親 バックドアパスの ない状態 X Y Z Zを“固定
66. 超ざっくり追記｜固定により”開く”こともある 合流点を固定すると、“親”の間に双⽅向 パスが開く ▪ バックドアパスが “開いた”状態 X Y Z （X→Yの） バックドアパス ⼦ 親 親 バックドアパスの ない状態 X Y Z Zを“固定 “⼦はかすがい”
67. 超ざっくり追記｜固定により”開く”こともある 参考：合流点となる「合否」を”固定” ▪ 実 技 試 験 Z:合否 Y:実技 試験 美⼤の⼊学試験における仮想例 合否Zで層別化すると、XとYの間に⾮因果的な連関 が⽣じ、X→Yの介⼊効果の推定にバイアスが⽣じる 合格者 Z=1 不合格者 Z=0 X:学⼒ 試験 100 100 学⼒試験 合 格 ラ イ ン
68. 内容の説明｜ バックドア基準とは？（概要版） “バックドア基準が満たされている” ▪ “バックドアパス”って何？ “開く／閉じる”ってどういうこと？ + (2) 処理X→結果Yの道がブロックされていない (1) 開きっぱなしのバックドアパスがない ≒ “ブロック”って何？
69. 超ざっくり説明｜ ”X→Yの道がブロックされてない” Zが中間点の場合は”固定”してはいけない ▪ Z Y X X Y Zを"固定"しなければ、何の問題も⽣じない
70. 超ざっくり説明｜ ”X→Yの道を閉じるな” 中間点をZを固定すると、Xからの「介⼊効 果」の流れ⾃体がブロックされてしまう ("overconditioning") ▪ 在来種Aを護るための 外来種Bの駆除の効果（仮想例） Z:外来種B の根絶 Y:在来種A の個体数 X:外来種B の駆除努⼒ 外来種Bの駆除努⼒ 在 来 種 A の 個 体 数 छBͷࠜઈʹ੒ޭ(Z=1) छBͷࠜઈʹࣦഊ(Z=0)
71. 内容の説明｜ バックドア基準とは？（概要版） “バックドア基準が満たされている” ▪ “バックドアパス”って何？ “開く／閉じる”ってどういうこと？ + (2) 処理X→結果Yの道がブロックされていない (1) 開きっぱなしのバックドアパスがない ≒ “ブロック”って何？
72. 具体例でおさらい｜３変量（X, Y, Z）の場合 3変量で⾒る”ほぼ バックドア基準” ▪ X Y Z ) X Y Z ) X Z Y ( ) Zを加えよ Zを加えるな ＊３変量の場合では、巷の”rule-of-thumb”的な変数選択 ルール以上のご利益はあまり感じられないかも
73. 具体例でおさらい｜ 4変量の場合 変数ベースではなく”道“の開閉がキモ ▪ このどの場合もバックドアパスは閉じられている ＊どのケースもXとYの"シンクロ"の原因となりうる上流側の流れが遮断されている X Y Z2 Z1 X Y X Y Z1 Z2 Z2 Z1
74. 4変量以上の場合｜本質は”道”にあり 4変量以上の”道”：合流点の例 ▪ Z1 X Z2 Y Z3 Z1 X Y Z3 X Z2 Y Z3
    X→Y   Z1 Z2 X Y Z3 Z1 Z2 バックドアパスが "開いている" バックドアパスは"閉まっている" バックドアパスの”道”がブロック されていることが肝要
75. 4変量以上の場合｜本質は”道”にあり 4変量以上の”道”：中間点の場合 ▪ X Y Z2 Z1 Z1, Z2は中間点 絶対にモデルに追加 したらダメ？ 中間点を⼊れるとXからの 因果効果⾃体が ブロックされてしまうのでダメ！ （これは3変量の場合と同様）
76. 4変量以上の場合｜本質は”道”にあり 4変量以上の”道”：中間点の場合 ▪ X Y Z1 Z2 全体としてX→Yの因果の"道"がブロックされないことが肝要 X Y X Y Z2 Z1 Z1 Z2
77. 4変量以上の場合｜まとめ X Y Z2 Z1 X Y X Y Z1 Z2 Z2 Z1 X Y Z1 Z2 X Y X Y Z2 Z1 Z1 Z2 バックドアパス をブロックしろ X→Yの道を ブロックするな もうほとんどバックドア基準 細かい論点をのぞけば 畢竟, 以下２つのメッセージに集約される ▪
78. "まとめ"としてのバックドア基準 教科書内の「バックドア基準」を解読してみよう ⿊⽊学(2017)『構造的因果モデルの基礎』 p99より引⽤ ▪
79. ステップ6 | "まとめ"としてのバックドア基準 教科書内の「バックドア基準」を解読してみよう ⿊⽊学(2017)『構造的因果モデルの基礎』 p99より引⽤ ▪ 1. XからZの任意の要素に有向道がない 2. GよりXから出る⽮線（X→）を除いたグラフにおいて、 ZがXとYを有向分離する 定義 3-5 バックドア基準（back door criterion） “逆”のケース（XがYの下流）ではない ⾮巡回的有向グラフGにおいてXはYの⾮⼦孫である。このとき、次 の２条件を満たす頂点集合Zは、（X, Y）についてバックドア基準を 満たすという 例：中間変数や下流の合流点 Xの下流の共変量をモデルに加えてはいけない バックドアパスがブロックされている
80. 4変量以上の場合｜まとめ（再掲） X Y Z2 Z1 X Y X Y Z1 Z2 Z2 Z1 X Y Z1 Z2 X Y X Y Z2 Z1 Z1 Z2 バックドアパス をブロックしろ X→Yの道を ブロックするな もうほとんどバックドア基準 ものすごい細かい論点をのぞけば 畢竟, 以下２つのメッセージに集約される ▪ 条件(1)に対応 条件(2)に対応
81. 例題その１｜どの変数を加えれば良い？ X Y Z1 Z2 Z3 バックドアパスの⾒極めに迷ったら、もし「その⼈⼯池にインクを ぶちまけたらどうなるか」を考えてみよう！
82. 例題その１｜どの変数を加えれば良い？ X Y Z1 Z2 Z3 Z1にインクをぶちまけると XとYの両⽅にインクは到達する バックドアパスの⾒極めに迷ったら、もし「その⼈⼯池にインクを ぶちまけたらどうなるか」を考えてみよう！
83. 例題その１｜どの変数を加えれば良い？ バックドアパスの⾒極めに迷ったら、もし「その⼈⼯池にインクを ぶちまけたらどうなるか」を考えてみよう！ X Y Z1 Z2 Z3 Z2にインクを ぶちまけると XとYの両⽅に インクは到達する
84. 例題その１｜どの変数を加えれば良い？ バックドアパスの⾒極めに迷ったら、もし「その⼈⼯池にインクを ぶちまけたらどうなるか」を考えてみよう！ X Y Z1 Z2 Z3 Z2からは こういう 経路もある
85. 例題その１｜どの変数を加えれば良い？ バックドアパスの⾒極めに迷ったら、もし「その⼈⼯池にインクを ぶちまけたらどうなるか」を考えてみよう！ X Y Z1 Z2 Z3 Z2からは この 経路もある
86. 例題その１｜どの変数を加えれば良い？ X Y Z1 Z2 Z3 バックドアパスの⾒極めに迷ったら、もし「その⼈⼯池にインクを ぶちまけたらどうなるか」を考えてみよう！ 全部で4つの バックドアパス がある
87. 例題その１｜どの変数を加えれば良い？ X Y Z3 Z2 Z1 Z1, Z2をブロックすると３つの バックドアパスが全て遮断される バックドア基準を満たす最⼩の変数セットは「Z1, Z2」 ＊「Z1, Z2, Z3」もバックドア基準を満たします
88. 例題その２｜どの変数を加えれば良い？ X Y Z1 Z4 Z3 Z2 Z5 Z6
89. 例題その２｜どの変数を加えれば良い？ X Y Z1 Z4 Z3 Z2 Z5 Z6
90. 例題その２｜どの変数を加えれば良い？ X Y Z4 Z3 Z2 Z5 Z6 Z1 Z1をブロックすると２つの バックドアパスが全て遮断される
91. 例題その２｜どの変数を加えれば良い？ X Y Z4 Z3 Z2 Z5 Z6 Z1 バックドア基準 を満たす最⼩の 変数セットは 「Z1」 ＊ここでZ5, Z6を⼊れると 新たなバックドアパスが 開いてしまう ＊Z4はそもそも 中間点なので⼊れてはダメ Z1をブロックすると２つの バックドアパスが全て遮断される
92. ⼀部が不明な場合｜どの変数を加えれば良い？ X Y Z1 Z2 この部分の 因果構造は不明
93. ⼀部が不明な場合｜どの変数を加えれば良い？ X Y Z1 Z2 この部分の 因果構造は不明 バックドアパスが ある場合、必ず Z1, Z2のどちらかを通る
94. ⼀部が不明な場合｜どの変数を加えれば良い？ X Z1 Z2 この部分の 因果構造は不明 Y 「Z1, Z2」を加えればバックドア基準を満たす Z1, Z2をブロックすれば バックドアパスは 遮断される ポイントとなる部分の適切な粒度の背景知識があればよい
95. ⼀部が不明な場合｜どの変数を加えれば良い？ X Z1 Z2 この部分の 因果構造は不明 Y Z1, Z2をブロックすれば バックドアパスは 遮断される ポイントとなる部分の適切な粒度の背景知識があればよい
96. image 例題その3| どの変数を加えれば良いのか？ • 構造の全てを知る必要は全くない：路線図の喩え つくば ⽴川 守⾕駅さえ封鎖すれば林は （電⾞では）家に帰れない
97. 例題その3| どの変数を加えれば良いのか？ • 構造の全てを知る必要は全くない：路線図の喩え image
98. 例題その3| どの変数を加えれば良いのか？ • 構造の全てを知る必要は全くない：路線図の喩え つくば 駅 ⽴川駅 守⾕ 駅 極度に 複雑な 鉄道網 必要なのは 「守⾕の先は⼀本線」 という"背景知識"のみ ポイントとなる部分の適切な粒度の背景知識があればよい TX
99. 例題その3| どの変数を加えれば良いのか？ 調整すべき変数の数は「処置Xを⽣じさせるシステム」と「結果 Yを⽣じさせるシステム」のデカップリングの程度に依存する （⼀般論として） X Y Z2 Z6 Z5 Z7 Z1 Z3 Z4 Z10 Z11 Z9 Z8
100. 例題その3| どの変数を加えれば良いのか？ 調整すべき変数の数は「処置Xを⽣じさせるシステム」と「結果 Yを⽣じさせるシステム」のデカップリングの程度に依存する （⼀般論として） X Y Z2 Z6 Z5 Z7 Z1 Z3 Z4 Z10 Z11 Z9 Z8 バックドアパスの成因となる変数=4 Xを⽣じさせる システム Yを⽣じさせる システム
101. 例題その3| どの変数を加えれば良いのか？ 調整すべき変数の数は「処置Xを⽣じさせるシステム」と「結果 Yを⽣じさせるシステム」のデカップリングの程度に依存する （⼀般論として） X Y Z2 Z6 Z5 Z7 Z1 Z3 Z4 Z10 Z11 Z9 Z8 バックドアパスの成因となる変数=11 Xを⽣じさせる システム Yを⽣じさせる システム
102. 例題その3| どの変数を加えれば良いのか？ 調整すべき変数の数は「処置Xを⽣じさせるシステム」と「結果 Yを⽣じさせるシステム」のデカップリングの程度に依存する （⼀般論として） X Y Z2 Z6 Z5 Z7 Z1 Z3 Z4 Z10 Z11 Z9 Z8 バックドアパスの成因となる変数=0 いわゆる”⾃然実験”のケース Xを⽣じさせる システム Yを⽣じさせる システム 調整すべき変数の数はかなり対象と⽂脈に依存する ＊あと実務上はバックドアパスが閉じているかはゼロイチの話ではない
103. 補⾜｜実験における無作為化との関係 X Y Z1 Z3 Z2 Z5 Z6 Z4 So many バックドアパス 実験はランダム⽣成機によりバックドアパスを閉じる ▪
104. 補⾜｜実験における無作為化との関係 X Y Z1 Z3 Z2 Z5 Z6 Z4 コイントス 実験は”ランダム⽣成機”によりバックドアパスを閉じる ▪
105. 補⾜｜傾向スコアとの関係 X Y Z1 Z3 Z2 Z5 Z6 Z4 So many バックドアパス 傾向スコアはバックドアパスをブロックする"合成変数" ▪
106. 補⾜ ｜傾向スコアとの関係 X Y Z1 Z3 Z2 Z5 Z6 Z4 e 傾向スコアで まとめてブロック 傾向スコアはバックドアパスをブロックする"合成変数" ▪ うまく適⽤できれば超強⼒な⼿法！
107. 補⾜ ｜差の差分析との関係 差の差分析は差分データ化してバックドアパスを消す ▪ X Y Z1 Z3 Z2 ΔX ΔY t t t t t X Y Z1 Z3 Z2 t' t' t' t' t' 時点t 時点tʼ 全要因の 時点間差分 をとる Z1, Z2, Z3が時間不変ならば ΔZ1, ΔZ2, ΔZ3はゼロ （項がキャンセルアウトされる） バックドアパス は消失
108. （再掲） ｜バックドア基準は「何について」の話？ ▪ 例えば、重回帰モデルの場合 モデルに追加された説明変数の組が「X→Y」についてバックドア基準を満 たすとき、重回帰分析から得られたXの偏回帰係数をそのまま「X→Yの介 ⼊効果」のバイアスのない推定量とみなせる 例えば、「シンプソンのパラドックス」(ここでは、どの変数で層別化す るかによって推定結果が変るケースの意味で⽤いる）が⽣じている場合 興味の対象となる「処理X→結果Y」についてバックドア基準を満たす変数 で層別化して解析すれば「X→Yの介⼊効果」をバイアスなく推定できる 具体的には、バックドア基準を満たすと： ＊データが適切に測定されており、かつ適切なモデルが適⽤されているという⼤前提での話です バックドア基準は介⼊効果推定での変数選択の基準 をもたらす
109. I. DAGとバックドア基準のまとめ “バックドア基準が満たされている” ▪ + (2) 処理Xの下流の変数が追加されていない (1) 開きっぱなしのバックドアパスがない ≒ 変数ベースではなく”道”の開閉がキモ ▪ ・たとえ森羅万象が因果関係で繋がっていても、因 果効果を識別可能とするために考えるべき局所的 な因果構造と変数の範囲を理論的に限定できる ・傾向スコアもバックドアパスを閉じるための強⼒ な⼿法（バックドア基準は理論的条件） 介⼊効果がバイアスなく推定できる ＊データが適切に測定されており、かつ適切なモデルが適⽤されているという⼤前提での話です
110. 林岳彦・⿊⽊学（2016） ｢相関関係」と「因果関係」は違います。これはよく知られています。ある要 因 X ともうひとつの要因 Y のあいだに高い相関が見られたからといって，それ らのあいだに因果的な関係があるとは限りません。一方で，そのような高い相関 を「因果関係」として解釈できる場合もたしかにあります。この辺りが難しいと ころです。もしあなたの同僚やクライアントが，あなたが作成した散布図を見て 「相関関係」と「因果関係」を明らかに混同した発言をしはじめたとしましょう。 このとき， 「この場合はこれこれこうだからこの相関関係は因果関係として解釈 できるんですよ/できないんですよ」と相手に向かって理路整然と説明するのは， それほど簡単なことではありません。こと因果関係の話になると，自分の頭の中 でその内容を整理するのも，その内容を相手に伝わるように説明するのも，なか なか難しいものです。 私たちの経験上，そんなときにとても役に立つのは，データの背後に想定して いる「因果構造(データ生成のメカニズム)」についての(分かる範囲での)簡単な ポンチ絵を丸と矢印で描いてみせることです[本稿ではそんなポンチ絵の例がたくさん 出てきます]。そして，そのようなポンチ絵を描いたあとに，その描かれた因果構 造が「あ ・ る ・ 特 ・ 定 ・ の ・ 条 ・ 件 ・ 群 ・ 」を満たしているかどうかを相手と共同で検討していき ます。多くの場合，その検討を通して「この相関関係は因果関係を示していると 解釈してよいのか？」や「相関関係を因果関係として解釈するためには本来はど のようなデータが必要なのか？」といった本質的な問いについて，より明確かつ 端的な議論ができるようになります。 相関と因果と丸と矢印のはなし はじめてのバックドア基準 林岳彦(国立環境研究所)･黒木学(統計数理研究所) [特集]因果推論 現実の課題に答える統計学 参考⽂献 | バックドア基準の解説原稿 (本⽇の元ネタ) 2016年6⽉発売「岩波データサイエンスvol. 3」因果推論特集号
111. 参考⽂献 | バックドア基準の解説 (ガチ勢向け)
112. 参考⽂献 | Pearlのやさしい本（昨年でた） Pearlが、⼀般向けに、 やさしく書いている！ （やさしく書かれているが けっこう深い／新しい 内容まで書いてある）
113. 参考⽂献 | 社会学者向け v IL Morgan and Winship (2014) 『Counterfactuals and Causal Inference: Methods and Principles for Social Research』 (2nd ed.) DAGと潜在結果モデル のhybrid
114. 1. 構造から⾒た統計的因果推論 ：DAGとバックドア基準の⼊⾨ 2. 差分から⾒た統計的因果推論 ：潜在結果モデル⼊⾨ 3. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在反応モデルのSCMによる表現 (3) 潜在反応モデルとDAGとの繋がり 4. 差分と構造の統合 --- Issueの総合的理解へ パートIの⽬次 0. 序：因果と相関と「対象の多様なありよう」
115. I. 解題 パートI：「構造」と「差分」の視点の統合 「対象を理解すること」にとって 統計的因果推論とは何か？ ⾮巡回有向グラフ(DAG) 潜在結果モデル 構造的因果モデル 𝒀 = 𝑿𝒀𝟏 + (𝟏 − 𝑿)𝒀𝟎 X Y Z1 Z2 Z3 Average Causal Effect = 𝑬[𝒀𝟏] − 𝑬[𝒀𝟎] 𝑋: = 𝑔: 𝑝𝑎 𝑋: , 𝜀: 𝑖 = 1, … , 𝑝 「差異からみた因果」 「構造からみた因果」
116. I. 解題 パートI：「構造」と「差分」の視点の統合 「対象を理解すること」にとって 統計的因果推論とは何か？ ⾮巡回有向グラフ(DAG) 潜在結果モデル 構造的因果モデル 𝒀 = 𝑿𝒀𝟏 + (𝟏 − 𝑿)𝒀𝟎 X Y Z1 Z2 Z3 Average Causal Effect = 𝑬[𝒀𝟏] − 𝑬[𝒀𝟎] 𝑋: = 𝑔: 𝑝𝑎 𝑋: , 𝜀: 𝑖 = 1, … , 𝑝 「差異からみた因果」 「構造からみた因果」
117. まず、個体レベルの因果効果を考える • 猫のノミの投薬治療の仮想例で考えていきます ノミに感染している猫へのノミ駆除剤の効果の例 ぴかそ 「駆除剤なし」 をX=なし、「駆除剤あり」 をX=ありとする 投薬により駆除までの⽇数Yを短くできるか？ ・ ・ “ぴかそ”における駆除剤の投薬の因果効果を 以下のように定義する “ぴかそ”に投薬した場合(X=1)の「ノミの駆除までの⽇数Y」 - “ぴかそ”に投薬しない場合(X=0)の「ノミの駆除までの⽇数Y」 たとえば投薬しない場合Y=14⽇、投薬した場合Y=4⽇なら、 投薬の因果効果は「-10⽇」となる
118. まず、個体レベルの因果効果を考える この差(-10⽇)が 個体の因果効果の定義 ぴかそ 投薬 ⽇数Y 14⽇ 4⽇ ॲཧT ݸମiͷ ൓ԠYi ͳ͠ ͋Γ • 実はこの因果効果の値を観測することは不可能である “ぴかそ”に投薬した場合(X=あり)の「ノミの駆除までの⽇数Y」 - “ぴかそ”に投薬しない場合(X=なし)の「ノミの駆除までの⽇数Y」 X
119. まず、個体レベルの因果効果を考える • 実はこの因果効果の値を観測することは不可能である “ぴかそ”に投薬した場合(X=あり)の「ノミの駆除までの⽇数Y」 - “ぴかそ”に投薬しない場合(X=なし)の「ノミの駆除までの⽇数Y」 同⼀の個体にはどちらか⼀⽅の処理しかできない 原理的に反応の「差」は観測できない→因果効果は観測できない “因果推論の根本問題” ぴかそに投薬した場合（X=あり） ぴかそに投薬した場合（X=なし） 観測 不可能 ↓ ↑ 観測 不可能 ॲஔT ݸମiͷ ൓ԠYi ͳ͠ ͋Γ ॲஔT ݸମiͷ ൓ԠYi ͳ͠ ͋Γ ぴかそ ぴかそ 14⽇ 4⽇ X X
120. 集団レベルの因果効果を考える ॲཧT ݸମiͷ ൓ԠYi ͳ͠ ͋Γ " # $ % & ' " # $ % & ' この差が 集団Aへの因果効果の定義 個体A, B, C, D, E, Fからなる集団Aへの平均効果を考える • 個体ではなくて集団なら”根本問題”を回避できる？ 集団Aに投薬した場合(X=あり)の「ノミの駆除までの平均⽇数Y」 - 集団Aに投薬しない場合(X=なし)の「ノミの駆除までの平均⽇数Y」 X
121. 集団レベルの因果効果を考える • 個体ではなくて集団なら”根本問題”を回避できる？ 集団Aに投薬した場合(X=あり)の「ノミの駆除までの平均⽇数Y」 - 集団Aに投薬しない場合(X=なし)の「ノミの駆除までの平均⽇数Y」 実際の処理では各個体にどちらか特定の処理を割り付けることなる 投薬なしのグループ（個体C, D, E) 投薬ありのグループ（個体A, B, D) ॲཧT ݸମiͷ ൓ԠYi ͳ͠ ͋Γ ॲཧT ͳ͠ ͋Γ " # $ % & ' " # $ % & ' " # $ % & ' " # $ % & ' " # $ % & ' 本来の因果効果 観察された 処理群間の差 各処理グループに含まれる個体は同⼀ではない →集団への因果効果も原理的に観測できない X X
122. 集団レベルの因果効果を考える • サブグループに分ければ”根本問題”を回避できる？ 集団A1に投薬した場合(X=あり)の「ノミの駆除までの平均⽇数Y」 - 集団A2に投薬しない場合(X=なし)の「ノミの駆除までの平均⽇数Y」 集団A 投薬なし（X=なし）のサブグループ （集団A1） 投薬あり（X=あり）のサブグループ （集団A2） ではどうしたら良いか？ 平均⽇数Yを観測 平均⽇数Yを観測 この差は観測可能であるが サブグループ間での「投薬の有無による差」と「背景因⼦の差」 が混ざってしまっていて区別がつかない 状況を潜在反応モデルで 整理してみよう
123. 潜在結果モデルという枠組み • 処置（投薬）×反応（⽇数）の組合せを整理してみる X=無のサブグループ X=有のサブグループ もしX=無 のときの 反応 (i)の「駆除までの⽇数」 (ii)の「駆除までの⽇数」 もしX=有 のときの 反応 (iii)の「駆除までの⽇数」 (iv)の「駆除までの⽇数」 ＠現実世界 ＠反事実的可能世界 ＠反事実的可能世界 ＠反事実的可能世界 因果推論とは での反応と ＠現実世界 観測可能な での反応の「差」の計算 ＠現実世界 観測不可能な ＠反事実的可能世界 の世界も扱える枠組みがほしい！
124. 潜在結果モデルという枠組み • Rubinの潜在結果モデルの考え⽅ 処理Tに対して、各個体𝑖は⼀意な（決定論的に定まる） 潜在結果𝑌(𝑖)!"($)をもつ 𝒀(ぴかそ)𝐢𝐟(𝑿% ) ぴかそ だり もしも「ぴかそ」が 「駆除剤あり」のときの 「⽇数」の潜在反応 ノミに感染している猫へのノミ駆除剤の効果の例 「駆除剤なし」 をX=無、「駆除剤あり」 をX=有、「駆除までの⽇数」をYとする 𝒀(だり)𝐢𝐟(𝑿% ) もしも「だり」が 「駆除剤なし」のときの 「⽇数」の潜在反応 有 無 ＊このifを明⽰する表記法は⼀般的ではないが 分かりやすさ優先して本⽇はこれを使⽤
125. 潜在結果モデルという枠組み • このモデルの何が嬉しいかと⾔うと、因果効果を 「処理群間での潜在反応の差」として定式化できる 潜在結果を⽤いた「ぴかそ」への投薬の因果効果の定義 𝒀 ぴかそ 𝐢𝐟(𝑿% ) − 𝒀 ぴかそ 𝐢𝐟(𝑿% ) ぴかそ 現実世界内ではこの両者を同時に観測できない “因果推論の根本問題” 同じ個体に同時に異なる処理はできない 有 無 ＊猫の名前が「しゅれーでぃんがー」の場合は別途要検討
126. ࣮ࡍͷॲஔT જࡏ൓Ԡ   Y(i) if(T=0) if(T=1) Y(i) if(T) Y 潜在結果モデルという枠組み • 潜在モデルの枠組みを⽤いると、”因果推論の根本 問題”を「⽋測」の枠組みで捉えることができる Y (i)if(T=0) if(T=1) ɿݸମ
     i ʹॲཧT=0Λ༩͑ͨͱ͖ͷજࡏ൓Ԡ ɿݸମ
     i ʹॲཧT =1Λ༩͑ͨͱ͖ͷજࡏ൓Ԡ Y (i) if(T=1) Y (i) Y (i)if(T=0) - ݸମ
     J
     ʹ͓͚ΔҼՌޮՌͷఆٛɿ この差が個体iの 因果効果の定義 ぴかそ 投薬 駆除までの⽇数 “⽋測” “⽋測” _ _ / 有 無 有 無 有 無 有 無 X X X X X X X X X X 結果 有 無
127. 潜在結果モデルという枠組み • ねこ集団Aへの「平均因果効果」を潜在反応を⽤い て定義すると： 𝑬[𝒀 集団𝑨 𝐢𝐟(𝑿% )] − 𝑬[𝒀(集団𝑨) 𝐢𝐟(𝑿% )] 集団レベルでもこの両者は同時に観測できない “因果推論の根本問題” 有 無
128. 潜在結果モデルという枠組み • 何を知る必要があるのかを潜在反応モデルの枠組み でもういちど整理してみる 計算したいのは（たとえば）因果効果 (iv) – (ii) Average Treatment Effect of the Treated (ATT) X=無のサブグループ X=有のサブグループ X=無の 潜在反応 𝑬[𝒀 集団𝑨|𝑿 = 𝐢𝐟(𝑿) )] 𝑬[𝒀 集団𝑨|𝑿 = 𝐢𝐟(𝑿) )] X=有の 潜在反応 𝑬[𝒀 集団𝑨|𝑿 = 𝐢𝐟(𝑿) )] 𝑬[𝒀 集団𝑨|𝑿 = 𝐢𝐟(𝑿) )] そもそも観測不能（反事実） そもそも観測不能（反事実） 観測可能 観測可能 (i) (ii) (iii) (iv) けれども、(ii)は反事実なので観測＆計算できない！ 反事実項の値を知るにはどうすればよい？ 無 無 有 有 無 無 有 有
129. 潜在結果モデルからみる無作為化 • 解決策：無作為化により反事実項を観察項で代替 X=無のサブグループ X=有のサブグループ X=無の 潜在反応 𝑬[𝒀 集団𝑨|𝑿 = 𝟎 𝐢𝐟(𝑿)𝟎)] 𝑬[𝒀 集団𝑨|𝑿 = 𝟏 𝐢𝐟(𝑿)𝟎)] X=有の 潜在反応 𝑬[𝒀 集団𝑨|𝑿 = 𝟎 𝐢𝐟(𝑿)𝟏)] 𝑬[𝒀 集団𝑨|𝑿 = 𝟏 𝐢𝐟(𝑿)𝟏)] そもそも観測不能（反事実） そもそも観測不能（反事実） 観測可能 観測可能 (i) (ii) (iii) (iv) サンプル集団 処理を無作為に割付 因果効果 (iv) – (ii) を (iv) – (i) として計算可能 無作為割付→サブグループの間で背景因⼦が揃うと期待できる = 反事実である(ii)を観察可能な(i)で代替可能 「同じ集団」ではないが「⼗分に似た集団」として =
130. 潜在結果モデルから⾒る無作為化 • 同じこと違う表現でもういちど説明すると： 𝑬[𝒀 集団𝑨 𝐢𝐟(𝑿A𝟏)] − 𝑬[𝒀(集団𝑨) 𝐢𝐟(𝑿A𝟎)] 同⼀の集団に異なる処理を与えることはできない 𝑬[𝒀 集団𝑨 𝐢𝐟 𝑿)𝟏 |𝑿 = 𝟏] − 𝑬[𝒀(集団𝑨) 𝐢𝐟 𝑿)𝟎 |𝑿 = 𝟎] 無作為に（たとえばコイン投げで）処理Tを割り付ける 潜在反応モデルによる集団Aへの因果効果の定義 「同⼀ではないけど⼗分に背景が揃っている集団」に 「異なる処理X=0/X=1」を与えることができる ＝因果効果の定義に⼗分に近い値を推定できると期待できる
131. 補⾜｜傾向スコア法 X=0を受けた群 X=1を受けた群 X=0の潜 在反応 𝑬[𝒀𝟎 ねこ集団𝑨|𝑿 = 𝟎 ] 𝑬[𝒀𝟎 ねこ集団𝑨|𝑿 = 𝟏 ] X=1の潜 在反応 𝑬[𝒀𝟏 ねこ集団𝑨|𝑿 = 𝟎 ] 𝑬[𝒀𝟏 ねこ集団𝑨|𝑿 = 𝟏 ] 共変量 の分布 そもそも観測不能（反事実） そもそも観測不能（反事実） 観測可能 観測可能 (i) (ii) (iii) (iv) = = 潜在反応と処理Xが条件付き独⽴→ 𝑌'(()) 𝑋 | 𝑒 ⊥ ⊥ eでマッチング・重みづけ 傾向スコアe 共変量 両群での共変量の分布のインバランスをpost-hocに補正するための バランシングスコアeを共変量から計算する 因果効果 (iv) – (ii) を (iv) – (i) として計算可能 条件付き無視 可能性条件
132. の条件を満たすことで、観察値 に基づき「反事実」の項も計算可能 （尚、バックドア基準を満たせば上記条件は満たされる） 潜在反応モデルと関数的因果モデル｜ Rubinの潜在結果モデルの特徴（まとめ） 因果効果を「異なる処理に対する潜在結果の 差分」として定義する 因果効果の定義式には「反事実」の項が含まれる 𝑌!"($) 𝑋|𝑍 ⊥ ⊥ 潜在反応の中⾝（”why it works”）はブラック ボックスでも構わない ＊潜在反応がブラックボックスであるのは理論的 な必然というわけではない（次の話につづく） 尚、処理Xが連続量のときの 表現が苦⼿という側⾯もある 中⾝に関わらず適⽤できるので 「推定プロトコル」への⾒通しがよい たとえばRCTとか傾向スコアなどで
133. 1. 構造から⾒た統計的因果推論 ：DAGとバックドア基準の⼊⾨ 2. 差分から⾒た統計的因果推論 ：潜在結果モデル⼊⾨ 3. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在結果モデルのSCMによる表現 (3) 潜在結果モデルとDAGとの繋がり 4. 差分と構造の統合 --- Issueの総合的理解へ パートIの⽬次 0. 序：因果と相関と「対象の多様なありよう」
134. 参考⽂献｜⿊⽊（2017） ここからの話はこの本の 第六章が元ネタです 本⽇の説明は やや「超訳」的に なりますので 正確な議論については ぜひ本書を ご参照ください
135. 因果グラフと関数的因果モデル｜ 構造的因果モデルSCM（a.k.a.関数因果モデル）とは 因果構造を関数の親⼦関係の連なりの形で記述したもの 関数因果モデルのフォーマルな定義 ⾮巡回的有向グラフGとその頂点に対応する確率変数の集合V = {X1,…Xp} が与えられている．グラフGが確率変数間の関数関係を ⿊⽊(2017), p70より引⽤ になる形に規定し、確率変数がこの関数関係にしたがって⾃律的でかつ 定常的に⽣成されるとき、Gを因果ダイアグラムという．ここに、錯乱 項ε1,…,εpは互いに独⽴であるとする．また、pa(Xi)は因果ダイアグラムG におけるXiの親全体からなる変数集合であり、GにおけるXiの直接的原因 (direct cause)と解釈される． 上式を広い意味で物理的なデータ⽣成過程のモデルとして捉えるとき、上式は 関数因果モデルあるいは、構造⽅程式モデル、構造的因果モデルと呼ばれる. 𝑋! = 𝑔! 𝑝𝑎 𝑋! , 𝜀! 𝑖 = 1, … , 𝑝 これは等式というより代⼊式 ベイジアンネットワークとは異なり「因果」概念を内包する Structural Causal Model
136. 因果グラフと関数的因果モデル｜ 構造的因果モデル（SCM）とは 因果ダイアグラムはSCMの視覚的表現といえる （実は後者の⽅が表現できるものの幅は広い） 因果構造を関数の親⼦関係で記述したもの GUIとCUIのようなイメージ（CUIが本体） Y X Z1 Z2 Z3 特定の関数型を特に 念頭に置く話ではない 𝑌 = 𝑔" (𝑋, 𝑍1) 𝑋 = 𝑔# (𝑍2, 𝑍3) 𝑍1 = 𝑔$% (𝑍2) 𝑌 = ℎ" (𝑍2, 𝑍3) 実は「Pearlの体系=グラフィカルモデル」と いうのは適切な理解ではない！ こうも書けるし このように 外⽣変数の関数 としても書ける 詳細は⿊⽊(2017) (錯乱項は省略）
137. 1. 構造から⾒た因果推論 ：DAGとバックドア基準の⼊⾨ 2. 差分から⾒た因果推論 ：潜在反応モデル⼊⾨ 3. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在反応モデルのSCMによる表現 (3) 潜在反応モデルとDAGとの繋がり 4. 差分と構造の統合 --- Issueの総合的理解へ パートIの⽬次
138. 潜在反応モデルと関数的因果モデル｜ 潜在結果モデルを関数的因果モデルで表現する 𝑌'(())(ぴかそ) = 𝑓(𝑖𝑓(𝑋), 𝑍1, 𝑍2, 𝑍3, . . , 𝑍𝑛) 「ぴかそ」の特性を表す全変数 𝑌'(()%*)(ぴかそ) = 𝑓(𝑋 = 0, 𝑍1, 𝑍2, . . , 𝑍𝑛) 𝑌'(()%+)(ぴかそ) = 𝑓(𝑋 = 1, 𝑍1, 𝑍2, . . , 𝑍𝑛) 潜在結果の”中⾝”を関数𝑓と変数𝑍で表現している ぴかそ 潜在結果は概念的にはその個体がもつ 諸特性の決定論的関数で表現できる 𝑖𝑓 𝑋 = 0, 𝑖𝑓 𝑋 = 1, 潜在結果の脱ブラックボックス化 SUTVA条件のより明⽰的な検討へ繋がる→パート2へ
139. 潜在反応モデルと関数的因果モデル｜ 例えば全特性が近似的に等しい個体がいれば… 𝑌'(()%*)(たつや) = 𝑓(𝑋 = 0, 𝑍1, 𝑍2, 𝑍3, . . , 𝑍𝑛) 𝑌'(()%*)(かずや) = 𝑓(𝑋 = 0, 𝑍1, 𝑍2, 𝑍3, . . , 𝑍𝑛) もし「たつや」「かずや」という双⼦がいて もし関数𝑓および変数𝑍が近似的に等しければ 𝑌#$(&'() たつや − 𝑌#$(&'*) たつや ≒ 𝑌#$(&'() たつや − 𝑌#$(&'*) かずや これは観測可能！ 「関数𝑓と特性𝑍」から識別可能性（バイアスの ない因果効果の推定が可能か）を議論しうる 観測不可能な因果効果 少数事例研究がやってることはこれに近い？
140. 潜在反応モデルと関数的因果モデル｜ 「たつや」が「かずや」と 代替可能とは⼀体 なにを意味しうるのか →第III部へ III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 ⼀般に統計解析／統計的因果推論はある種の”交換可能性”を想定しているが、 現実にはどこまでいっても完全な代替などできない ＊可能世界と反事実と代替の（不）可能性を巡る物語の例 image
141. 潜在反応モデルと関数因果モデルと因果グラフ｜ 集団Aへの平均因果効果を考える 𝐸[𝑌'(())(集団𝐴)] = 𝐸[𝑓(𝑖𝑓(𝑋), 𝑍1, 𝑍2, . . . , 𝑍𝑛)] 平均因果効果＝𝐸[𝑌,-($).) ねこ集団𝐴 − 𝐸[𝑌,-($)/) ねこ集団𝐴 ] = 𝐸[𝑌,-($).) ねこ集団𝐴|𝑋 = 1 ] − 𝐸[𝑌,-($)/) ねこ集団𝐴|𝑋 = 0 ] 𝐸[𝑌:C(DA=) 集団𝐴|𝑋 = 1 ] = 𝐸[𝑓 𝑋, 𝑍1, 𝑍2, . . , 𝑍𝑛 𝑋 = 1 ] 𝐸[𝑌:C(DAE) 集団𝐴|𝑋 = 0 ] = 𝐸[𝑓(𝑋, 𝑍1, 𝑍2, . . , 𝑍𝑛|𝑋 = 0)] もし異なる処理Xを受けたサンプル集団間で 関数fおよび特性Zの分布が同等とみなせる場合 異なる処理を受けた集団間で 潜在反応の分布が異ならない 観測値から因果効果を計算可能
142. 潜在反応モデルと関数因果モデルと因果グラフ｜ 「集団間（もしくは個体間）での 関数fおよび特性Zの分布の同等性」 構造的 類似性 特性的 類似性 ⼀般に：潜在結果の類似性を決める⼆つの要素 𝑌'(()%,) 集団𝐴 = 𝑓(𝑖𝑓(𝑋 = 𝑥), 𝑍1, 𝑍2, . . , 𝑍𝑛) 𝑌'(()%,) 集団𝐵 = 𝑓(𝑖𝑓(𝑋 = 𝑥), 𝑍1, 𝑍2, . . , 𝑍𝑛) 両者が⼗分に近ければ集団Aの結果を集団Bの結果で代替できる （と⼀般に解釈される）
143. 1. 構造から⾒た因果推論 ：DAGとバックドア基準の⼊⾨ 2. 差分から⾒た因果推論 ：潜在反応モデル⼊⾨ 3. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在反応モデルのSCMによる表現 (3) 潜在反応モデルとDAGとの繋がり 4. 差分と構造の統合 --- Issueの総合的理解へ パートIの⽬次
144. 潜在反応モデルと関数因果モデルと因果グラフ｜ 潜在反応モデルをDAGで表現してみると 畢竟、潜在反応モデルと関数因果モデルと バックドア基準の話はシームレスに繋がっている Y X Z1 Z2 Z3 バックドア基準に基づき、たとえばこの潜在反応モデ ルの例で調整する必要があるのはZ1だけ 潜在反応に実際のXの値が与えられて観測値Yとなる 𝑌'(()) 𝑌 = 𝑓 𝑌,-($), 𝑋 = 𝑋𝑌,-($).) + (1 − 𝑋)𝑌,-($)/) たとえばXが0/1のときはこう書ける 𝑌#$(&) = 𝑓(𝑖𝑓 𝑋 , 𝑍1, 𝑍3) 𝑌:C(D) 𝑋|𝑍1 ⊥ ⊥ 条件付き無視可能性
145. 本⽇の内容 I. 構造から⾒た因果推論 ：DAGとバックドア基準の⼊⾨ II. 差分から⾒た因果推論 ：潜在反応モデル⼊⾨ III. 構造的因果モデル(SCM)による統合 (1) 因果グラフのSCMによる表現 (2) 潜在反応モデルのSCMによる表現 (3) 潜在反応モデルとDAGとの繋がり IV. 差分と構造の統合 --- Issueの総合的理解へ
146. Y X Z1 Z2 Z3 𝑌 = 𝑔" (𝑋, 𝑍1) 𝑋 = 𝑔# (𝑍2, 𝑍3) 𝑍1 = 𝑔$% (𝑍2) 𝑌 = ℎ" (𝑍2, 𝑍3) こうも書けるし こうも書ける (1) 構造的因果モデルから帰納的／演繹的因果 推論の全体を⾒る 構造的因果モデルは⾮常に広いクラスのモデルを含む
147. 構造的因果モデルから帰納的／演繹的因果推論の全体を⾒る 𝑌 = ℎ" (𝑍2, 𝑍3) こうも書ける この形、⾮常に⼀般的である 𝑇𝐸𝑀𝑃 = 𝑚𝑎𝑖𝑛 (モデルパラメータベクトル) たとえばスパコンで回すような気候変動シミュレー ションモデル（What if/反事実的計算に⽤いられる）もこんな 形で書きうる 構造的因果モデルは⾮常に広いクラスのモデルを含む
148. 潜在反応 モデル メカニズムが ブラックボックス であるたぐいの 物理シミュレー ションモデル 因果効果推定アプローチの連続的な理論観（私⾒） 中間変数を介 す２段階回帰 under フロントドア 基準 構造推定＠ 経済学 X以外の共変量を バランシングして いくたぐいの 𝑌 = 𝑓 (処理𝑋, 外⽣変数𝑧) 重回帰 under バックドア 基準 RCT, 準実験 例えばスパコンで回す ような気候変動モデル 関数fと変数zについての モデル式の正当性に依存しない ただし、関数fと変数zの理論が なければ移設可能性の把握が困難 関数fと変数zについての モデル式の正当性に依存する 関数fと変数zの同等性の検討 により移設可能性を議論可能 祖先変数を介 す操作変数・ ⼆段階SLS 構造的因果モデルから帰納的／演繹的因果推論の全体を⾒る Full-of-Model Model-free X→Yの祖先-⼦孫 系列の効果伝搬を 演繹計算・集計し ていくたぐいの
149. (2) 「Issueの総合的理解」における因果推論 ༗ػਫۜ ͷੜ੒ ਫ࢈෺ ࢢຽ ࢢ৔ 例：有機⽔銀汚染に対する環境リスク研究者の"頭の中” ޻৔Ͱͷ ੜ࢈աఔ ړࢣ ഉਫ ւͰͷ ੜ෺ೱॖ ւͰͷ Ԛછ஝ੵ “Issue/problem-driven”の解析者と、”Method- driven”の解析者の”頭の中”は異なる Issue-drivenの解析者の 多くは常に頭の中に 何らかの "因果モデル"をもつ
150. (2) 「Issueの総合的理解」における因果推論 150 どこに介⼊しうるのか／もし介⼊したらどうなるのかについてのアイデ アが常に頭の中をぐるぐるしている ༗ػਫۜ ͷੜ੒ ਫ࢈෺ ࢢຽ ࢢ৔ ޻৔Ͱͷ ੜ࢈աఔ ړࢣ ഉਫ ւͰͷ ੜ෺ೱॖ ւͰͷ Ԛછ஝ੵ ⼯場の⽣産過程を 変えれば有機⽔銀 の⽣成を抑えられ るのでは？ 排⽔処理を変えれ ば有機⽔銀を除去 できるのでは？ 禁漁により⼈への 曝露を防げるのは？ 浚渫により汚染を 低減できるので は？ 市場への流通を⽌めれば 市⺠への曝露は防げるのでは？ 妊婦や⼦供などへの注意喚起により ハイリスクグループへの曝露を減ら せるのでは？ 市場への流通を⽌める と⾃家消費により曝露 が増えるかも？ 例：有機⽔銀汚染に対する環境リスク研究者の"頭の中”
151. 「差異」と「構造」の繋がりが⾒えないと・・ 151 • 頭の中の"因果モデル"と「介⼊効果推定」の関係性は？ 頭の中の"因果モデル" ༗ػਫۜ ͷੜ੒ ਫ࢈෺ ࢢຽ ࢢ৔ ޻৔Ͱͷ ੜ࢈աఔ ړࢣ ഉਫ ւͰͷ ੜ෺ೱॖ ւͰͷ Ԛછ஝ੵ （政策的対応を念頭においた） 介⼊効果推定 市場= 流通禁⽌ 市場= 流通あり 市⺠の 健康 市⺠の 健康 ここの差分 ＝介⼊効果 どちらかに思考のピントが合っているとき、もう⼀⽅はぼやけてしまう
152. 「差異」と「構造」の繋がりが⾒えないと・・ 152 • 頭の中の"因果モデル"と「介⼊効果推定」の関係性は？ 頭の中の"因果モデル" ༗ػਫۜ ͷੜ੒ ਫ࢈෺ ࢢຽ ࢢ৔ ޻৔Ͱͷ ੜ࢈աఔ ړࢣ ഉਫ ւͰͷ ੜ෺ೱॖ ւͰͷ Ԛછ஝ੵ （政策的対応を念頭においた） 介⼊効果推定 市場= 流通禁⽌ 市場= 流通あり 市⺠の 健康 市⺠の 健康 ここの差分 ＝介⼊効果 どちらかに思考のピントが合っているとき、もう⼀⽅はぼやけてしまう
153. 153 両者にピントを合わせて思考できるようになってくる 頭の中の"因果モデル" ༗ػਫۜ ͷੜ੒ ਫ࢈෺ ࢢຽ ࢢ৔ ޻৔Ͱͷ ੜ࢈աఔ ړࢣ ഉਫ ւͰͷ ੜ෺ೱॖ ւͰͷ Ԛછ஝ੵ （政策的規制措置等を念頭においた） 介⼊効果推定 do(市場= 流通禁⽌） do(市場= 流通あり） 市⺠の 健康 市⺠の 健康 ここの差分 ＝介⼊効果 両者の関係が⾒えるぞ・・・！ 「差異」と「構造」の繋がりが⾒えてくると・・・
154. バックドア基準 構造的因果モデル そのときふと ⾃分の⾜元を ⾒てみてみると・・ 両者にピントを合わせて思考できるように 「差異」と「構造」の繋がりが⾒えてくると・・・ DAG
155. 変数を概念的に吟味するとはどういうことか｜ 統計的因果推論と諸研究の関連性のイメージ 因果 モデル 関数fと変数Zのあり よう⾃体を分析対象 とした質的/量的研究 介⼊効果 推定 統計的 因果推論 分析概念の／関数fと変数Zの モデルの⾼解像度化 バイアスの減少 Target/SUTVAの明確化 深化 効果の定量的な実証 量的な重要度の明確化 介⼊効果がいかに／どの程度⽣じるかについての “Sufficiently deep explanation”を社会に提供できる Morgan and Winship (2014) 「差がありました」の先の 「どうすればいいの？」に答えたい！ 例えば、個別ケースのプロセストレーシング
156. 変数を概念的に吟味するとはどういうことか｜ 統計的因果推論で因果グラフを⽤いることの効⽤ ・ 調整すべき変数についての論理が得られる ・ 量的（統計的）研究と質的研究の橋渡しとなりうる ・ 調整すべき変数を絞り込むことができる ・ 感度分析・媒介効果分析への橋渡しとなりうる ・ 調整すべき変数の抜け漏れを防げる ・ 因果の図解⾃体が関係者とのやりとりに役⽴つ ・ 移設可能性についての議論も深まりやすくなる（移設可能性は⽣ 成メカニズムの安定性や類似性にも⼤きく依存するため） ・ 識別可能性の前提条件の成⽴を検証するための視覚化 “Assumptions visualization tool” ・ 結果・処理・共変量の変数の概念的妥当性の吟味に役⽴つ “Issueの総合的理解”への⾒通しが開きやすい
157. バックドア基準は、「バイアスのない因果効果推 定」のためにはどの変数を考慮するべきかの理論的 条件を背景にある因果構造を基に与える パートIの話の全体まとめ ▪ 潜在反応モデルは「潜在反応の差分」として因果効 果を定式化し、因果効果推定の解析プロトコルへの 明瞭な⾒通しを与える ▪ 構造的因果モデルは⾮常に広いクラスの因果効果の 推定アプローチを包含しうるものであり、issueの総 合的理解へ向けた統計的／理論的／理念的な枠組み の基盤となりうる→混合研究指向との親和性 ▪
158. 導⼊｜統計的推測の主要な3⽬的（再掲） ⽬的の 分類 やりたいことの内容 事例のイメージ （⼤学の留年の例で） 理念的な理想をなす 解析法のイメージ しばしば 現実的な 解析法 予測 現在得られている情報 から、将来／未知の値 ／状況を予測したい （要因への介⼊は伴わない場合） 多数の変数の情報から、留年 リスクの⾼い学⽣を⾒つけた い（→その後⽀援につなげる など） 機械学習 回帰分析！ 説明 要素間の関連性を説明 したい なぜ学⽣が留年するかを説明 したい（要因を探索し理解／ 対策につなげたい） フルスタックの 構造⽅程式 （含む因⼦分析等） 回帰分析！ 因果／介⼊ 効果の推定 ある特定の要因に介⼊ したときの、ある特定 の要因への効果を知り たい 特定のプログラムの導⼊が留 年を防⽌する効果を知りたい ランダム化⽐較試験 回帰分析！ ⼀⾒同じ「回帰分析」でも、⽬的に応じて係数の 解釈のあり⽅も、重要となる仮定群も異なる （異論は認める） ＊ ＊互いに排反な分類という意味ではないので注意 ↑ パート1では「因果効果の推定」が主テーマ ▪
159. 導⼊｜統計的推測の主要な3⽬的（再掲） ⽬的の 分類 やりたいことの内容 事例のイメージ （⼤学の留年の例で） 理念的な理想をなす 解析法のイメージ しばしば 現実的な 解析法 予測 現在得られている情報 から、将来／未知の値 ／状況を予測したい （要因への介⼊は伴わない場合） 多数の変数の情報から、留年 リスクの⾼い学⽣を⾒つけた い（→その後⽀援につなげる など） 機械学習 回帰分析！ 説明 要素間の関連性を説明 したい なぜ学⽣が留年するかを説明 したい（要因を探索し理解／ 対策につなげたい） フルスタックの 構造⽅程式 （含む因⼦分析等） 回帰分析！ 因果／介⼊ 効果の推定 ある特定の要因に介⼊ したときの、ある特定 の要因への効果を知り たい 特定のプログラムの導⼊が留 年を防⽌する効果を知りたい ランダム化⽐較試験 回帰分析！ ⼀⾒同じ「回帰分析」でも、⽬的に応じて係数の 解釈のあり⽅も、重要となる仮定群も異なる （異論は認める） ＊ ＊互いに排反な分類という意味ではないので注意 ↑ パート1では「因果効果の推定」が主テーマ ▪ 構造的因果モデルはこれらの⽅向性の 異なる世界線を架橋しうる点が嬉しい
160. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
161. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
162. 2. そのXは本当にXか：概念的妥当性を吟味する 3. その因果効果はどこまで”拡⼤解釈”できるか パートIIの⽬次 ・ その1：媒介経路の観点から ・ その2：要因間の依存性の観点から ・ その3：⽂脈依存性の観点から ・ その4：測定の観点から ・ ターゲット妥当性の分解 1. 可能世界への”チャネル”としてのSUTVA条件 ・ 事例：壊⾎病と”シトラス” ・ 事例：マシュマロテストの”拡⼤解釈”
163. 2. そのXは本当にXか：概念的妥当性を吟味する 3. その因果効果はどこまで”拡⼤解釈”できるか パートIIの⽬次 ・ その1：媒介経路の観点から ・ その2：要因間の依存性の観点から ・ その3：⽂脈依存性の観点から ・ その4：測定の観点から ・ ターゲット妥当性の分解 1. 可能世界への”チャネル”としてのSUTVA条件 ・ 事例：壊⾎病と”シトラス” ・ 事例：マシュマロテストの”拡⼤解釈”
164. SUTVA条件とは SUTVA条件の定義 個体 i を特徴づける潜在的な結果{Yi(1), Yi(0)}が次の2条件を満⾜するとき、 それはSUTVA（Stable Unit Treatment Value Assumption）条件を満⾜する という 岩崎学『統計的因果推論』, p78より引⽤ （表記法を⼀部編集） ▪ Rubinの潜在結果モデルでの因果効果の定義 はSUTVA条件を前提とする (a) 個体 i の潜在的な結果は、他の個体の受ける処置に依存しない (b) 個体 i に対する処置は1通りに定まる 潜在結果を⽤いた「ぴかそ」への投薬Xの因果効果の定義 𝒀 ぴかそ 𝐢𝐟(𝑿% ) − 𝒀 ぴかそ 𝐢𝐟(𝑿% ) 有 無 Stable Unit Treatment Value Assumption
165. (hidden versionsがあると) SUTVA条件にあてはまらない SUTVA条件とは ▪ 「個体 i に対する処置は1通りに定まる」とは？ （no hidden version of treatment） “処置＝ワクチン接種” アストラゼネカ？ モデルナ？ “処置”を明確に定義できないので 因果効果がill-definedになる シノファーム？ スプートニクV？ ファイザー？
166. SUTVA条件とは ▪ 「個体 i の潜在的な結果は、他の個体の受 ける処置に依存しない」とは？ 例えば、他の個体のワクチンの効果はその個体の 感染を防ぐだけではなく、（その個体がキャリア にならないことで）⾃分の感染も防ぐ効果をもつ 「個体iの潜在的な結果が、他の個体の受ける処置に 依存する」→ 個体iへの処置の効果が⼀義に定まらない （他の個体が受けた処置によって⾃分の結果が変わるため） SUTVA条件にあてはまらない
167. SUTVA条件とは ▪ なぜSUTVA条件が重要なのか？ 個体i@反事実世界A どの反事実世界との コントラストを取れば良い？ 内部機構を⼀切もたないRubinの潜在結果モデルの枠組み では、ある特定の可能世界とのコントラストを取ること （なぜこの世界がその可能世界と接続されるのか）を正当 化する内在的な論理がない （観測不可能） 個体i@この世界 （観測可能） 𝒀 𝒊 個体i@反事実世界B 個体i@反事実世界C 個体i@反事実世界D 3ヶ⽉⽔を誰もやらず、花が枯れた世界 林が⽔を毎⽇あげた世界？ or 林が⽔を⽉⼀であげた世界？ orトランプが毎⽇来⽇して⽔をあげた世界？
168. SUTVA条件とは ▪ “まずSTUVAあれ！” 個体i@この世界 （観測不可能） “処置X”という チャネルが ⼀意に定まる ことで 特定の可能世界との 接続がdefineされる Stable Unit Treatment Value 𝒀 𝒊 𝐢𝐟(𝑿) 「なぜ」とかじゃなくて 個体i@反事実世界A 個体i@反事実世界B 個体i@反事実世界C 個体i@反事実世界D （観測可能） 𝒀 𝒊 ⽔がなく、花が枯れた世界 X = 林が週⼀で⽔をあげる 林が⽔を週⼀で⽔をあげた世界 Rubinの潜在結果モデルではSUTVAが基底的な （それ以上その体系内では深掘りできない） 仮定として要請される
169. SUTVA条件とは ▪ STUVA条件はどう吟味するのか “概念を質的に深掘りしていくこと”と 因果効果推定の関係を考えてみる 「概念に(hidden) multiple versionがない」 かどうかを判断するには、ほんらい、質的な ドメイン知識が必須である Rubinの潜在結果モデルでは、その基底的 な仮定の妥当性を質的知識に完全に外部化 （丸投げ）している
170. 2. そのXは本当にXか：概念的妥当性を吟味する 3. その因果効果はどこまで”拡⼤解釈”できるか パートIIの⽬次 ・ その1：媒介経路の観点から ・ その2：要因間の依存性の観点から ・ その3：⽂脈依存性の観点から ・ その4：測定の観点から ・ ターゲット妥当性の分解 1. 可能世界への”チャネル”としてのSUTVA条件 ・ 事例：壊⾎病と”シトラス” ・ 事例：マシュマロテストの”拡⼤解釈”
171. 概念の”質的”吟味：媒介経路の観点から 薬の処⽅X→健康指標Yの因果効果は？ ▪ ॲํ9 ݈߁ ࢦඪ: いや、そもそも 「処⽅」って何？ ॲํ9͋Γͳ͠ ΛϥϯμϜʹׂ෇ ॲํ9ͷҼՌޮՌΛ όΠΞεͳ͘ਪఆՄೳ
172. 概念の”質的”吟味：媒介経路の観点から 薬の処⽅X→健康指標Yの因果効果は？ ▪ ॲํ9 ݈߁ ࢦඪ: ༀͷ੒෼ͷ ༀཧ࡞༻ 元々知りたいのは 薬理作⽤の効果
173. 概念の”質的”吟味：媒介経路の観点から ▪ ॲํ9 ݈߁ ࢦඪ: ෰ༀʹΑΔ ඃݧऀͷ ৺ཧతมԽ 「処⽅」には ⼼理的効果もある ʮϓϥηϘޮՌʯ 薬の処⽅X→健康指標Yの因果効果は？
174. 概念の”質的”吟味：媒介経路の観点から ▪ ॲํ9 ݈߁ ࢦඪ: ॲํʹΑΔ ؍࡯ऀͷ ৺ཧతมԽ ʮ؍࡯ऀޮՌʯ こちらの経路からの ⼼理的効果もありうる 薬の処⽅X→健康指標Yの因果効果は？
175. 概念の”質的”吟味：媒介経路の観点から ▪ ॲํ9 ݈߁ ࢦඪ: ༀͷ੒෼ͷ ༀཧ࡞༻ ෰ༀͷ༗ແ 薬の処⽅X→健康指標Yの因果効果は？ そもそも 「処⽅」しても 「服薬」してない ときもある ʮॲஔͷෆ९कʯ
176. 概念の”質的”吟味：媒介経路の観点から ▪ 統合: 処⽅X→健康指標Yの「より”厚い”記述」 ॲํ9 ݈߁ ࢦඪ: ॲํ9͋Γͳ͠ ΛϥϯμϜʹׂ෇
177. 概念の”質的”吟味：媒介経路の観点から ▪ ॲํ9 ݈߁ ࢦඪ: ༀͷ੒෼ͷ ༀཧ࡞༻ ෰ༀͷ༗ແ ෰ༀʹΑΔ ඃݧऀͷ ৺ཧతมԽ ॲํʹΑΔ ؍࡯ऀͷ ৺ཧతมԽ 統合: 処⽅X→健康指標Yの「より”厚い”記述」 ແ࡞ҝׂ෇͞Εͨॲํ9ͷޮՌ͸ɺશͯͷ ഔհܦ࿏Λ૯߹ͨ͠஋Ͱ͋Δʢ㱠ༀཧޮՌͰ͸ͳ͍ʂʣ
178. 概念の”質的”吟味：媒介経路の観点から ▪ ِༀPSຊༀ ͷॲํ9 ݈߁ ࢦඪ: ༀͷ੒෼ͷ ༀཧ࡞༻ ෰ༀͷ༗ແ ෰ༀʹΑΔ ඃݧऀͷ ৺ཧతมԽ ॲํʹΑΔ ؍࡯ऀͷ ৺ཧతมԽ 薬理作⽤の因果効果を知るには？
179. 概念の”質的”吟味：媒介経路の観点から ▪ ِༀPSຊༀ ͷॲํ9 ݈߁ ࢦඪ: ༀͷ੒෼ͷ ༀཧ࡞༻ ෰ༀͷ༗ແ ෰ༀʹΑΔ ඃݧऀͷ ৺ཧతมԽ ॲํʹΑΔ ؍࡯ऀͷ ৺ཧతมԽ 薬理作⽤の因果効果を知るには？ શһ෰ༀ ९कͷపఈ ϒϥΠϯυʹΑΔ੾அ ϒϥΠϯυ ʹΑΔ੾அ
180. 概念の”質的”吟味：媒介経路の観点から ▪ ِༀPSຊༀ ͷॲํ9 ݈߁ ࢦඪ: ༀͷ੒෼ͷ ༀཧ࡞༻ ෰ༀͷ༗ແ ෰ༀʹΑΔ ඃݧऀͷ ৺ཧతมԽ ॲํʹΑΔ ؍࡯ऀͷ ৺ཧతมԽ 薬理作⽤の因果効果を知るには？ શһ෰ༀ ९कͷపఈ ղܾ๏ɿೋॏ໡ݕ๏ʴ९कͷపఈ
181. 概念の”質的”吟味：媒介経路の観点から 処置（処⽅）の概念的吟味で何が変わったか？ ▪ 最初は薬Aの影響を⾒る際に、無作為化された「薬の処⽅の有無」が 「健康指標」に影響する単純な図式で捉えていた →「処⽅」が健康に与える影響が何に媒介されうるかを質的な知識に基 づき吟味した →「処⽅」の概念的意味がより⾼い解像度で理解された →「処⽅の因果効果」は「薬Aの⽣理的作⽤による因果効果」と同義でな いことが明⽰化された →「薬の⽣理的作⽤による効果」以外の経路を遮断するという実験デザ インの必要性が理解された →「本当に知りたい因果効果」にフォーカスするための具体的なデータ 取得プロトコル（⼆重盲検法）が検討された 「処理」概念と媒介経路の吟味により、因果効果の 推定・解釈がより適切に実施できるようになった （質的知⾒による概念の吟味は重要）
182. SUTVAを考える：要因間の依存性の観点から 効果の修飾：集団内での特性に応じて効果が異なる ▪ ༀ.ͷ ෰ༀ9 ೣͷ಺෼ൻ ࣬ױ,ͷ ࢦඪ: Φεೣʗ ϝεೣ 9ˠ:ͷฏۉҼՌޮՌ͸ ΄΅θϩͱਪఆ͞Εͨ ෰ༀ9ͷҼՌޮՌ͸ͳ͍ʁ ෰ༀ͸ແ࡞ҝʴೋॏ໡ݕ ˎԾ૝ͷྫͰ͢
183. SUTVAを考える：要因間の依存性の観点から 効果の修飾：集団内での特性に応じて効果が異なる ▪ ༀ.ͷ ෰ༀ9 ࣬ױ,ͷ ࢦඪ: Φεೣ ༀ.ͷ ෰ༀ9 ࣬ױ,ͷ ࢦඪ: ϝεೣ ૚ผʹղੳͯ͠ΈΔͱɾɾɾ  ʮ9ˠ:ͷฏۉҼՌޮՌʯ͸ʮ9ˠ:ͷҼՌޮՌʯͷ͋ΓΑ͏Λ ͏·͘൓ө͍ͯ͠ͳ͍ʢz9ͷҼՌޮՌzͷ֓೦త༳Εʣ ͦ΋ͦ΋Φεೣʹ͓͚Δʮ࣬ױ,ʯͱϝεೣʹ͓͚Δʮ࣬ױ,ʯ͸ ྨࣅͷපଶΛ΋ͭผͷ࣬ױ͔΋ʁʢҼՌਪ࿦ˠ࣭త֓೦΁ͷϑΟʔυόοΫʣ ˎԾ૝ͷྫͰ͢
184. SUTVAを考える：要因間の依存性の観点から 交互作⽤：処置後に決まる要因により効果が異なる ▪ ݈߁ ࢦඪ: ྍ๏$͋Γ ෰ༀ9 ʮ9ˠ:ͷҼՌޮՌʯ͸ྍ๏$ͷঢ়ଶʹ΋ґଘ͢Δ ʮ෰ༀ9ͷޮՌʯΛͲ͏ఆٛ͢Δ͔ʁ ʢͲͷੈքʹ͓͚ΔҼՌޮՌͳͷ͔ʹؔ͢Δෆఆੑ໰୊ͷ࠶೩ʣ ݈߁ ࢦඪ: ྍ๏$ͳ͠ ෰ༀ9 ˎԾ૝ͷྫͰ͢
185. SUTVAを考える：⽂脈依存性の観点から 処置の⽂脈依存性の存在 ▪ 服薬Xの因果効果 ＝𝑌'( )% だだ − 𝑌'( )% だだ これ、断薬の⽂脈でも同じ値になる？ 服薬あり 服薬なし ʮ9ˠ:ͷҼՌޮՌʯ͸ॲஔ9ͷจ຺ʹ΋ґଘ͢Δ 他の例だと例えば、同じ教育プログラムのRCTでも、⾃発的参加と、 謝⾦による参加と、強制参加の場合の処置効果は異なる場合など SUTVA条件により潜在結果は⼀義的な値を取ることが想定されているが、 0→1の効果と、1→0の効果は同じとは限らない
186. 因果推論と「質的な研究」｜ 概念吟味の重要性の歴史的な例（壊⾎病） 1747年 対照実験によりシトラスが船員の壊⾎病を防ぐ 効果が確認 シトラス→”酸味”→壊⾎病（の防⽌） Pearl and Mackenzie (2018)の記述に基づく 1800’s 英国海軍ではシトラス積載により壊⾎病は過去 のものに しかしmediatorの取り違えにより悲劇が・・・ シトラス→レモン→ライム→加熱濃縮ライムジュース へと代替されていった（代替のたびにビタミンが減っていく・・） 1900年前後の極地探検にて壊⾎病による死者が多発 本当のmediatorはビタミンC 柑橘の効果⾃体が疑問視され無視されるように・・ ▪
187. 因果推論と「質的な研究」｜ Pearl and Mackenzie (2018)の記述に基づく シトラス→”酸味”→壊⾎病（の防⽌） しかしmediatorの取り違えにより悲劇が・・・ シトラス→レモン→ライム→加熱濃縮ライムジュース へと代替されていった（代替のたびにビタミンが減っていく・・） 現実の系において”処置”の同⼀性の担保は必ずし も簡単な話ではない e.g., RCTが⾏われた”意識の⾼い”学校と同⼀の ”処置”を普通の学校で実現できるのか？ 異なる対象／⽂脈における“処置”の同⼀性の担保 は質的な理解度に依存する 概念吟味の重要性の歴史的な例（壊⾎病） ▪
188. II.2 概念の吟味の⼩まとめ Rubinの因果モデルはSUTVA条件を所与の前提とす ることで、理論的な定義と解析プロトコルの⾒通し を明確にしている ⼀⽅、現実の系では(hidden) multiple versionsが あることは普通である、あるいは、 概念の”厚い”記述への志向はSUTVA条件の成⽴の 難しさを明らかにしがち →この話は第IV部へつづく IV. さて、われわれは何を知りたいのか？ 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
189. 2. そのXは本当にXか：概念的妥当性を吟味する 3. その因果効果はどこまで”拡⼤解釈”できるか パートIIの⽬次 ・ その1：媒介経路の観点から ・ その2：要因間の依存性の観点から ・ その3：⽂脈依存性の観点から ・ その4：測定の観点から ・ ターゲット妥当性の分解 1. 可能世界への”チャネル”としてのSUTVA条件
190. 解釈の”質的”吟味：ターゲット妥当性の分解 「サンプル集団」と「ターゲット集団」の関係の類型 ▪ (A) サンプル元集団が ターゲット集団 (B) サンプル元集団の 中の⼀部の特性をもっ たサブ集団がターゲッ ト集団 (C) サンプル元集団と は異なる集団が ターゲット集団 サンプル元集団 ＝ターゲット集団 サンプル 集団 ターゲット妥当性＝サンプル集団の推定値をターゲット 集団に適応することの妥当性（内的／⼀般化／外的妥当性等を包含する概念） サンプル元集団 サンプル 集団 ターゲット集団 サンプル元集団 サンプル 集団 ターゲット 集団 ⼀橋⼤⽣ ⼀橋⼤⽣ ⼀橋⼤⽣中の 不登校リスク 集団 ⼀橋⼤⽣ ハーバード⼤⽣ ⼀橋⼤⽣ ハーバード⼤⽣
191. 解釈の”質的”吟味：ターゲット妥当性の分解 “⼀般化可能性”を巡るバイアスの分解と整理 ▪ サンプル元集団= ターゲット集団に おける真の因果効果 サンプル集団に おける真の因果効果 サンプル集団における 処置群間の差 BIAS於サンプリング BIAS於⼀般化 ＝BIAS於サンプリング +BIAS於交絡 サンプル元集団 ＝ターゲット集団 サンプル 集団 ⼀橋⼤⽣ ⼀橋⼤⽣ BIAS於交絡 共変量は バランシング しているか？ サンプリングは ランダムか？
192. 解釈の”質的”吟味：ターゲット妥当性の分解 “移設（外的）可能性”を巡るバイアスの分解と整理 ▪ サンプル元集団Gに おける真の因果効果 サンプル集団に おける真の因果効果 ターゲット集団Hに おける真の因果効果 BIAS於サンプリング@ 集団G BIAS於移設 ＝BIAS於潜在結果 +BIAS於サンプリング +BIAS於交絡 BIAS於交絡 @集団G サンプル集団における 処置群間の差 サンプル元集団 サンプル 集団 ターゲット 集団 ハーバード⼤⽣ ⼀橋⼤⽣ ハーバード⼤⽣ サンプリングは ランダムか？ 共変量は バランシング しているか？ BIAS於移設@潜在結果 潜在反応は代替可能か？ 「集団間（もしくは個体間）での 関数fおよび特性Zの分布の同等性」 構造的 類似性 特性的 類似性
193. 解釈の”質的”吟味：ターゲット妥当性の分解 “移設（外的）可能性”を巡るバイアスの分解と整理 ▪ サンプル元集団Gに おける真の因果効果 サンプル集団に おける真の因果効果 ターゲット集団Hに おける真の因果効果 BIAS於サンプリング@ 集団G BIAS於移設 ＝BIAS於潜在結果 +BIAS於サンプリング +BIAS於交絡 BIAS於交絡 @集団G サンプル集団における 処置群間の差 サンプル元集団 サンプル 集団 ターゲット 集団 ハーバード⼤⽣ ⼀橋⼤⽣ ハーバード⼤⽣ サンプリングは ランダムか？ 共変量は バランシング しているか？ BIAS於移設@潜在結果 潜在反応は代替可能か？ 「集団間（もしくは個体間）での 関数fおよび特性Zの分布の同等性」 構造的 類似性 特性的 類似性 ここのバイアスのチェックに 関して”量的”のアドバンテージ がどの程度あるかは実際問題と してかなり怪しいのではないか （c.f., RCT）
194. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性の例 Shoda et al. (1990) マシュマロ実験 「マシュマロを⾷べるのを我慢できた（7分待てた）⼦／ できなかった⼦」の⻑期追跡により、⼦供の”⾃制⼼ （delay of gratification）”が、その後の学業や社会⾏動にお ける成功において重要であることを⽰唆した 近年の追試的研究により、 - 概念解釈の妥当性 - 推定結果の移設可能性 の両⾯において結論（の捉え⽅）が疑問視されて きている
195. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性の例 Watts et al. (2018)マシュマロ実験の再試験 より⼤きくかつ多様なsample of childrenを調査 より洗練された統計モデルを使⽤ Shodaらで⽰されていた”delay of gratification”と後 のachievementの関係は社会経済的要因と認知能⼒ で調整すると⼤幅に減少し、統計的調整⽅法によ り有意になったりならなったりする程度となった 尚、Shodaらの研究はいわゆる「再現性が怪しい」という話とは少し違って、 「研究／実験そのものの不備」ではなく（ 解釈はともかく結果⾃体には⼀定の 再現性は認められる）、hypeがひどい（概念と移設可能性を真剣に吟味せずに ⼀般書で広めた点がひどい）例と思われる ←マシュマロ・テストとその後続研究については森⼝（2021）に現 在の知⾒がまとまっている
196. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性の例 再試験では”最初の20秒”でほぼ決まっていた 「測定しているのものは”何か”」「媒介しているも のは”何か” 」を質的に吟味しないと解釈を誤りうる The relationship between delay of gratification and academic achievement might be driven by the ability to generate useful metacognitive strategies that will influence self-regulation throughout one’s life. Shodaら (1990)での「7分待てること」の概念解釈 + 認知能⼒で調整すると効果量は⼤幅減 測定されていたのは”そもそも何だった“のか？ 後続研究からむしろ「他者への信頼度」に近い指標であることが⽰唆されている
197. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性の例 ターゲット妥当性の観点から Shodaらの研究ではサンプル元 集団は⽐較的均⼀性が⾼い Shodaらの研究ではサンプル元 集団は多様性ありの設計 ଌఆ͞Εͨ l੍ࣗ৺z ޙͷ BDIJFWFNFOU ଌఆ͞Εͨ l੍ࣗ৺z ޙͷ BDIJFWFNFOU ଞऀ΁ͷ ৴པ౓ 4&4 എܠཁҼͷ ෼ࢄ͕খ͍͞ 移設可能性はかなり不明 ⽐較的に移設可能性の⾒通しあり 作業的 概念図の例
198. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性を⽰す事例 マシュマロ実験の真の教訓は「われわれ研究者に は”エビデンス”の安易な⼀般化への⾃制⼼がしば しば著しく⽋けている」という事実である 概念と移設可能性の質的＆量 的な吟味の不⼗分さにより、 社会的・経済的要因の重要性 を後景化し、⾃⼰責任論を強 化しうる⽅向への解釈に陥っ た感は否めない マシュマロ・テストから何を学ぶべきか？
199. II.3 “拡⼤解釈可能性”の⼩まとめ ターゲット妥当性＝ ＋サンプリングの妥当性 ＋共変量のバランシングの妥当性 ＋構造的類似性×特性分布の類似性 構造＋特性的類似性が本質的な場合 には量的分析のアドバンテージはあるのか？ そもそも質的な知⾒がなければ類似性をみるための 量的アプローチ⾃体が成⽴しなそうではある
200. パートIIの話の全体まとめ ▪ ▪ ▪ Rubinの因果モデルはSUTVA条件を所与の前提と することで、理論的な定義と解析プロトコルの⾒ 通しを明確にしている ⼀⽅、現実の系では(hidden) multiple versions があることは普通である、あるいは、概念の”厚 い”記述への志向はSUTVA条件の成⽴の難しさを 明らかにしがち ⼿元のサンプルからのターゲット妥当性を考える 際に、構造＋特性的類似性が本質的な場合には量 的分析のアドバンテージはあるのか？ 概念やターゲット妥当性の吟味に⽋けた量的分析 はしばしば実害があるレベルで解釈を誤りうる
201. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
202. 1. EBM/EBPsの「往路」と「復路」 2. EBM/EBPsにおける量的/質的の主従を考える ：斎藤清⼆（2011,2012,2018）を⼿がかりに パートIIIの⽬次 Evidence-Based Medicine/ Evidence-Based Practices
203. 1. EBM/EBPsの「往路」と「復路」 2. EBM/EBPsにおける量的/質的の主従を考える ：斎藤清⼆（2011,2012,2018）を⼿がかりに パートIIIの⽬次 Evidence-Based Medicine/ Evidence-Based Practices
204. 因果推論と「質的な研究」｜ Evidence-Based Practicesの”往路”と”復路” 個人 集団A 集団A’ 集団Σ 集団因果効果の 推定量 法則的（タイプ的） 個別的（トークン的） “往路” “復路” 論⽂を書くことを⽣業とする「研究者」 という⼈種はしばしば”往路”までのこと しか真剣に考えていない 𝛼, 𝛽, 𝛾 … エビデンス の⽣産 個人 個人 エビデンス の利⽤ 𝑎, 𝑏, 𝑐 … 𝑠, 𝑡, 𝑢 … トークンto タイプ タイプto トークン 研究者は”丘の上”で論⽂を書く 復路では移設可能性／外的 妥当性が問題となる
205. 因果推論と「質的な研究」｜ 個人 集団A 集団A’ 集団Σ 集団因果効果の 推定量 “往路” “復路” 𝛼, 𝛽, 𝛾 … エビデンス の⽣産 個人 個人 エビデンス の利⽤ 𝑎, 𝑏, 𝑐 … 𝑠, 𝑡, 𝑢 … トークンto タイプ EBPsにおける「RCT最強論」の半可通性 元ネタのEBMedicineは⼀定の ⽣物学的⻫⼀性を想定しうる 点で、EBPMと⼤きく異なる 「往路最強」かもだが、その強みのモデルフリー性 ⾃体が”復路”での⼤きな脆弱性を⽣む ＊ Cf., Deaton and Cartwright 2017 ⼀般にRCTの内部には異なる状況 間の”距離”に対する情報が皆無 ? ? ? ? ?
206. 因果推論と「質的な研究」｜ 質的な”⼿触り”の重要性：介⼊の⽂脈的な意味 𝑌D (ばんとらいん) = 𝑓(𝑋, 観測可能＆数量化可能なもの, 数量化不可能なもの, 観測不可能なもの, 固有なもの, ⾮固有なもの … ) そもそもを⾔えば、個体の潜在反応は観測可能な もの／不可能なもの、数量化可能なもの／不可能 なもの／、固有なもの／⾮固有なものの全てに よって規定されている 観測可能＆数量化可能なもの以外の個体の諸々は ”往路”にて”誤差”として⽚付けられてしまう 本当の「ばんとらいん」 の全て 半可通にとっての「ばんとらいん」 ぼくの名前は ばんとらいん ばんとらいんにとって”ばんとらいんでないもの”は⼀つもないのに！
207. 因果推論と「質的な研究」｜ 𝑌D ばんとらいん = 𝑓(𝑋, 𝐶1, 𝐶2, 𝐶3, 𝐶4, 𝐶5, 𝐶6, 𝐶7, 𝐶8, 𝐶9, 𝐶10 … ) 「（やや／概ね／完全に）断⽚的なもの」 ばんと らいん X→Yの法則性を 考える上で 交絡要因として 重要な特性 （バックドア基準を 満たす変数セット） (X→)Yの中範囲 の⼀般性を 考える上で ⽂脈上 重要な特性 特定の ⽂脈に 回収され ない特性 (X→)Yの固有性 を考える上で ⽂脈上 重要な特性 集団の因果効果の 推定の観点からは “断⽚的なノイズ” 集団の因果効果の 推定のために 考慮が必要 ばんとらいんにおけるX→Yの全体像は ここの特性も全て関わる
208. 因果推論と「質的な研究」｜ 𝑌D ばんとらいん = 𝑓(𝑋, 𝐶1, 𝐶2, 𝐶3, 𝐶4, 𝐶5, 𝐶6, 𝐶7, 𝐶8, 𝐶9, 𝐶10 … ) 「（やや／概ね／完全に）断⽚的なもの」 ばんと らいん X→Yの法則性を 考える上で 交絡要因として 重要な特性 （バックドア基準を 満たす変数セット） X→Yの中範囲の ⼀般性を 考える上で ⽂脈上 重要な特性 特定の ⽂脈に 回収され ない特性 Y→Xの固有性を 考える上で ⽂脈上 重要な特性 RCTではX以外の 全ての要因を “誤差”とみなす （断⽚的なノイズの集積） ばんとらいんの 固有性はほぼ 消去される
209. 因果推論と「質的な研究」｜ 個人 集団A 集団A’ 集団Σ 集団因果効果の 推定量 “往路” “復路” 𝛼, 𝛽, 𝛾 … エビデンス の⽣産 個人 個人 エビデンス の利⽤ 𝑎, 𝑏, 𝑐 … 𝑠, 𝑡, 𝑢 … トークンto タイプ タイプto トークン 質的な”⼿触り”の重要性：介⼊の⽂脈的な意味 例えばリスク・コミュニケーションの 現場では「数値」のナラティブ内での 位置づけまでの考慮が必要となる 往路で”誤差”として⽚付けた「断⽚的なもの」に ”復路”でふたたび向き合うべきときがある こういう感覚が分からない⼈をコミュニケーター役にすると⼤抵ろくなことにならない
210. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 ▪ 臨床⼼理学/内科学/⼼⾝医学の⽴場からエビデンスと ナラティブの関係について⻑年論考が重ねられている 斎藤(2011) ＊尚、林は⾃分のEBPM/リスクコミュケーションへの興味と重ねあわせて読んでます 斎藤(2012=2016) 斎藤(2018) 会話形式で⼀⾒カジュアルだが、 実は⼀番踏み込んだ議論をしてるかも EBM+NBM論を⽐較的フォーマル に議論した内容 EBM+NBM、量的×質的論の 現状での到達点
211. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 EBMの概念解釈の3類型（斎藤2012/松下2021より引⽤） ▪ 松下(2021)『⼼理療法における「エビデンス」とは何か』より引⽤ 1) EBM 正統派(批判的吟味派) EBM の定義と⽅法論を遵守する⽴場。エビデンスを疫学研究に よる経験的実証的情報と定義 し、EBM の本質をエビデンスおよ びエビデンスを⽤いた臨床実践の批判的吟味におく。 2) ガイドライン派(グローバル・スタンダード派) エビデンスの定義は 1)と共通するが、EBM を「エビデンスに基 づいた診療ガイドラインを作成し、それを医療実践に普及させる こと」と捉える⽴場。エビデンスの格付けを重視する。 3) 伝統科学派(病態⽣理派) エビデンスを臨床疫学的な情報に限定せず、⽣物化学的な理論や 病態⽣理を推定する実験的 研究の成果を重視しようとする⽴場。 臨床実践が出発点 Treatment/治療法が出発点 本来のEBMは1)であるが、残念ながら2)として喧伝されている場合が多い
212. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 Narrative-Based Medicineとは（斎藤2012=2016より） ▪ 『患者が主観的に体験する物語』を全⾯的に尊重し、医療 者と患者との対話を通じて、新しい物語を共同構成してい くことを重視する医療 「病いを、患者の⼈⽣という⼤きな物語の中で展開する⼀ つの『物語』であるとみなし、患者を『物語を語る主体』 として尊重する⼀⽅で、医学的な疾患概念や治療法もあく まで⼀つの『医療側の物語』と捉え、さらに治療とは両者 の物語をすり合わせる中から『新たな物語』を作り出して いくプロセスである、と考えるような医療」 斎藤(2016) より引⽤ （主張者の中でもグラデーションはあるが）基本的に解釈主義的な⽴場
213. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 ⼆重構造をもつ実践科学としてのEBMの構造+NBM ▪ エビデンスの⽣成プロセス 厳密に統計的な実験デザイン 論⽂化 エビデンスの評価プロセス ⽅法論的妥当性の吟味 臨床的妥当性の吟味 エビ デンス EBMの実践プロセス 患者の問題の定式化 問題についての情報収集 得られた情報の批判的吟味 得られた情報の患者への適⽤ これまでの実践の評価 NBMによる 対話的臨床実践 ⽣成 利⽤ 吟味 吟味 合法則的合理性 合⽬的的合理性 合⽬的的合理性 実証主義的 世界観 解釈主義的世界観 客観的に存在する疾患 を治療/除去する≈ 患者の抱える「病いの物語」を 編み直す ＊この図は⻫藤(2018)を元に林が勝⼿にNBM部分を増築
214. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 EBMとNBMの調和/包含関係を考える ▪ 斎藤（2011）で⽰された３つの考え⽅の現状分類 ① EBMとNBMは相互に補完的であり、NBMを加えることに よってEBMの体系は完成するという楽観的な考え⽅ ② EBMとNBMは異なる２つの世界観であるが、患者と医師 の出会いの場⾯において共存しうるとする慎重な考え⽅ ③ EBMとNBMは異なる２つの世界観であるが、患者と医師 の対話の現場において、NBMはEBMを包摂／統合するとい う⼤胆な考え⽅ →⽇本の多くのEBMの専⾨家はこの考え⽅。臨床実践の実感に合う →NBMの主張者の⼤部分はこの考え⽅（認識論の違いに対して慎重） →斎藤(2011)の⽴場。EBMはそれ⾃体が「有効に活⽤すべき⼀つの物語」 プラグマティズム的折衷／統合ともいえるかも 斎藤(2011)は「ナラティブ三年エビ⼋年」、つまり、EBMをしっかりと勉強して理解した⼈はNBMをすぐ理解するとも書いている 「偽権威」や「トンデモ医療」から患者を護るためにはEBMのしっかりした理解は必須
215. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 （補⾜）権威主義との闘争としての初期EBM c.f., https://togetter.com/li/1746410 image
216. 斎藤清⼆（2011,2012,2018）のEBM/NBM論 （補⾜）権威主義との闘争としての初期EBM c.f., h<ps://toge<er.com/li/1746410 近藤誠 外科医たちの 権威主義的な 癌治療実践 権威を笠にきた 苛烈な⼈格攻撃 闇落ち ＊「ナラティブ三年エビ⼋年」を理解する上での背景のひとつ 「海外の論⽂がどうした、今ま でこうしてたんだからいいんだ、 俺の勘ではこうだからごちゃご ちゃいうな」 近藤⽒が最初期に⾔っていたこと は今では「当たり前のこと」に 海外の論⽂や エビデンスを もとにした批判
217. 量的が”主”で、質的が”従”なのか？ ▪ 個別の”datum”によって”更新”されることで初めて完成 されるのがEBPsであれば、量的が”従”で質的が“主”と も⾔えるのかも？ 階層ベイズ的なアナロジーで EBPsを“個別→法則→個別”の往復運動であるとする ならば、法則的エビデンスは”事前分布（個体の データが無い場合のデフォルト情報）”を提供する ものとしてみなせるようにも⾒える 施策の有効性に関する 法則的エビデンス 個別の対話や （やや／概ね／完全に） ”断⽚的”な情報 個別の施策・ 意思決定 重み付け・ 更新 前菜 主菜 “事前分布” 個別の”datum”
218. III.のまとめ EBPsはしばしば「個別→法則→個別」の往復運動 最後の「個別」こそ最終的な⽬的であるならば量的 分析は”事前分布”を提供し、それを質的知⾒が”更 新”する過程である、という⾒⽅もあり？ EBMとNBMの折衷／統合は、プラグマティズム （合⽬的性志向）の半ば必然的帰結であるやも？ “ナラティブ3年エビ8年”→(初期)EBMは既存の権 威主義との闘争でもあったが、質的研究は「権威 主義」や「トンデモ」とどう戦いうるのか？ ▪ ▪ ▪ ▪ 良くも悪くも英語論⽂IF⾄上主義とエビデンスヒエラルキー的有効性評価体系が昭和的医学 権威主義を打破した側⾯はある（それがまたある種のエビデンス／グローバル権威主義へ と繋がるのはまた別の話として）。
219. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
220. 統計学因果推論と社会学：筒井(2019)を元に 筒井(2019）『計量社会学と因果推論』 ▪ ཧ࿦ͱํ๏ 34ר1߸ 2019೥
     ਺ཧࣾձֶձ Sociological Theory and Methods Vol.34 No.1 2019 Japanese Association for Mathematical Sociology ಛू࿦จ ܭྔࣾձֶͱҼՌਪ࿦ɿ ؍࡯σʔλʹج͍ͮͨࣾձͷཧղʹ޲͚ͯ ౵Ҫ ३໵ʢ໋ཱؗେֶʣ ʦཁ໿ʧ ౷ܭతҼՌਪ࿦͸ܭྔ෼ੳͷओྲྀͱͳ͍ͬͯΔ͕ɼܭྔࣾձֶʹ͓͚Δͦͷҙຯ΍Πϯ ύΫτʹ͍ͭͯମܥతʹ࿦ͨ͡ݚڀ͸͍·ͩʹগͳ͍ɽຊ࿦จͰ͸ɼҼՌਪ࿦ϞσϧΛؚ Ήܭྔ෼ੳͷख๏ʹ͍ͭͯɼҟ࣭ੑͱ͍͏֓೦Λ࣠ʹ੔ཧ͠ɼͦͷ্Ͱܭྔࣾձֶ͕ҟ࣭ ੑʹରͯ͠ଞͷ෼໺ͷख๏ͱ͸ҟͳͬͨΞϓϩʔνΛͱΔ܏޲͕͋Δ͜ͱΛࣔ͢ɽ͜ͷ͜ ͱ͸ɼϚϧνϨϕϧ෼ੳͱ΋ݺ͹ΕΔࠞ߹ޮՌϞσϧͷ׆༻ʹ͓͍ͯ໌Β͔Ͱ͋Δɽ͞Β ʹɼհೖ΍੾அΛ༻͍ΔҼՌਪ࿦ΞϓϩʔνͱɼཁҼؒͷؔ࿈ੑΛڧௐ͢Δܭྔࣾձֶత Ξϓϩʔνͷҧ͍Λઆ໌͠ɼͦΕ͕ਓʑͷ֓೦࿈ؔΛࢀর͢ΔࣾձֶͷಛੑͷݱΕͰ͋Δɼ ͱ͍͏͜ͱΛ࿦͡Δɽ
221. 統計学因果推論と社会学：筒井(2019)を元に 個体の異質性を巡る観点の違いの議論（筒井2019） ▪ ただ，ここではもうひとつのインプリケーションの⽅がより重要であると思わ れる.それは，計量社会学がいわゆるデモグラフィック要因(個体の異質性)に強 い関⼼を払うがゆえに，個体の異質性の効果を除去することを⽬的とするルー ビン・モデルやその他の因果推論志向が強い⼿法を忌避したか，あるいは導⼊ の必要性が強く意識されなかった，という可能性である(Xie 2007). 社会学で統計的因果推論⼿法の導⼊が起きていないことに関して: 筒井(2019)より引⽤)
222. 統計学因果推論と社会学：筒井(2019)を元に 混合効果モデルを軸に異質性への態度を議論 ▪ 以上のように，混合効果モデルの発想は，個体の異質性の効果を除去した上で 固定効果 の推定値を得るという⽅針とは逆に，固定効果の影響を除去した上で 個体特有の効果の⼤ きさを追究する，という点にある. 因果推論においては異 質性がノイズであるのに対して， 混合効果モデルにおいては固定効果こそがノ イズなのである. 混合効果モデルは y = X β + Zu + ε という基本形で表現できる.ここにおいて，Z は 説明 変数であり，たいてい X(切⽚を含む)の部分集合である.u は個体(クラス ター)間の 変量効果ベクトルであり，個体(クラスター)数× Z の⾏列である.ここ で個体とは，デー タの最⼩単位が集合して作られる集まり(クラスター)を指す. ⽣徒個⼈が最⼩単位であ るとき，その集まりである学級などが想定される. 混合効果モデルと固定効果モデルの⽐較の議論: 筒井(2019)より引⽤)
223. 因果推論と「質的な研究」｜ 𝑌D ばんとらいん = 𝑓(𝑋, 𝐶1, 𝐶2, 𝐶3, 𝐶4, 𝐶5, 𝐶6, 𝐶7, 𝐶8, 𝐶9, 𝐶10 … ) 「（やや／概ね／完全に）断⽚的なもの」 ばんと らいん X→Yの法則性を 考える上で 交絡要因として 重要な特性 （バックドア基準を 満たす変数セット） (X→)Yの中範囲 の⼀般性を 考える上で ⽂脈上 重要な特性 特定の ⽂脈に 回収され ない特性 (X→)Yの固有性 を考える上で ⽂脈上 重要な特性 集団の因果効果の 推定の観点からは “断⽚的なノイズ” 集団の因果効果の 推定のために 考慮が必要 ばんとらいんにおけるX→Yの全体像は ここの特性も全て関わる 個体の異質性 処置X （再掲スライド） 固定効果
224. 統計学因果推論と社会学：筒井(2019)を元に ルービンの枠組みと「概念の厚み」 ▪ デモグラフィック変数が，少なくとも単純なかたちでは因果推論における「措置/ 介 ⼊」として考えにくくなるのは，さしあたりはそれが「動かない/動きにくい」 からだと 考えることができる.ただし，因果推論における介⼊の特徴は，それだけ だと捉えられな い，ということを強調しておきたい.というのは，しばしば注⽬さ れる⾃⼰選択との違い が，個体内変動がないことだけでは説明できないからだ.よ り根本的な違いは，デモグラ フィック変数が通常，介⼊変数と違って，《意味的 な厚みのある概念を反映》した変数で あるということにある.たとえば「学歴」に ついて考えてみよう.学歴は，家庭的背景， 能⼒，シグナル機能，アスピレーショ ンなど，さまざまな概念との関連において有意味な要因として理解される.関連す る要素(「他の条件」)をすべて剥ぎ取ってしまうと，もは や私たちは学歴という概 念について理解することが難しくなる.というのは，パールの語彙を使うと，ある 選択を介⼊として捉えるということは，因果ダイアグラムにおける「切 除」を意 味しているからだ. 因果推論の枠内においても，いわゆるメカニズムの分析におい ては，措置の効果の有意 味な解釈が可能である.たとえば学歴を割当することは原 理的には可能であるが，その結 果の解釈においては，学歴と関連する概念(たとえ ば教育を通じて取得された能⼒やスキ ル)を考慮した上で，効果を及ぼすメカニズ ムの分析・解釈がなされることになる.他⽅で， ルービン・モデルを使⽤した場 合，焦点を当てる要因が社会的なプロセスにおいてどのよ うな重要性を持ってい るのかについては，全体的な推論や解釈を⾏うことはできなくなる. 筒井(2019)より引⽤)
225. 統計学因果推論と社会学：筒井(2019)を元に 社会学は何が知りたいのか？ ▪ 例外は，シエや⽯⽥が記録している，ダンカンが「遺産」として残したとされる ⼀ 連の知⾒(Xie 2007; Goodman 2007)や，最近ではゴールドソープによる精⼒的な 発信で あろう.ダンカンは異質性こそが社会学の研究対象であると主張し(⽯⽥ 2012: 11)，異質性をノイズとして捉える因果推論とは距離を置いた.ゴールドソー プは，(計量)社会 学を population science として特徴づける.population science とは， 統計データを通じて⼈ ⼝集団の規則性を確率論的に把握するものである 12) .ゴール ドソープによれば，統計学 的に発⾒された集団の規則性(たとえば出⽣率の低下)を もたらした諸原因について考 察する際に，ルービンの potential outcomes アプロー チの活動度が下がる(Goldthorpe 2016: 107). 筒井(2019)より引⽤ 観察データにおける変数は，それぞれ⼀定の「概念の厚み」を持っている.これと 対局 なのが実験における介⼊で，介⼊は⽂脈から切断された操作として，概念的 には空洞で あっても実施し，結果を観察することができる.もちろん，観察のレベ ルにおいても外⽣ 性のショックや制度的な断絶におけるように，⽂脈からの切断 は⽣じる.むしろ，だから こそ因果推論ではこれらのショックや切断を⾃然実験と して，介⼊の代わりに⽤いるので ある.他⽅で，社会的なプロセスの記述や理解 は，通常はこういった断絶を含まない，要因間の関連を想定してなされる.その意 味で，社会的プロセスの記述が社会学的な課題で あるとすれば，それは観察デー タに基づいて⾏われるしかない.
226. 統計学因果推論と社会学：筒井(2019)を元に 研究課題に対する依存性 ▪ ただし，要素間の連関と断絶とは，そのあいだにグラデーションを含んだもので あり， ある要因の概念をどこまで分節化するのかは，当該概念の⼀般的理解のみ ならず，研究課 題に依存する.たとえば「学歴」と関連する要素(家庭的背景，ア スピレーション，能⼒， シグナリング等々)を切り離していき，学歴の⼀般的な概 念理解を掘り崩すに⾄ったと しても，介⼊プログラムとしてそれが有意味である 場合には，問題がない.あるいは観察 データの分析において，出⾝家庭要因を切り 離した上で学歴による違いをみるといった場 合，進学が⾃⼰選択であることを勘 案しても，出⾝家庭要因とは独⽴した学歴の影響の⼤ きさをある程度推察するこ とができる.肝⼼なのは，因果の捉え⽅の⽔準にそった分析を ⾏うことである. 筒井(2019)より引⽤
227. ҼՌޮՌਪఆͷجຊతͳߟ͑ํʢ౵Ҫ͔Βͷ຋Ҋʣ 9 : ; 9 : ; ここの関係だけ⾒たい！ （それ以外はキャンセルしたい） ランダム化 均⼀化 統計学的調整 学歴 （⼤学） 親の 社会階層 収⼊ 『計量社会学と因果推論』 統計学因果推論と社会学：筒井(2019)を元に
228. ࣾձֶͰͷߟ͑ํʢ౵Ҫ͔Βͷ຋Ҋʣ 9 : ; 9 : ; 社会学は ここまでスコープに ⼊れてナンボでしょ Zの存在を キャンセルするのは ないわー 収⼊ 親の 社会階層 学歴 （⼤学） →より”説明”志向の⼿法選択 『計量社会学と因果推論』 統計学因果推論と社会学：筒井(2019)を元に 安易に“キャンセル”してはいけないものがあるのでは
229. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性の例 ターゲット妥当性の観点から Shodaらの研究ではサンプル元 集団は⽐較的均⼀性が⾼い Shodaらの研究ではサンプル元 集団は多様性ありの設計 ଌఆ͞Εͨ l੍ࣗ৺z ޙͷ BDIJFWFNFOU ଌఆ͞Εͨ l੍ࣗ৺z ޙͷ BDIJFWFNFOU ଞऀ΁ͷ ৴པ౓ 4&4 എܠཁҼͷ ෼ࢄ͕খ͍͞ 移設可能性はかなり不明 ⽐較的に移設可能性の⾒通しあり 作業的 概念図の例 （再掲スライド）
230. 因果推論と「質的な研究」｜ 概念と移設可能性の吟味の重要性を⽰す事例 マシュマロ実験の真の教訓は「われわれ研究者に は”エビデンス”の安易な⼀般化への⾃制⼼がしば しば著しく⽋けている」という事実である 概念と移設可能性の質的＆量 的な吟味の不⼗分さにより、 社会的・経済的要因の重要性 を後景化し、⾃⼰責任論を強 化しうる⽅向への解釈に陥っ た感は否めない マシュマロ・テストから何を学ぶべきか？ ֶज़ํ๏࿦తzΩϟϯηϧzͷةݥੑ （再掲スライド）
231. 因果推論のはじまりの場所｜ 初⼼に戻り「対象のありよう」を丁寧に考える 10個のリンゴがあります 統計の教科書が想定してる「ありよう」は Exchangeable! （再掲スライド）
232. 因果推論のはじまりの場所｜ 因果推論でまず重要なのは「対象集団のありよ う」について丁寧に省察することである ある意味、統計的因果推論とは集団内の異質 性／⾮均⼀性を巡る体系である（後述） 「ありよう」の質的な省察は推論の⼟台をなす Exchangeable? （再掲スライド）
233. 因果推論のはじまりの場所｜ 因果推論でまず重要なのは「対象集団のありよ う」について丁寧に省察することである ある意味、統計的因果推論とは集団内の異質 性／⾮均⼀性を巡る体系である（後述） 「ありよう」の質的な省察は推論の⼟台をなす Exchangeable? （再掲スライド）
234. 導⼊｜統計的推測の主要な3⽬的（再掲） ⽬的の 分類 やりたいことの内容 事例のイメージ （⼤学の留年の例で） 理念的な理想をなす 解析法のイメージ しばしば 現実的な 解析法 予測 現在得られている情報 から、将来／未知の値 ／状況を予測したい （要因への介⼊は伴わない場合） 多数の変数の情報から、留年 リスクの⾼い学⽣を⾒つけた い（→その後⽀援につなげる など） 機械学習 回帰分析！ 説明 要素間の関連性を説明 したい なぜ学⽣が留年するかを説明 したい（要因を探索し理解／ 対策につなげたい） フルスタックの 構造⽅程式 （含む因⼦分析等） 回帰分析！ 因果／介⼊ 効果の推定 ある特定の要因に介⼊ したときの、ある特定 の要因への効果を知り たい 特定のプログラムの導⼊が留 年を防⽌する効果を知りたい ランダム化⽐較試験 回帰分析！ （異論は認める） ＊ ＊互いに排反な分類という意味ではないので注意 パート1では「因果効果の推定」が主テーマ ▪ 構造的因果モデルはこれらの⽅向性の 異なる世界線を架橋しうる点が嬉しい ʮհೖޮՌʯͱʮઆ໌ʯ͕੾Εͯ ͠·͏͜ͱͷةݥੑ （再掲スライド）
235. 「説明モデルの良さ」の 基準は明確ではないところが どうしても残る 実態としては、分野内で蓄積されてきた既存知⾒との ある種のcoherencyにより評価されているところもあり ありようの「説明／分析」につきまとう悩み ▪ ⼀⽅、Rubinの因果モデルはSUTVA条件を所与の 前提とすることで、理論的な定義と解析プロトコ ルの⾒通しを明確にしている 査読共同体の中で、何をしたら何が⽰されたこと になるかのコンセンサスが明確である （そうした中で研究を蓄積していくことのメリットはそれ はそれとして⾊々とあるだろう） 統計学因果推論と社会学：筒井(2019)を元に
236. IV.のまとめ 「ありよう」の質的な省察は推論の⼟台をなす 背景要因を キャンセルアウト 背景要因との関連や “（やや／概ね／完全に） 断⽚的な もの”⾃体が興味の対象 尚、近年はRubin系の因果推論でも異質性への関⼼が⾼まっているが、主にマーケティング的な関⼼であり 固有性に対する興味／リスペクトなどは基本的に無い さて、われわれは何を知りたいのか？
237. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
238. I. 因果推論、構造から⾒るか？差分から⾒るか？ II. 統計的因果推論における概念/解釈妥当性の 検討、そのドメイン知識との繋がり III. 法則性と固有性：EBM/EBPsにおける 「（やや／概ね／完全に）断⽚的なもの」 本⽇の構成（パートI-IV） IV. さて、われわれは何を知りたいのか？ 〜⾮巡回有向グラフ(DAG)、潜在反応モデル、そして構造的因果 モデル、加えて反事実シミュレーションについて 〜ドメイン知識で「SUTVAの井⼾の底」を掘る 〜斎藤清⼆(2011, 2012, 2018)を⼿がかりに量的/質的を考える 〜筒井淳也(2019)『計量社会学と因果推論』を⼿がかりに 「交絡要因の”キャンセルアウト”」の学問的含意を考える
239. 統計的因果推論のフレームワークとは いかなるものかを理解することで 統計的因果推論がいかに質的な知⾒を （しばしば暗黙に）その基盤としているか 統計的因果推論のスコープがしばしば いかに”狭い＆薄い”のか を議論したい 本⽇のねらい （再掲スライド）
240. • 統計的因果推論では陰に陽に対象についての質的な知⾒を その基盤としており、その概念的および構造的仮定の質的 な吟味は重要である 全体のまとめ・雑感 • 潜在結果モデルは「潜在結果の差分」として因果効果を定 式化し、因果効果推定の解析プロトコルへの明瞭な⾒通し を与える • 構造的因果モデルはissueの総合的理解へ向けた統計的／理 論的／理念的な枠組みの基盤となりうる • EBPsはしばしば「固有→法則→固有」の往復運動であり、 合⽬的的合理性の追求は半ば必然的に量的・質的の統合を 指向するのでは？ • 統計的因果推論における⽅法論的”キャンセル”の学術的含 意と（その裏⾯の「説明」の学術的捉えがたさ）をどう考えるか