Slide13. Logistic Regression

Multivariable Binary Logistic Regression
Overview:
Basic Concept of subgroup analysis assessing role of 3rd variable
Detecting confounding and interaction based on sub-group analyses
Detecting confounding and interaction using SPSS
1
22.3 Conducting subgroup analysis assessing the role of the 3rd factor
in SPSS
The VA randomized trial of coronary artery bypass surgery (CABG) vs.
standard medical therapy in patients with stable angina pectoris.
Dataset CABG.sav contains the following variables:
•
•
•
Primary study endpoint:
– All cause mortality at 6 years (Status6) (dead=1 / alive=0)
Primary study exposure
– Treatment (Surgery=1; Medical Treatment=0)
Covariates
– Left main disease (LMD) (yes=1/no=0)
– Angiographic risk (Angrisk) (low=0/high=1)
– Clinical risk (Clinrisk) (mild=1/moderate=2/ high=3)
– Diabetes (Diabetes) (yes=1/ no=0)
Peduzzi P, Kamina A, Detre K. Twenty-two-year follow-up in the VA Cooperative Study of Coronary
Artery Bypass Surgery for Stable Angina. Am J Cardiol (1999) 83 (2) 301-4
2
Assessing over all effect of surgery (all patients).
3
Crude (unadjusted) effect of CABG (All patients)
Over-all (unadjusted) OR for surgery and mortality = 0.73, p=0.081
indicating that the effect of surgery is marginally significant to improve
mortality.
4
Sub-group analysis by LMD
(with or without left main disease)
5
Among patients with
LMD (LMD=1)
Among patients with
No LMD (LMD=0)
a
Mortality at 6 years * CABG Crosstabulation
Count
Mortality at
6 years
Total
Count
Treatment
CABG
Alive
Died
Non CABG
235
76
311
a. Left Main Disease = Absent
CABG
223
61
284
a
Mortality at 6 years * CABG Crosstabulation
Total
458
137
595
Mortality at
6 years
Total
Treatment
CABG
Alive
Died
Non CABG
23
20
43
CABG
38
10
48
Total
61
30
91
a. Left Main Disease = Present
Different?
It appears that surgery improved mortality only among patients with LMD vs those without LMD
So let’s seek an evidence for the interaction (p-value for interaction) to numerically assess
ORLMD = 0.846 differs from ORno-LMD =0.303.
6
Testing for interaction using a binary logistic regression
ORno-LMD
ORint
P=0.043 for the interaction
indicates that LMD seems to modify association of surgery and
Mortality. Let’s calculate two stratified OR’s of mortality and surgery
by LMD based on this result.
7
How to interpret Exp(B) of a logistic regression with an interaction term
ORno-LMD
ORint
Exp(B) for Treatment among patients without LMD= 0.846 =ORno-LMD
Exp(B) for Treatment among patients with lmd = ORLMD
= ORno-MD X ORint
= 0.846 * 0.358 = 0.303. (this matches with the result of the sub-group
analysis)
8
Now let’s see what happen, if you adjust OR for LMD instead of
detecting the interaction.
This is LDM
adjusted OR of
death by
treatment.
9
Checking for confounding
•
Definition of confounder: Association of outcome and exposure
differs when another variable (covariate) is considered (i.e.,
adjusted).
Thus in order to assess if the third factor is confounding factor, we will
compare adjusted OR and un-adjusted OR.
10
Using Logistic Regression to detect confounding by the effect of
LMD
Model excluding LMD
Unadjusted OR of death by treatment = 0.731
Model adjusting for LMD
OR of death by treatment adjusted for LMD= 0.720
11
Does Diabetes Confound the Effect of Mortality and Surgery by
comparing crude and M-H adjusted OR’s?
Two OR’s are similar, thus we conclude that LMD is not confounding
the effect of treatment (surgery vs non-surgery).
Note: This is in fact expected, because in RCT, randomization provided
similar proportion of patients with LMD between surgery and control
groups, thus, LMD is not associated with exposure. In order to be a
confounder, it must associates with both exposure and outcome.
Confounding is not much of an issue in RCT, but interaction can well
exist.
12
Example of assessing an interaction
13