Transcript Clash of Causal Inference Techniques

CLASH OF CAUSAL
INFERENCE TECHNIQUES
TERRENCE WILLETT AND NATHAN PELLEGRIN
RP GROUP CONFERENCE
KELLOGG WEST APRIL 2014
OUTCOMES
• Describe purpose of regression and propensity score
matching (PSM)
• Explain data requirements and basic procedure of
regression and PSM
• Compare and contrast regression and PSM
• Identify additional resources for further exploration
WHY CAUSAL INFERENCE
• If you need to use statistics, then you should design a
better experiment. –attributed to Rutherford
• Most education research is observational/correlational,
not experimental.
COMMON SCENARIO
•
Did participation in an activity, class or support result in higher outcomes
for students than would have happened had they not participated?
•
Students self-selected to participate and/or were recruited to
participate.
•
When participants are compared to non-participants, differences in
outcomes can be attributed to differences in background variables or
motivation.
•
•
Can we determine if the participation caused a change in outcomes?
No, but…
REGRESSION
𝑦 = 𝑎0 + 𝑎1 𝑥1 + 𝑎2 𝑥2 …+𝑎𝑛 𝑥𝑛 + 𝑒
• Classic correlational technique
• Covariates used in model to attempt to control for differences in
background variables or motivation
• Background variables can include measures of or proxies for skill
level, social capital, or socio-economic status
• Measures of self-motivation often unavailable
• Models are imperfect and generally must be combined with other
evidence to more completely describe the possible influence of an
intervention, program or strategy
PROPENSITY SCORE MATCHING
• One of several ways to create a matched comparison
group of non-participants intended to be similar to
participant group for a valid comparison
• Likelihood or other techniques used to create a score
indicating the likelihood that a particular non-participant
would have been a participant based on similarity to one
or more participants
THE COUNTERFACTUAL
(POTENTIAL OUTCOMES) FRAMEWORK
∆ = 𝑌𝑡 − 𝑌𝑐
This shows the causal effect of some event (participation,
treatment, intervention) on an individual.
𝑌 𝑡 is the potential outcome under treatment.
𝑌 𝑐 is the potential outcome under non-treatment (control,
counterfactual) and thus cannot be observed in our universe*.
*In light of recently published evidence for inflationary theories of the
cosmos, we note that we may be living in a multiverse (Alan Guth, 2014).
Counterfactual conditions may obtain in alternate universes.
THE POTENTIAL OUTCOMES MATRIX
Potential Outcome
Yt
Actual
Treatment
Status
(T)
1
0
= observable
= not observable
Yc
AVERAGE TREATMENT EFFECT (ATE)
•
Participants and non-participants differ systematically (w.r.t.
demographics, trajectories, risk profiles, self-selection, etc.)
•
Different people respond differently to treatment (differential
response)
•
These facts must be taken into account when
modeling/computing treatment effects.
•
This means all four cells of the matrix must be estimated in
order to obtain an average treatment effect !
•
How ?
SYMBOLIC DERIVATION OF
AVERAGE EFFECTS
ATE
1
∆ = 𝑌𝑖 𝑡 − 𝑌𝑖 𝑐
2
𝑌𝑖 = 𝑇𝑖 𝑌𝑖 𝑡 − (1 − 𝑇𝑖 )𝑌𝑖 𝑐
3
𝐸[∆] = 𝐸[𝑌𝑖 𝑡] − 𝐸[𝑌𝑖 𝑐]
4
𝐸𝑁 [𝑌𝑖 |𝑇𝑖 = 1] − 𝐸𝑁 [𝑌𝑖 |𝑇𝑖 = 0]
5
𝐸 ∆ = 𝜌𝐸 𝑌𝑖 𝑡 𝑇𝑖 = 1 + 1 − 𝜌 𝐸 𝑌𝑖 𝑡 𝑇𝑖 = 0
−{𝜌𝐸 𝑌𝑖 𝑐 𝑇𝑖 = 1 + 1 − 𝜌 𝐸 𝑌𝑖 𝑐 𝑇𝑖 = 0 }
ATT
6
𝐸 ∆ 𝑇𝑖 = 1 = 𝐸 𝑌𝑖 𝑡 𝑇𝑖 = 1 − 𝐸 𝑌𝑖 𝑐 𝑇𝑖 = 1
ATU
7
𝐸 ∆ 𝑇𝑖 = 0 = 𝐸 𝑌𝑖 𝑡 𝑇𝑖 = 0 − 𝐸 𝑌𝑖 𝑐 𝑇𝑖 = 0
CONDITIONAL INDEPENDENCE
(PERFECT STRATIFICATION)
(SELECTION ON OBSERVABLES)
• If our observations include information on every one of
the variables influencing likelihood of participation or
differential responses*, then it is possible to avoid
omitted variable bias and so achieve CI (PS, SOO). And
if we have CI, this is like randomized assignment ….
𝐸 𝑌 𝑡 𝑇 = 1, 𝑿 = 𝐸 𝑌 𝑡 𝑇 = 0, 𝑿
𝐸 𝑌 𝑐 𝑇 = 1, 𝑿 = 𝐸 𝑌 𝑐 𝑇 = 0, 𝑿
* How often do we encounter such datasets in institutional research?
ASSUMING CI…
Potential Outcome
Yt
Actual
Treatment
Status
(T)
Y(t)
1
0
Yc
Y(c)
= observable
= not observable
AVERAGE TREATMENT EFFECT (ATE)
• PSM offers a way to plug values into (5), (6) and (7) to
obtain unbiased estimates.
• However, if we have CI, then why not just something
like….
𝑌𝑖 = 𝛼 + 𝛾𝑇𝑖 + 𝜷′𝑿𝑖 + 𝜀𝑖 ?
(and many other regression techniques)
There is no firm and fast answer to this question. Decisions
are based on pragmatic considerations. However…..
EXAMPLE
PROS AND CONS
• Regressions can be “easier” to run but harder to explain to a general audience
• PSM can be more time consuming to conduct but easier to explain to a general
audience
• Regressions tend to perform better with large data sets while PSM tends to
perform better with few observations provided the non-participant group has
sufficient numbers of individuals with the key confounding variables
• Regressions have been used for many years and are well described
mathematically with broad consensus on proper error terms
• PSM is newer and there is not consensus on optimal matching procedures or
proper error terms
• Regression will use all cases with non-missing data while PSM may only a subset
of cases from the pool of non-participants
• All analytic methods suffer if key variables are not available
• Conclusions can often be the same with either method
HOW TO RUN PSM
•
•
Create data file (95% of effort)
Match participants and non-participants on a set of control variables to
create a comparison group with similar proportions on all characteristics
(i.e. comparison group would have a similar percent female, Hispanic,
low income, etc. as compared to the participant group)
•
This step is referred to as “balancing” and generally must be repeated several
times to obtain balance on all variables of interest either by adjusting matching
criteria or removing variables
•
Run comparative analyses, which can include simple t-tests, post-PSM
regressions, or other techniques
•
Major packages that conduct PSM include STATA, R, and SAS
•
•
STATA version 12 and older have psmatch2, v13 has teffects psmatch
Note SPSS/PASW does not do PSM directly but there is an R plugin for SPSS
• http://arxiv.org/ftp/arxiv/papers/1201/1201.6385.pdf
AN ALTERNATIVE STRATEGY TO DECIDING
THIS BATTLE
(WARNING: EMPIRICISTS APPROACHING)
ALTERNATIVE PERSPECTIVES
• Estimating program effects based on observational data can also be understood as…
•
•
•
•
Delimiting the error “built-into” inductive reasoning about causes (the problem of induction)
An “inverse problem” in the study of social dynamics.
A missing data problem
An optimization problem
• From these alternate perspectives there is a large menu of methods and extensions, including Euclidean,
Mahalanobis and Gower’s distance, nearest neighbor, cosine similarity, kernal functions, genetic algorithms,
imputation, dimensionality reduction, and other supervised/unsupervised learning algorithms
• A criterion that can be applied to regression, PSM and other methods is this: how do they perform at
predicting new (future) observations? (false positives, false negatives, correlated errors)
• In empirical investigations we do not answer this question just once using new replications or sampling
under varying conditions. It is applicable in any situation where predictions are made concerning new
(possibly future) observations. We update our models based on new evidence. And in this respect,
regression and PSM methods can both be used as tools of discovery; as ways to extend our understanding
of the processes producing patterns we find in sets of observations (and in streams of information,
generally).
• SO: CHOOSE THE METHOD/MODEL WHICH YIELDS THE SMALLEST PREDICTION ERROR.
• That may decide a battle in a particular setting (or occasion), but the war between methods will go on ….
“The inability to predict outliers implies the
inability to predict the course of history”
― Nassim Nicholas Taleb, The Black Swan: The
Impact of the Highly Improbable
“If you insist on strict proof (or strict disproof) in the
empirical sciences, you will never benefit from
experience, and never learn from it how wrong you
are.“
— Karl Raimund Popper, The Logic of Scientific
Discovery: Logik Der Forschung (2002), 28.
FURTHER READING
•
Angrist, J. D., & Pischke, J. (2008). Mostly Harmless Econometrics: An Empiricists
Companion
•
Morgan, S., Harding, D. (2006) Matching Estimators of Causal Effects: From
Stratification and Weighting to Practical Data Analysis Routines
•
Caliendo and Kopeinig. 2005. Practical Guide for PSM
•
www.caliendo.de/Papers/practical_revised_DP.pdf
•
Rosenbaum, P. R., & Rubin, D. B. (1983). The central role of the propensity score in
observational studies for causal effects. Biometrika, 70, 41-55.
•
Padgett, R.; Salisbury, M.; An, B.; & Pascarella, E. (2010). Required, Practical, or
Unnecessary? An Examination and Demonstration of Propensity Score Matching Using
Longitudinal Secondary Data. New Directions for Institutional Research – Assessment
Supplement (pp. 29-42). San Francisco, CA: Jossey-Bass.
•
Soledad Cepeda, M.; Boston, R.; Farrar, J., & Strom, B. (2003). Comparison of Logistic
Regression versus Propensity Score When the Number of Events Is Low and There Are
Multiple Confounders. American Journal of Epidemiology, 158, 280-287.
•
http://www.biostat.jhsph.edu/~estuart/propensityscoresoftware.html
THANK YOU
Terrence Willett
Director of Planning, Research, and Knowledge Systems
Cabrillo College
[email protected]
Nathan Pellegrin
Data Processing Specialist
Peralta District
[email protected]