novel approaches to adjusting for confounding: propensity scores, instrumental variables and msms...
TRANSCRIPT
![Page 1: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/1.jpg)
Novel Approaches to Adjusting for Confounding:
Propensity Scores, Instrumental Variables and MSMs
Matthew Fox
Advanced Epidemiology
![Page 2: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/2.jpg)
What are the exposures you are interested in studying?
![Page 3: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/3.jpg)
Assuming I could guarantee you that you would not create bias,
which approach is better: randomization or adjustment for
every known variable?
![Page 4: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/4.jpg)
What is intention to treat analysis?
![Page 5: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/5.jpg)
Yesterday
Causal diagrams (DAGS)– Discussed rules of DAGS– Goes beyond statistical methods and
forces us to use prior causal knowledge– Teaches us adjustment can CREATE bias
Helps identify a sufficient set of confounders– Not how to adjust for them
![Page 6: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/6.jpg)
This week
Beyond stratification and regression– New approaches to adjusting for (not
“controlling” ) confounding– Instrumental variables– Propensity scores (Confounder scores)– Marginal structural models
Time dependent confounding
![Page 7: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/7.jpg)
Given the problems with the odds ratio, why does everyone use it?
![Page 8: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/8.jpg)
Non-collapsibility of the OR (Excel; SAS)
Odds ratio collapsibility but confounding
C+ C- Total
E+ E- E+ E- E+ E-
Disease+ 400 300 240 180 640 480
Disease- 100 200 360 720 460 920
Total 500 500 600 900 1100 1400
Risk 0.80 0.60 0.40 0.20 0.58 0.34
Odds 4.00 1.50 0.67 0.25 1.39 0.52
Crude Adj
RR 1.33 2.00 1.6968 1.5496
OR 2.67 2.67 2.67 2.67
RD 0.2 0.2 0.23896 0.2
![Page 9: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/9.jpg)
![Page 10: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/10.jpg)
Solution: SAS Code
title "Crude relative risk model";proc genmod data=odds descending;
model d = e/link=log dist=bin;run;
title "Adjusted relative risk model";proc genmod data=odds descending;
model d = e c/link=log dist=bin;run;
![Page 11: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/11.jpg)
Model crude:
Exp(0.5288) = 1.6968Crude RR was 1.6968
Results
![Page 12: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/12.jpg)
Results
Model adjusted: Exp(0.3794) = 1.461
MH RR was 1.55
![Page 13: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/13.jpg)
STATA
glm d e, family(binomial) link(log) glm d e c, family(binomial) link(log)
![Page 14: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/14.jpg)
What about risk differences?
![Page 15: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/15.jpg)
![Page 16: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/16.jpg)
Solution: SAS Code
title "Crude risk differences model";proc genmod data=odds descending;
model d = e/link=bin dist=identity;run;
title "Adjusted risk differences model";proc genmod data=odds descending;
model d = e c/link=bin dist=identity;run;
![Page 17: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/17.jpg)
Model crude:
Exp(0.5288) = 1.6968Crude RR was 1.7
Results
Model crude: 0.239
Crude RD = 0.23896
![Page 18: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/18.jpg)
Results
Adjusted model : 0.20
MH RD = 0.20
![Page 19: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/19.jpg)
STATA
glm d e, family(binomial) link(identity) glm d e c, family(binomial) link(identity)
glm d e c c*e, family(binomial) link(identity)
![Page 20: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/20.jpg)
Novel approaches to controlling confounding
![Page 21: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/21.jpg)
Limitations of Stratification and Regression
Stratification/regression work well with point exposures with complete follow up and sufficient data to adjust– Limited data on confounders or small cells– No counterfactual for some people in our dataset
Regression often estimates parameters– Time dependent exposures and confounding
A common situationWith time dependence, DAGs gets complex
![Page 22: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/22.jpg)
Randomization and Counterfactuals
Ideally, evidence comes from RCTs– Randomization gives expectation unexposed can
stand in for the counterfactual ideal Full exchangeability: E(p1=q1, p2=q2, p3=q3, p4=q4)
– In expectation, assuming no other bias [Pr(Ya=1=1) - Pr(Ya=0=1)] = [Pr(Y=1|A=1) - Pr(Y=1|A=0)]
Since we assign A, RRAC = 1– If we can’t randomize, what can we do to
approximate randomization?
![Page 23: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/23.jpg)
How randomization works
Randomized Controlled Trial
Randomization strongly predicts exposure (ITT)
C2
C1
C3
A DRandomization
![Page 24: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/24.jpg)
A typical observational study
Observational Study
C2
C1
C3
A D?
![Page 25: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/25.jpg)
A typical observational study
Observational Study
C2
C1
C3
A D
Regression/stratification seeks to block backdoor path from A to D by averaging A-D associations within levels of Cx
![Page 26: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/26.jpg)
Approach 1: Instrumental Variables
![Page 27: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/27.jpg)
Intention to treat analysis
In an RCT we assign the exposure– e.g. assign people to take an aspirin a day vs. not– But not all will take aspirin when told to and others
will take it even if told not to
What to do with those who don’t “obey”?– The paradigm of intention to treat analysis says
analyze subject in the group they are assigned Maintains the benefits of randomization Biases towards the null at worst
![Page 28: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/28.jpg)
Instrumental variables
An approach to dealing with confounding using a single variable– Works along the same lines as randomization
Commonly used approach in economics, yet rarely used in medical research– Suggests we are either behind the times or they are
hard to find– Party privileged in economics because little
adjustment data exists
![Page 29: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/29.jpg)
Instrumental variables
An instrument (I):– A variable that satisfies 3 conditions:
Strongly associated with exposureHas no effect on outcome except through A (E)Shares no common causes with outcome
Ignore E-D relationship– Measure association between I and D
This is not confounded
– Approximates an ITT approach
![Page 30: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/30.jpg)
Adjust the IV estimate
Can optionally adjust IV estimate to estimate the effect of A (exposure)– But differs from randomization
If an instrument can be found, has the advantage we can adjust for unknown confounders– This is the benefit we get from randomization?
![Page 31: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/31.jpg)
Intention to Treat (IV Ex 1)
A(Exposure): Aspirin vs. Placebo Outcome: First MI Instrument: Randomized assignment
Therapy MI
Confounders
Randomization
Condition 1: Predictor of A ?
Condition 2: no direct effect on the outcome?
Condition 3: No common causes with outcome?
![Page 32: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/32.jpg)
Confounding by indication (IV Ex 2)
A(Exposure): COX2 inhibitor vs NSAID Outcome: GI complications Instrument: Physician’s previous prescription
COX2/NSAID GI comp
Indications
Previous Px
Regression (17 confounders), no effect RD: -0.06/100; 95% CI -0.26 to 0.14
IV: Protective effect of COX-2 RD: -1.31/100; -2.42 to -0.20
Compatible with trial results RD: -0.65/100; -1.08 to -0.22
![Page 33: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/33.jpg)
Unknown confounders (IV Ex 3)
A(Exposure): Childhood dehydration Outcome: Adult high blood pressure Instrument: 1st year summer climate
dehydration High BP
SES
1st year climate
Hypothesized hottest/driest summers in infancy would be associated with severe infant diarrhea/dehydration, and consequently higher blood pressure in adulthood.
For 3,964 women born 1919-1940, a 1 SD (1.3 ºC) > mean summer temp in 1st year life associated with 1.12-mmHg (95% CI: 0.33, 1.91) > adult systolic blood pressure, and 1 SD > mean summer rainfall (33.9 mm) associated with < systolic blood pressure (-1.65 mmHg, 95% CI: -2.44, -0.85).
![Page 34: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/34.jpg)
Optionally we can adjust for “non-compliance”
Optionally if we want to estimate A-D relationship, not I-D, we can adjust:– RDID / RDIE
– Inflate the IV estimator to adjust for the lack of perfect correlation between I and E
– If I perfectly predicts E then RDIE = 1, so adjustment does nothing
Like per protocol analysis– But adjusted for confounders
![Page 35: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/35.jpg)
To good to be true?
Maybe The assumptions needed for an
instrument are un-testable from the data– Can only determine if I is associated with A
Failure to meet the assumptions can cause strong bias– Particularly if we have a “weak” instrument
![Page 36: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/36.jpg)
Approach 2: Propensity Scores
![Page 37: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/37.jpg)
Comes out of a world of large datasets (Health insurance data)
Cases where we have a small (relative to the size of the dataset) exposed population and lots and lots of potential comparisons in the unexposed group– And lots of covariate data to adjust for
Then we have luxury of deciding who to include in study as a comparison group based on a counterfactual definition
![Page 38: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/38.jpg)
Propensity Score
Model each subject’s propensity to receive the index condition as a function of confounders– Model is independent of outcomes, so good for
rare disease, common exposure Use the propensity score to balance
assignment to index or reference by:– Matching– Stratification– Modeling
![Page 39: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/39.jpg)
Propensity Scores
The propensity score for subject i is:– Probability of being assigned to treatment A = 1 vs. reference
A = 0 given a vector xi of observed covariates:
In other words, the propensity score is:– Probability that the person got the exposure given anything
else we know about them
)|1(Pr iiiA xX
![Page 40: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/40.jpg)
Why estimate the probability a subject receives a certain treatment when it is known what treatment they received?
![Page 41: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/41.jpg)
How Propensity Scores Work
Quasi-experiment – Using probability a subject would have been treated
(propensity score) to adjust estimate of the treatment effect, we simulate a RCT
2 subjects with = propensity, one E+, one E- – We can think of these two as “randomly assigned”
to groups, since they have the same probability of being treated, given their covariates
– Assumes we have enough observed data that within levels of propensity E is truly random
![Page 42: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/42.jpg)
Propensity Scores:Smoking and Colon Cancer
Have info on people’s covariates:– Alcohol use, sex, weight, age, exercise, etc:
Person A is a smoker, B is not– Both had 85% predicted probability of smoking
If “propensity” to smoke is same, only difference is 1 smoked and 1 didn’t
– This is essentially what randomization does– B is the counterfactual for A assuming a correct
model for predicting smoking
![Page 43: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/43.jpg)
Obtaining Propensity Scores in SAS
Calculate propensity scoreproc logistic;
model exposure = cov_1 cov_2 … cov_n; output out = pscoredat pred = pscore;
run;
Either match subjects on propensity score or adjust for propensity score
proc logistic;model outcome = exposure pscore;
run;
![Page 44: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/44.jpg)
Pros and Cons of PS
Pros– Adjustment for 1 confounder– Allows estimation of the exposure and fitting a final
model without ever seeing the outcome– Allows us to see parts of data we really should not
be drawing conclusions on b/c no counterfactual Cons
– Only works if have good overlap in pscores– Does not fix conditioning on a collider problem– Doesn’t deal with unmeasured confounders
![Page 45: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/45.jpg)
Study of effect of neighborhood segregation on IMR
![Page 46: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/46.jpg)
Approach 3: Marginal Structural Models
![Page 47: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/47.jpg)
Time Dependent Confounding
Time dependent confounding:1) Time dependent covariate that is a risk
factor for or predictive of the outcome and also predicts subsequent exposure
Problematic if also:2) Past exposure history predicts
subsequent level of covariate
![Page 48: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/48.jpg)
Example
Observational study of subjects infected with HIV– E = HAART therapy– D = All cause mortality– C = CD4 count
![Page 49: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/49.jpg)
Time Dependent Confounding
A0 A1
C0 C1
D
A0 A1
C0 C1
D
1)
2)
![Page 50: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/50.jpg)
Failure of Traditional Methods
Want to estimate causal effect of A on D– Can’t stratify on C (it’s an intermediate)– Can’t ignore C (it’s a confounder)
Solution – rather than stratify, weight– Equivalent to standardization
Create pseudo-population where RRCE = 1– Weight each person by “inverse probability of
treatment” they actually received– Weighting doesn’t cause problems pooling did– In DAG, remove arrow C to A, don’t box
![Page 51: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/51.jpg)
Remember back to the SMRCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
0.2
100050
*50020020
*1500
50050
*5001500300
*1500
*
*
*
*
ˆ
01
11
0
1
Nb
N
Na
N
WNb
W
WNa
W
RSM
![Page 52: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/52.jpg)
The SMR asks, what if the exposed had also been unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 D+ 300 D+ 50D- 1650 D- 1200 D- 450Total 2000 Total 1500 Total 500Risk 0.18 0.2 0.1RR RR RR
SM
R
![Page 53: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/53.jpg)
The SMR asks, what if the exposed had also been unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 D+ 300 D+ 50D- 1650 D- 1200 D- 450Total 2000 Total 1500 1500Total 500 500Risk 0.18 0.2 0.1RR RR RR
SM
R
![Page 54: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/54.jpg)
The SMR asks, what if the exposed had also been unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 D+ 300 D+ 50D- 1650 D- 1200 D- 450Total 2000 Total 1500 1500Total 500 500Risk 0.18 0.2 0.1 0.1 0.05RR RR RR
SM
R
![Page 55: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/55.jpg)
The SMR asks, what if the exposed had also been unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 D+ 300 150 D+ 50 25D- 1650 D- 1200 1350D- 450 475Total 2000 Total 1500 1500Total 500 500Risk 0.18 0.2 0.1 0.1 0.05RR RR 2.0 RR 2.0
SM
R
![Page 56: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/56.jpg)
The SMR asks, what if the exposed had also been unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 175 D+ 300 150 D+ 50 25D- 1650 1825 D- 1200 1350D- 450 475Total 2000 2000 Total 1500 1500Total 500 500Risk 0.175 0.875 0.2 0.1 0.1 0.05RR 2.0 RR 2.0 RR 2.0
Crude now equals the adjusted. No need to adjust.
SM
R
![Page 57: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/57.jpg)
Could also ask, what if everyone was both exposed, unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0
![Page 58: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/58.jpg)
Could also ask, what if everyone was both exposed, unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ D+ D+D- D- D-Total Total 1700 1700Total 1500 1500Risk 0.2 0.1 0.1 0.05RR RR 2.0 RR 2.0
![Page 59: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/59.jpg)
Could also ask, what if everyone was both exposed, unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ D+ 340 170 D+ 150 75D- D- 1360 1530D- 1350 1425Total Total 1700 1700Total 1500 1500Risk 0.2 0.1 0.1 0.05RR RR 2.0 RR 2.0
![Page 60: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/60.jpg)
Could also ask, what if everyone was both exposed, unexposed?
Crude C1 C0
E+ E- E+ E- E+ E-D+ 490 245 D+ 340 170 D+ 150 75D- 2710 2955 D- 1360 1530D- 1350 1425Total 3200 3200 Total 1700 1700Total 1500 1500Risk 0.153 0.077 0.2 0.1 0.1 0.05RR 2.0 RR 2.0 RR 2.0
![Page 61: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/61.jpg)
What is Inverse Probability Weighting (IPW)?
Weight each subject by inverse probability of treatment received
Probability of treatment is:– p(receiving treatment received| covariates)– Adjust # of E+ and E- subjects in C strata
Weighting breaks E-C link only– Now Marginal (Crude) = Causal Effect
But that’s what we just did
![Page 62: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/62.jpg)
Calculate the weightsCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0 PTIPTW Calculate p(receiving treatment received|C) For C=1, E=1
– PT = 1500/1700 = 0.88– IPTW = 1/0.88 = 1.13
![Page 63: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/63.jpg)
Calculate the weightsCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0 PT 0.88IPTW 1.13 Calculate p(receiving treatment received|C) For C=1, E=1
– PT = 1500/1700 = 0.88– IPTW = 1/0.88 = 1.13
![Page 64: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/64.jpg)
Calculate the weightsCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0 PT 0.88 0.12IPTW 1.13 8.50 Calculate p(receiving treatment received|C) For C=1, E=0
– PT = 200/1700 = 0.12– IPTW = 1/0.12 = 8.50
![Page 65: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/65.jpg)
Calculate the weightsCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0 PT 0.88 0.12 0.33 0.67IPTW 1.13 8.50 3.00 1.50 Calculate p(receiving treatment received|C) For C=1, E=0
– PT = 200/1700 = 0.12– IPTW = 1/0.12 = 8.50
Multiply cell number by the weights
![Page 66: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/66.jpg)
Apply the weightsCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0 PT 0.88 0.12 0.33 0.67IPTW 1.13 8.50 3.00 1.50Pseudo population
Crude C1 C0
E+ E- E+ E- E+ E-D+ D+ 340 170 D+ 150 75D- D- 1360 1530 D- 1350 1425Total Total 1700 1700 Total 1500 1500Risk 0.2 0.1 0.1 0.05RR RR 2.0 RR 2.0
![Page 67: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/67.jpg)
CollapseCrude C1 C0
E+ E- E+ E- E+ E-D+ 350 70 D+ 300 20 D+ 50 50D- 1650 1130 D- 1200 180 D- 450 950Total 2000 1200 Total 1500 200 Total 500 1000Risk 0.18 0.06 0.2 0.1 0.1 0.05RR 3.0 RR 2.0 RR 2.0 PT 0.88 0.12 0.33 0.67IPTW 1.13 8.50 3.00 1.50Pseudo population
Crude C1 C0
E+ E- E+ E- E+ E-D+ 490 245 D+ 340 170 D+ 150 75D- 2710 2955 D- 1360 1530 D- 1350 1425Total 3200 3200 Total 1700 1700 Total 1500 1500Risk 0.153 0.077 0.2 0.1 0.1 0.05RR 2.0 RR 2.0 RR 2.0
Broke link between C and E without stratification, so no problem of conditioning on collider
![Page 68: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/68.jpg)
Pseudo-population
The “pseudo-population” breaks the link between the exposure and the outcome without stratification– Note this is different from stratifying – Create a standard population without confounding
By creating multiple copies of people, standard errors will be biased– Use robust standard errors to adjust
![Page 69: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/69.jpg)
Robins and Hernán
“The IPTW method effectively simulates the data that would be observed had, contrary to fact,
exposure been conditionally randomized”
![Page 70: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/70.jpg)
Time Dependent Confounding
Extend method to time dependent confounders– Predict p(receiving treatment actually received at time t1|
covariates, treatment at t0)
Probability of treatment at t1 is:
p(receiving treatment received at t0) *
p(receiving treatment received at t1)
See Hernán for SAS code, not hard scwgt command, robust SE (repeated statement)
![Page 71: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/71.jpg)
Time Dependent Confounding
E0 E1
C0 C1
D
E0 E1
C0 C1
D
1) Before IPTW
2) After IPTW
![Page 72: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/72.jpg)
![Page 73: Novel Approaches to Adjusting for Confounding: Propensity Scores, Instrumental Variables and MSMs Matthew Fox Advanced Epidemiology](https://reader035.vdocument.in/reader035/viewer/2022081512/56649e2d5503460f94b1c7d7/html5/thumbnails/73.jpg)
Limitations of MSMs
Very sensitive to weights Still need to be able to be able to predict
the exposure– The methods solves the structural problem, but
we still need the data to be able to accurately predict exposure
Still have to get the model right