1 General Purpose
In the Chaps. 36 and 37 logistic regression with a single
binary or continuous predictor was explained. Just like linear
regression, logistic regression can also be performed on data with
multiple predictors. In this way the effects on the outcome of not
only treatment modalities, but also of additional predictors like
age, gender, comorbidities etc. can be tested simultaneously.
2 Schematic Overview of Type of Data File
3 Primary Scientific Question
Do all of the predictors independently
of one another predict the outcome.
4 Data Example
The example of Chap. 35 is used once more. In 55
hospitalized patients the risk of falling out of bed was assessed.
The question to be answered was: is there a significant difference
between the risk of falling out of bed at the departments of
surgery and internal medicine. The first 10 patients of the 55
patient file is underneath.
Fall
|
Dept
|
Age
|
Gender
|
Lett of complaint
|
1,00
|
,00
|
60,00
|
,00
|
1,00
|
1,00
|
,00
|
86,00
|
,00
|
1,00
|
1,00
|
,00
|
67,00
|
1,00
|
1,00
|
1,00
|
,00
|
75,00
|
,00
|
1,00
|
1,00
|
,00
|
56,00
|
1,00
|
1,00
|
1,00
|
,00
|
46,00
|
1,00
|
1,00
|
1,00
|
,00
|
98,00
|
,00
|
,00
|
1,00
|
,00
|
66,00
|
1,00
|
,00
|
1,00
|
,00
|
54,00
|
,00
|
,00
|
1,00
|
,00
|
86,00
|
1,00
|
1,00
|
5 Multiple Logistic Regression
The entire data file is entitled
“chapter35unpairedbinary” and is in extras.springer.com. We will
start by opening the data file in SPSS. First, simple logistic
regression with department as predictor and falloutofbed as outcome
will be performed. For analysis the statistical model Binary
Logistic Regression in the module Regression is required.
Command:
-
Analyze....Regression....Binary Logistic Regression....Dependent: enter falloutofbed....Covariates: enter department....click OK.
Variables in the equation
B
|
S.E.
|
Wald
|
df
|
Sig.
|
Exp(B)
|
||
---|---|---|---|---|---|---|---|
Step 1a
|
Department
|
1,386
|
,619
|
5,013
|
1
|
,025
|
4,000
|
Constant
|
−,288
|
,342
|
,709
|
1
|
,400
|
,750
|
The above results table of the logistic
regression shows that the department is a significant predictor at
p = 0,025.
Next, we will test whether age is a
significant predictor of falloutofbed.
Variables in the equation
B
|
S.E.
|
Wald
|
df
|
Sig.
|
Exp(B)
|
||
---|---|---|---|---|---|---|---|
Step 1a
|
Age
|
,106
|
,027
|
15,363
|
1
|
,000
|
1,112
|
Constant
|
−6,442
|
1,718
|
14,068
|
1
|
,000
|
,002
|
Also age is a significant predictor of
falling out of bed at p < 0,0001.
Subsequently, we will test all of the
predictors simultaneously, and, in addition, will test the
possibility of interaction between age and department on the
outcome. Clinically, this could very well exist. Therefore, we will
add an interaction-variable of the two as an additional predictor.
Command:
-
Analyze....Regression....Binary Logistic Regression....Dependent: falloutofbed.... Covariates: age, department, gender, lettereof complaint, and interaction variable “age by department” (click for that “ > a*b > ”in the dialog window)....click OK.
Variables in the equation
B
|
S.E.
|
Wald
|
df
|
Sig.
|
Exp(B)
|
||
---|---|---|---|---|---|---|---|
Step 1a
|
Age
|
,067
|
,028
|
5,830
|
1
|
,016
|
1,069
|
Department
|
−276,305
|
43760,659
|
,000
|
1
|
,995
|
,000
|
|
Gender
|
,235
|
1,031
|
,052
|
1
|
,819
|
1,265
|
|
Letter complaint
|
1,582
|
1,036
|
2,331
|
1
|
,127
|
4,862
|
|
Age by department
|
4,579
|
720,744
|
,000
|
1
|
,995
|
97,447
|
|
Constant
|
−4,971
|
1,891
|
6,909
|
1
|
,009
|
,007
|
The above table shows the output of
the multiple logistic regression. Interaction is not observed, and
the significant effect of the department has disappeared, while age
as single variable is a statistically significant predictor of
falling out of bed with a p-value of 0,016 and an odds ratio of
1,069 per year.
The initial significant effect of the
difference in department is, obviously, not caused by a real
difference, but rather by the fact that at one department many more
elderly patients had been admitted than those at the other
department. After adjustment for age the significant effect of the
department had disappeared.
6 Conclusion
In the Chaps. 36 and 37 logistic regression with a single
binary or continuous predictor was explained. Just like linear
regression, logistic regression can also be performed on data with
multiple predictors. In this way the effects on the outcome of not
only treatment modalities, but also of additional predictors like
age, gender, comorbidities etc. can be tested simultaneously. If
you have clinical arguments for interactions, then interaction
variables can be added to the data. The above analysis shows that
department was a confounder rather than a real effect (Confounding
is reviewed in the Chap. 22).
7 Note
More background, theoretical, and
mathematical information about logistic regression is given in
Statistics applied to clinical studies 5th edition, Chaps. 17 and
65, Springer Heidelberg Germany, 2012, from the same authors.