Mediation for the 21st Century

Similar documents
Revision list for Pearl s THE FOUNDATIONS OF CAUSAL INFERENCE

Causality II: How does causal inference fit into public health and what it is the role of statistics?

Causal mediation analysis: Definition of effects and common identification assumptions

Binary Logistic Regression

STA441: Spring Multiple Regression. This slide show is a free open source document. See the last slide for copyright information.

Marginal versus conditional effects: does it make a difference? Mireille Schnitzer, PhD Université de Montréal

Chapter 6. Logistic Regression. 6.1 A linear model for the log odds

Statistics in medicine

Statistics Boot Camp. Dr. Stephanie Lane Institute for Defense Analyses DATAWorks 2018

y ˆ i = ˆ " T u i ( i th fitted value or i th fit)

Statistical Analysis of Causal Mechanisms

Business Statistics. Lecture 10: Correlation and Linear Regression

Logistic regression: Miscellaneous topics

Module 03 Lecture 14 Inferential Statistics ANOVA and TOI

Mediation Analysis for Health Disparities Research

Hypothesis testing, part 2. With some material from Howard Seltman, Blase Ur, Bilge Mutlu, Vibha Sazawal

Ignoring the matching variables in cohort studies - when is it valid, and why?

Mediation analyses. Advanced Psychometrics Methods in Cognitive Aging Research Workshop. June 6, 2016

An Introduction to Causal Mediation Analysis. Xu Qin University of Chicago Presented at the Central Iowa R User Group Meetup Aug 10, 2016

CAUSAL INFERENCE IN THE EMPIRICAL SCIENCES. Judea Pearl University of California Los Angeles (

Poisson regression: Further topics

AP Statistics. Chapter 6 Scatterplots, Association, and Correlation

Investigating mediation when counterfactuals are not metaphysical: Does sunlight exposure mediate the effect of eye-glasses on cataracts?

Hypothesis testing. Data to decisions

Empirical Bayes Moderation of Asymptotically Linear Parameters

Introduction to Statistical Analysis

Two-sample Categorical data: Testing

An introduction to biostatistics: part 1

OUTLINE CAUSAL INFERENCE: LOGICAL FOUNDATION AND NEW RESULTS. Judea Pearl University of California Los Angeles (

Casual Mediation Analysis

One-sample categorical data: approximate inference

Lecture 10: F -Tests, ANOVA and R 2

Introduction to Statistical Inference

Chapter 1 Statistical Inference

PS2.1 & 2.2: Linear Correlations PS2: Bivariate Statistics

STA441: Spring Multiple Regression. More than one explanatory variable at the same time

LECTURE 5. Introduction to Econometrics. Hypothesis testing

appstats27.notebook April 06, 2017

Chapter 19: Logistic regression

Flexible mediation analysis in the presence of non-linear relations: beyond the mediation formula.

SC705: Advanced Statistics Instructor: Natasha Sarkisian Class notes: Introduction to Structural Equation Modeling (SEM)

Do students sleep the recommended 8 hours a night on average?

1 Least Squares Estimation - multiple regression.

Lecture 4: Training a Classifier

f rot (Hz) L x (max)(erg s 1 )

An Introduction to Mplus and Path Analysis

Chapter 27 Summary Inferences for Regression

SPECIAL TOPICS IN REGRESSION ANALYSIS

Overview. Overview. Overview. Specific Examples. General Examples. Bivariate Regression & Correlation

OUTLINE THE MATHEMATICS OF CAUSAL INFERENCE IN STATISTICS. Judea Pearl University of California Los Angeles (

Computational Systems Biology: Biology X

Linear Regression. In this lecture we will study a particular type of regression model: the linear regression model

Comments on The Role of Large Scale Assessments in Research on Educational Effectiveness and School Development by Eckhard Klieme, Ph.D.

Logistic regression: Why we often can do what we think we can do. Maarten Buis 19 th UK Stata Users Group meeting, 10 Sept. 2015

Hierarchical Generalized Linear Models. ERSH 8990 REMS Seminar on HLM Last Lecture!

Causal Mechanisms Short Course Part II:

Stat 642, Lecture notes for 04/12/05 96

Final Exam. Name: Solution:

CS 147: Computer Systems Performance Analysis

General Regression Model

Statistical Methods for Causal Mediation Analysis

Lecture 14: Introduction to Poisson Regression

Modelling counts. Lecture 14: Introduction to Poisson Regression. Overview

Statistics: revision

1 A simple example. A short introduction to Bayesian statistics, part I Math 217 Probability and Statistics Prof. D.

Approximate analysis of covariance in trials in rare diseases, in particular rare cancers

Review of Statistics 101

Bounding the Probability of Causation in Mediation Analysis

where Female = 0 for males, = 1 for females Age is measured in years (22, 23, ) GPA is measured in units on a four-point scale (0, 1.22, 3.45, etc.

Earthquakes in Ohio? Teacher Directions and Lesson

Bag RED ORANGE GREEN YELLOW PURPLE Candies per Bag

Dynamics in Social Networks and Causality

DETERMINE whether the conditions for performing inference are met. CONSTRUCT and INTERPRET a confidence interval to compare two proportions.

Path Analysis. PRE 906: Structural Equation Modeling Lecture #5 February 18, PRE 906, SEM: Lecture 5 - Path Analysis

Lecture 9: Linear Regression

Correlation and regression

Lecture 4: Testing Stuff

regression analysis is a type of inferential statistics which tells us whether relationships between two or more variables exist

Difference Between Pair Differences v. 2 Samples

An Introduction to Path Analysis

Review of probability and statistics 1 / 31

STA6938-Logistic Regression Model

Introduction. Consider a variable X that is assumed to affect another variable Y. The variable X is called the causal variable and the

Linear Models for Regression CS534

13.1 Causal effects with continuous mediator and. predictors in their equations. The definitions for the direct, total indirect,

Variance. Standard deviation VAR = = value. Unbiased SD = SD = 10/23/2011. Functional Connectivity Correlation and Regression.

THE PEARSON CORRELATION COEFFICIENT

Basic Medical Statistics Course

2. Two binary operations (addition, denoted + and multiplication, denoted

REVIEW 8/2/2017 陈芳华东师大英语系

Lecture Discussion. Confounding, Non-Collapsibility, Precision, and Power Statistics Statistical Methods II. Presented February 27, 2018

Statistical Analysis of List Experiments

Math101, Sections 2 and 3, Spring 2008 Review Sheet for Exam #2:

HST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007

Confounding, mediation and colliding

Harvard University. Rigorous Research in Engineering Education

Course Introduction and Overview Descriptive Statistics Conceptualizations of Variance Review of the General Linear Model

Unpacking the Black-Box: Learning about Causal Mechanisms from Experimental and Observational Studies

Correlation and regression

Counterfactual Reasoning in Algorithmic Fairness

Transcription:

Mediation for the 21st Century Ross Boylan ross@biostat.ucsf.edu Center for Aids Prevention Studies and Division of Biostatistics University of California, San Francisco Mediation for the 21st Century p. 1

Fixed formula and page 27 since the presentation. 1-1

Motivating Problem Ethnic Group HIV+ Jail AA 51% 41% API 14% 13% Latino 43% 35% Logistic regression indicates that Jail raises the odds ratio of being HIV+ by 50%. How much of the group differences in HIV+ rates are accounted for by differences in incarceration? Binary outcome, 3-category predictor, binary mediator, categorical and continuous covariates. Mediation for the 21st Century p. 2

Overview Classical Mediation and Moderation Some Updates Limitations Structural Causal Models Causal Modeling Approach to Mediation Mediation for the 21st Century p. 3

Scope This not a talk about determining true causal relationships or estimating them from data. It is about how to interpret models once they are in hand. Mediation for the 21st Century p. 4

Classical Mediation Started with Baron & Kenny 1986: In general, a given variable may be said to function as a mediator to the extent that it accounts for the relation between the predictor and the criterion. Mediation for the 21st Century p. 5

4 Step Test For Mediation X c.. Y X c.... Y a b.. M 1. c is significant. Mediation for the 21st Century p. 6

4 Step Test For Mediation X c.. Y X c.... Y a b.. M 1. c is significant. 2. a is significant. Mediation for the 21st Century p. 6

4 Step Test For Mediation X c.. Y X c.... Y a b.. M 1. c is significant. 2. a is significant. 3. b is significant. Mediation for the 21st Century p. 6

4 Step Test For Mediation X c.. Y X c.... Y a b.. M 1. c is significant. 2. a is significant. 3. b is significant. 4. c < c. If c is not significant, total mediation. Mediation for the 21st Century p. 6

Modern Reconsiderations X c.. Y X c.... Y a.. M b Mediation for the 21st Century p. 7

Modern Reconsiderations X c.. Y X c.... Y a.. M b 1. ab is significant. Mediation for the 21st Century p. 7

Modern Reconsiderations X c.. Y X c.... Y a.. M b 1. ab is significant. That s it! Mediation for the 21st Century p. 7

Modern Reconsiderations X c.. Y X c.... Y a.. M b 1. ab is significant. That s it! c doesn t matter because of inconsistent mediation, e.g. c and c have opposite signs. Mediation for the 21st Century p. 7

Modern Reconsiderations X c.. Y X c.... Y a.. M b 1. ab is significant. That s it! c doesn t matter because of inconsistent mediation, e.g. c and c have opposite signs. Focus on ab rather than separate tests on each. Mediation for the 21st Century p. 7

Moderation In general terms, a moderator is a qualitative (e.g., sex, race, class) or quantitative (e.g., level of reward) variable that affects the direction and/or strength of the relation between an independent or predictor variable and a dependent or criterion variable. Baron and Kenny (1986) X c.. Y.... moderates..... M Mediation for the 21st Century p. 8

Limitations Big mess if M is a mediator and a moderator Little guidance in handling categorical variables and non-linear models. Abuse of the null hypothesis. Mediation for the 21st Century p. 9

The Null in Practice If any of a, b or c is not statistically significant, conclude there is no mediation. In more modern forms, if ab is not statistically significant, conclude there is no mediation. If the interaction term d in Y = c X + bm + dmx + ǫ is not statistically significant, conclude there is no moderation. Mediation for the 21st Century p. 10

Some Smaller Problems The relation might not be linear, and you probably didn t check. Even the tests of ab require huge sample sizes to have a reasonable chance of rejecting the null when the effects are small. Mediation for the 21st Century p. 11

The Bigger Problem Absence of evidence is not evidence of absence. Don Rumsfeld. Mediation for the 21st Century p. 12

The Bigger Problem Absence of evidence is not evidence of absence. Don Rumsfeld. The fact that the data are consistent with null does not mean the null is true. Mediation for the 21st Century p. 12

The Bigger Problem Absence of evidence is not evidence of absence. Don Rumsfeld. The fact that the data are consistent with null does not mean the null is true. Fighting this is an uphill battle. Mediation for the 21st Century p. 12

Recommendations Say what you know, not what you don t know. Give the range of values covered by the confidence intervals and their substantive interpretations. In this context, that means focusing on how much mediation there is, rather than whether there is mediation. If there really is no mediation, no amount of data will ever suffice to tell us the effect is 0. Data can tell us the effect is almost certainly trivial. Mediation for the 21st Century p. 13

An Example Conference Abstract, first draft: Background The deleterious effects of racism on a wide range of health outcomes including HIV risk is well documented among racial and ethnic minority groups in the United States. However, little is known about how MSM of color cope with stress from racism and whether coping with racism moderates the association between stress from racism and HIV risk among these men. Methods... Mediation for the 21st Century p. 14

Example, cont d Results... None of the interactions of stress with race/ethnicity, the four coping measures with race/ethnicity, and stress with the four coping measures was statistically significant. Conclusions Stress from racism within the gay community increased the likelihood of engaging in unprotected anal intercourse among U.S. MSM of color, but this association was not moderated by coping responses to racism.... Mediation for the 21st Century p. 15

Example, cont d Results... None of the interactions of stress with race/ethnicity, the four coping measures with race/ethnicity, and stress with the four coping measures was statistically significant. Conclusions Stress from racism within the gay community increased the likelihood of engaging in unprotected anal intercourse among U.S. MSM of color, but this association was not moderated by coping responses to racism.... 1. Strike all references to moderation. Mediation for the 21st Century p. 15

Example, cont d Results... None of the interactions of stress with race/ethnicity, the four coping measures with race/ethnicity, and stress with the four coping measures was statistically significant. Conclusions Stress from racism within the gay community increased the likelihood of engaging in unprotected anal intercourse among U.S. MSM of color, but this association was not moderated by coping responses to racism.... 1. Strike all references to moderation. 2. We can not draw any firm conclusions about moderation. Mediation for the 21st Century p. 15

Example concluded The wording as submitted: Background... However, little is known about how MSM of color cope with stress from racism and whether coping with racism buffers the impact of stress from racism on HIV risk among these men. Conclusions... However, we found little evidence that coping responses to racism buffered stress from racism.... Mediation for the 21st Century p. 16

Discussion How do you think statistical significance should be handled? Mediation for the 21st Century p. 17

Moderation: A Slippery Concept When the model is non-linear, it s not clear moderation is a meaningful concept. Suppose the true model for a binary outcome is the logistic: E [ log ( p 1 p )] = 1 + 3X + 4Y + 0XY. Since the interaction term is 0, there appears to be no moderation. Mediation for the 21st Century p. 18

Moderation: A Slippery Concept Y X 0 1 Effect of X 0 27% 50 23 1 50 88 38 Mediation for the 21st Century p. 19

Moderation: A Slippery Concept Y X 0 1 Effect of X 0 27% 50 23 1 50 88 38 The effect of X on the outcome varies depending on the level of Y, and we have moderation by the usual definition, even though the XY coefficient is 0. Mediation for the 21st Century p. 19

Moderation and Nonlinear Models This line of argument implies that any conceivable coefficients on X and Y, except ones that are exactly 0, produce moderation. This is likely to be true for most nonlinear models. So it does not seem useful to be worrying about whether there is or is not moderation. Mediation for the 21st Century p. 20

Structural Causal Models Recent work in Causal Inference (e.g., Pearl, 2010) has sought to clarify and generalize previous work on causation. Stochastic Effects operate by changing the probability distributions of outcomes. Non-Parametric The critical assumptions concern the absence of certain causal relationships. No functional forms are assumed. General Works with non-linear relations and discrete and continuous variables, both as effects and causes. Causal Inference Concerned with conditions under which we can extract causal relations from the world. Vocabulary Either graphs or functions. Mediation for the 21st Century p. 21

News You Can Use Good News We can think about and calculate mediation for any kind of model. Mediation for the 21st Century p. 22

News You Can Use Good News We can think about and calculate mediation for any kind of model. Bad News We have to be much more precise about exactly what questions we are asking. Mediation for the 21st Century p. 22

News You Can Use Good News We can think about and calculate mediation for any kind of model. Bad News We have to be much more precise about exactly what questions we are asking. News No magic: identification of true causal effects remains challenging to impossible. Mediation for the 21st Century p. 22

SCM Vocabulary and Concepts m = f M (x,µ M ) is a claim that X and a random term determine the value of M. Crucially, it claims that Y does not determine the value of M. f M (do(x 0 ),µ M ) is a statement about the distribution of M if we intervened and set X to x 0. An alternate notation is the potential outcome notation: M x0 (u) refers to the value individual u would have had if X had been x 0. Mediation for the 21st Century p. 23

Controlled Direct Effects CDE(m) = E(Y X = 1,M = m) E(Y X = 0,M = m). This is the effect of a 1 unit change in X on Y when the mediator is fixed at m. The expression assumes independence of the error terms for X and M. Without that assumption the more general expression is CDE(m) = E[Y do(x = 1,M = m)] E(Y do(x = 0,M = m)). This is the effect of X if we fix the mediator at m. Mediation for the 21st Century p. 24

Natural Direct Effects NDE x,x (Y ) = m { E(Y X = x,m = m) or, more compactly, E(Y X = x,m = m)} Pr(M = m X = x), NDE x,x = m { E(Y x,m) E(Y x,m) } Pr(m x). The NDE is the weighted average of the CDE s using the baseline X = x distribution of M as the weights. Mediation for the 21st Century p. 25

Indirect Effects IE x,x = m E(Y x,m) [ Pr(m x ) Pr(m x) ] X operates here only through its effects on M. Note that we evaluate those effects at the baseline level of X. Mediation for the 21st Century p. 26

Indirect Effects IE x,x = m E(Y x,m) [ Pr(m x ) Pr(m x) ] X operates here only through its effects on M. Note that we evaluate those effects at the baseline level of X. m E(Y x,m) Pr(m x ) is the expected value of Y for group x when the mediator has been changed to the levels for group x. Mediation for the 21st Century p. 26

Indirect Effects IE x,x = m E(Y x,m) [ Pr(m x ) Pr(m x) ] X operates here only through its effects on M. Note that we evaluate those effects at the baseline level of X. m E(Y x,m) Pr(m x ) is the expected value of Y for group x when the mediator has been changed to the levels for group x. m E(Y x,m) Pr(m x) is the expected value of Y for group x when the mediator has been left at the levels for group x. Mediation for the 21st Century p. 26

Total Effects Both the NDE and IE evaluate from a baseline of X = x, but the total effect has 3 sources: 1. The direct effect of raising X, holding M constant (NDE) 2. The indirect of effect of raising X on M, holding baseline X constant (IE). 3. The interaction of the changed M values with the changed X values. So we conclude that in general TE NDE + IE. Mediation for the 21st Century p. 27

Relation Between Effects AA F J API X AA IE(AA, API) 1 2 TE NDE(AA, API) API 3 4 Mediation for the 21st Century p. 28

Total and Reverse Effect Surprisingly, the total effect appears in an expression that includes travel backward along the paths. By going backward we get a term that starts from the interaction and steps down from it. TE x,x (Y ) = NDE x,x (Y ) IE x,x(y ). Mediation for the 21st Century p. 29

Caution: Multivariate Dangers AA F J API X AA 1 2 API 3 4 In a multivariate model, going from 1 to 2 to 4 does NOT turn AA s into API s. It turns AA s into API s who still have the AA distribution of other variables. The total effect is not the same as observed or modeled group differences, since those depend on the other covariates as well. Mediation for the 21st Century p. 30

No Estimation of Damaged Models None of the previous analysis relied on a statistical estimate of the gross effect of race on HIV status. A lot of the traditional procedures rely on estimating models with and without mediators. The SCM approach does not; it draws out the implications of the true model by manipulating it. Mediation for the 21st Century p. 31

Not Quite Home The SCM tools let us compare the direct and indirect effects of, e.g., being AA on an API. But we wanted to know what share of the differences between the three groups owed to indirect effects of jail. Mediation for the 21st Century p. 32

An Extension Define a target measure G(F) = g i,g j Ȳg i Ȳg j. This is the sum of the the absolute differences in means between all pairs of groups in a population with distribution F. Mediation for the 21st Century p. 33

An Extension Define a target measure G(F) = g i,g j Ȳg i Ȳg j. This is the sum of the the absolute differences in means between all pairs of groups in a population with distribution F. T is the distribution of the population under the model m = F M (x), y = F Y (x,m). This is the Theoretical distribution under the model. Mediation for the 21st Century p. 33

An Extension Define a target measure G(F) = g i,g j Ȳg i Ȳg j. This is the sum of the the absolute differences in means between all pairs of groups in a population with distribution F. T is the distribution of the population under the model m = F M (x), y = F Y (x,m). This is the Theoretical distribution under the model. A is an alternate distribution of the population under the model m = F m (do( )), y = F y (x,m). In this model all groups have the same distribution of Jail. do( ) means we use the overall population distribution, disregarding race. Mediation for the 21st Century p. 33

The Answer G(T) G(A) G(T) is the fraction of group differences mediated by Jail. Mediation for the 21st Century p. 34

Conclusion This is a useful, general purpose approach. Try it! Mediation for the 21st Century p. 35