Study population
A total of 1,236 patients with their spouses as unit were included in the analysis, 359 (29%) ICU survivor's patients whose spouses were trained as ICU-doulas, 461 (37.3%) ICU survivor's patients without a trained spouse, and 416 (33.7%) patients discharged from general wards of hospital without a trained spouse. Comparison of demographics and clinical characteristics between three study groups are presented in Table 1. Demographic characteristics such as age (P = 0.605), gender (P = 0.849), comorbidities (P = 0.221) were similar between three groups of study. Also, in terms of some variables such as high alarm silence strategy (P = 0.261), physical restraint (P = 0.196), sleep disturbance (P = 0.068) and the mean total of hospital length of stay [LOS] (P = 0.370) were not significant difference between two groups of study. However, in terms of severity of illness based on Acute Physiology and Chronic Health Evaluation (APACHE) IV score (P < 0.001), dose of sedation (P = 0.016) and nursing care (P = 0.002) a significant difference was observed between the three groups.
Table 1
Comparison of demographics and clinical characteristics between three study groups
Characteristics
|
ICU survivors with doulas (n = 359)
|
ICU survivors without doulas
(n = 461)
|
No ICU patients without doulas
(n = 416)
|
P-value
|
Age (years)
|
Mean ± SD
|
67.11 ± 11.08
|
67.70 ± 11.31
|
67.90 ± 11.51
|
0.605
|
Gender
|
Male (%)
|
159 (44.3)
|
212 (46)
|
192 (46.2)
|
0.849
|
|
Female (%)
|
200 (55.7)
|
249 (54)
|
224 (53.8)
|
|
Comorbidities
|
Yes (%)
|
85 (23.7)
|
87 (18.9)
|
83 (20)
|
0.221
|
|
No (%)
|
274 (76.3)
|
374 (81.1)
|
333 (80)
|
|
Comorbidities
types
|
Cardiovascular disease (%)
|
5 (1.4)
|
12 (2.6)
|
6 (1.4)
|
|
CHF/cardiomyopathy (%)
|
14 (3.9)
|
10 (2.2)
|
10 (2.4)
|
0.856
|
|
Diabetes (%)
|
14 (3.9)
|
13 (2.8)
|
10 (2.4)
|
|
|
Chronic kidney disease (%)
|
19 (5.3)
|
11 (2.4)
|
14 (3.4)
|
|
|
Chronic liver disease (%)
|
7 (1.9)
|
14 (3.0)
|
14 (3.4)
|
|
|
Chronic pulmonary disease (%)
|
6 (1.7)
|
3 (0.7)
|
8 (1.9)
|
|
|
Liver cancer (%)
|
8 (2.2)
|
4 (0.9)
|
5 (1.2)
|
|
|
Gastrointestinal cancer (%)
|
5 (1.4)
|
4 (0.9)
|
7 (1.7)
|
|
|
Brain cancer (%)
|
4 (1.1)
|
11 (2.4)
|
3 (0.7)
|
|
|
Hypertension (%)
|
2 (0.6)
|
3 (0.7)
|
3 (0.7)
|
|
|
Arthritis (%)
|
1 (0.3)
|
2 (0.4)
|
3 (0.7)
|
|
APACHE IV
|
Mean ± SD
|
15.69 ± 2.48
|
16.79 ± 2.45
|
14.72 ± 2.07
|
< 0.001*
|
|
(Range)
|
(11–22)
|
(15–22)
|
(11–19)
|
|
Physical restraint
|
Yes (%)
|
142 (39.6)
|
210 (45.6)
|
172 (41.3)
|
0.196
|
|
No (%)
|
217 (60.4)
|
251 (54.4)
|
244 (58.7)
|
|
Sedation
|
Low (%)
|
124 (34.5)
|
120 (26)
|
90 (21.6)
|
0.016*
|
|
Moderate (%)
|
187 (52.1)
|
249 (54)
|
240 (57.7)
|
|
|
High (%)
|
48 (13.4)
|
92 (20)
|
86 (20.7)
|
|
Nursing Care
|
Low (%)
|
132 (36.8)
|
125 (27.1)
|
118 (28.4)
|
0.002*
|
|
Moderate (%)
|
175 (48.7)
|
230 (49.9)
|
198 (47.6)
|
|
|
High (%)
|
52 (14.5)
|
106 (23)
|
100 (24)
|
|
Sleep disturbance
|
Yes (%)
|
256 (71.3)
|
300 (65.1)
|
266 (63.9)
|
0.068
|
|
No (%)
|
103 (28.7)
|
161 (34.9)
|
150 (36.1)
|
|
Alarm silence
|
Low (%)
|
33 (9.2)
|
33 (7.2)
|
35 (8.4)
|
0.261
|
|
Moderate (%)
|
289 (80.5)
|
363 (78.7)
|
330 (79.3)
|
|
|
High (%)
|
37 (10.3)
|
65 (14.1)
|
51 (12.3)
|
|
Hospital LOS
|
Mean ± SD
|
17.57 ± 11.17
|
16.70 ± 10.17
|
16.61 ± 9.92
|
0.370
|
|
(Range)
|
(3–95)
|
(2–92)
|
(2–70)
|
|
* P-value < 0.05 considered as significant, Abbreviations: APACHE IV means Acute Physiology and Chronic Health Evaluation IV, LOS means length of stay |
Comparison of outcomes between two groups of study
Comparison of outcomes based on WHO-5 index, HADS-depression/anxiety scores and IES-R score between the three study groups are presented in Table 2. According to the results, the mean WHO-5 index in the ICU survivors with doulas was significantly higher than the ICU survivors (68.57 ± 27.47 vs. 50.80 ± 15.75, P < 0.001) and no ICU patients without doulas (68.57 ± 27.47 vs. 50.42 ± 16.20, P < 0.001), which the higher scores indicate better well-being. The mean HADS-anxiety scores and HADS-depression scores in the ICU survivors with doulas was significantly lower than the both groups without doulas (P < 0.001), so that lower scores indicate normality. The mean scores of three subscales IES-R score; IES-R intrusion, IES-R avoidance and IES-R hyperarousal scores was significantly higher in both groups without doulas than the ICU survivors with doulas (P < 0.001). The mean total IES-R score was 41.58 ± 5.82, 50.48 ± 10.75 and 50.06 ± 10.33 in the ICU survivors with doulas, the ICU survivors without doulas and no ICU patients without doulas, respectively (≥ 33 signifies distress). IES-R score was significantly lower in the ICU survivors with doulas than the both groups without doulas (P < 0.001). However, there was no statistically significant difference in terms of WHO-5 index, HADS-depression/anxiety and IES-R scores between ICU survivors without doulas and no ICU patients without doulas (P > 0.05).
Table 2
Comparison of outcomes between three study groups
Outcomes
|
ICU survivors with doulas (n = 359)
|
ICU survivors without doulas
(n = 461)
|
No ICU patients without doulas
(n = 416)
|
P-value
|
WHO-5 index score
|
Mean ± SD
|
68.57 ± 27.47
|
50.80 ± 15.75
|
50.42 ± 16.20
|
< 0.001*
|
Median (IQR)
|
84 (40–92)
|
52 (36–64)
|
48 (36–64)
|
|
(Range)
|
(20–100)
|
(24–80)
|
(24–80)
|
|
WHO-5 index categorized
|
Low and moderate (%)
|
117 (32.6)
|
276 (59.9)
|
248 (59.6)
|
< 0.001*
|
High (%)
|
242 (67.4)
|
185 (40.1)
|
168 (40.4)
|
|
HADS depression score
|
Mean ± SD
|
7.18 ± 2.60
|
9.05 ± 4.30*
|
8.60 ± 3.95*
|
< 0.001*
|
Median (IQR)
|
7 (5–9)
|
8 (5.5–12)
|
8 (5–11)
|
|
(Range)
|
(3–14)
|
(3–21)
|
(3–21)
|
|
HADS depression categorized
|
Normal (< 8 score)
|
198 (55.2)
|
207 (44.9)
|
204 (49)
|
< 0.001*
|
Abnormal (≥ 8 score)
|
161 (44.8)
|
254 (55.1)
|
212 (51)
|
|
HADS anxiety score
|
Mean ± SD
|
8.96 ± 4.48
|
11.20 ± 5.13*
|
11.06 ± 4.98*
|
< 0.001*
|
Median (IQR)
|
8 (6–12)
|
11 (7–15)
|
11 (7.25-14)
|
|
(Range)
|
(3–21)
|
(3–21)
|
(3–21)
|
|
HADS anxiety categorized
|
Normal (< 8 score)
|
171 (47.6)
|
116 (25.2)
|
104 (25)
|
< 0.001*
|
Abnormal (≥ 8 score)
|
188 (52.4)
|
354 (78.4)
|
312 (75)
|
|
IES-R intrusion score
|
Mean ± SD
|
13.07 ± 3.19
|
14.75 ± 5.33
|
14.64 ± 5.10
|
< 0.001*
|
Median (IQR)
|
13 (11–15)
|
13 (11–18)
|
13 (11–18)
|
|
(Range)
|
(8–21)
|
(8–32)
|
(8–32)
|
|
IES-R avoidance score
|
Mean ± SD
|
15.37 ± 3.79
|
18.62 ± 6.42
|
18.39 ± 6.04
|
< 0.001*
|
Median (IQR)
|
15 (12–18)
|
17 (14–23)
|
17 (14–22)
|
|
(Range)
|
(8–32)
|
(8–32)
|
(8–32)
|
|
IES-R hyperarousal score
|
Mean ± SD
|
13.14 ± 3.43
|
17.10 ± 7.05
|
17.02 ± 7.02
|
< 0.001*
|
Median (IQR)
|
13 (11–16)
|
16 (11-23.5)
|
16 (11–23)
|
|
(Range)
|
(8–31)
|
(8–32)
|
(8–32)
|
|
IES-R total score
|
Mean ± SD
|
41.58 ± 5.82
|
50.48 ± 10.75
|
50.06 ± 10.33
|
< 0.001*
|
Median (IQR)
|
42 (38–46)
|
50 (42–58)
|
50 (42–57)
|
|
(Range)
|
(25–58)
|
(29–84)
|
(27–84)
|
|
IES-R total score categorized
|
Normal (< 33 score)
|
22 (6.1)
|
8 (1.7)
|
16 (3.8)
|
< 0.001*
|
Abnormal (≥ 33 score)
|
337 (93.9)
|
453 (98.3
|
400 (96.2)
|
|
* P-value < 0.05 considered as significant, World Health Organization- Five Well-Being (WHO-5) index that low scores (20–53) indicate distress and high scores (54–100) indicate well-being, Hospital Anxiety and Depression Scale (HADS) scores that cutoff < 8 scores indicates normality and ≥ 8 scores shows anxiety and depression disorders, Event Scale-Revised (IES-R) score that cutoff < 33 scores showed normality and ≥ 33 scores showed Post-traumatic stress disorder (PTSD) |
Multinomial logistic regression finding
Univariate and multivariate multinomial logistic regression analysis to find effects of variables on WHO-5 index, IER-S and HADS-depression/anxiety scores are presented in Table 3–6. Multivariate analysis showed that the risk of well-being over the mental distress were significantly decreased in patients with comorbidities (OR: 0.661, 95%CI: 0.46–0.94, P = 0.024), those required more nursing care (OR: 0.701, 95%CI: 0.491–0.99, P = 0.050), and who received the high dose of sedations (OR: 0.513, 95%CI: 0.35–0.74, P < 0.001). However, the OR well-being was significantly increased in the ICU survivors with doulas than the ICU survivors without doulas (OR: 3.614, 95% CI: 2.65–4.91, P < 0.001) and no ICU patients without doulas (OR: 2.231, 95% CI: 1.49–3.33, P < 0.001). However, the risk of medium score of WHO-5index over the low score of WHO-5 index was decreased in the ICU survivors with doulas compared with ICU survivors without doulas (OR: 0.451, 95% CI: 0.308–0.661, P < 0.001) and no ICU patients without doulas (OR: 0.244, 95% CI: 0.172–0.347, P < 0.001) (Table 3).
Table 3
Univariate and multivariate multinomial logistic regression analysis to find effects of variables on WHO-5 index
Variables
|
Univariate
|
Multivariate
|
High scores
|
Moderate scores
|
High scores
|
Moderate scores
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
WHO-5 Index
|
|
|
|
|
|
|
|
|
Age ≥ 66
|
0.953
(0.72–1.25)
|
0.727
|
0.947
(0.72–1.24)
|
0.697
|
-
|
-
|
-
|
-
|
Sex
(male vs. female)
|
0.847
(0.64–1.11)
|
0.235
|
0.808
(0.61–1.06)
|
0.128
|
-
|
-
|
-
|
-
|
Comorbidities
(yes vs. no)
|
0.622
(0.44–0.87)
|
0.007*
|
0.713
(0.50–1.01)
|
0.059
|
0.661
(0.46–0.94)
|
0.024*
|
0.712
(0.49–1.01)
|
0.062
|
APACHE II ≥ 15
|
0.882
(0.66–1.17)
|
0.392
|
0.787
(0.59–1.03)
|
0.335
|
-
|
-
|
-
|
-
|
Physical restraint
(yes vs. no)
|
0.727
(0.55–0.95)
|
0.024*
|
0.830
(0.63–1.09)
|
0.185
|
0.751
(0.56-1.00)
|
0.054
|
0.821
(0.62–1.08)
|
0.166
|
Sedation
(high vs. low/moderate)
|
0.457
(0.32–0.65)
|
< 0.001*
|
0.843
(0.57–1.24)
|
0.389
|
0.513
(0.35–0.74)
|
< 0.001*
|
0.515
(0.36–0.73)
|
< 0.001*
|
Nursing care
(high vs. low/moderate)
|
0.623
(0.44–0.87)
|
0.006*
|
0.852
(0.6–1.21)
|
0.374
|
0.701
(0.491–0.99)
|
0.05*
|
0.664
(0.47–0.92)
|
0.071*
|
Sleep disturbance
(yes vs. no)
|
0.91
(0.67–1.22)
|
0.529
|
1.326
(0.99–1.76)
|
0.054
|
-
|
-
|
-
|
-
|
Alarm silence strategy
(high vs. low and moderate)
|
0.882
(0.58–1.33)
|
0.553
|
0.983
(0.64–1.49)
|
0.935
|
-
|
-
|
-
|
-
|
Groups
(ICU survivors with doulas vs. ICU survivors without doulas)
|
4.092
(2.88–5.81)
|
< 0.001*
|
0.448
(0.30–0.67)
|
< 0.001*
|
3.614
(2.65–4.91)
|
< 0.001*
|
0.451
(0.308–0.661)
|
< 0.001*
|
Groups
(ICU survivors with doulas vs. no ICU patients without doulas)
|
3.872
(2.72–5.49)
|
< 0.001*
|
0.560
(0.37–0.84)
|
0.005*
|
2.231
(1.49–3.33)
|
< 0.001*
|
0.244
(0.172–0.347)
|
0.001*
|
Groups
(ICU survivors without doulas vs. no ICU patients without doulas)
|
0.946
(0.66–1.35)
|
0.760
|
1.249
(0.92–1.68)
|
0.149
|
-
|
-
|
-
|
-
|
* P-value < 0.05 considered as significant, for the WHO-5 index; high scores from 65 to 100 (represents of the best well-being) and moderate scores (44–64) compared to low scores between 20–43 (represents of distress) as a reference group, |
In terms of total IER-S score, multivariate analysis showed that the risk of partial PSTD over the high risk of PSTD was significantly increased in the ICU survivors with doulas compared with ICU survivors without doulas (OR: 2.650, 95% CI: 1.45–4.81, P < 0.001). In addition, the risk of moderate IER-S score was significantly increased in the ICU survivors with doulas compared with ICU survivors without doulas (OR: 2.703, 95% CI: 1.961–3.72, P < 0.001) and no ICU patients without doulas (OR: 3.103, 95% CI: 2.07–4.64, P < 0.001) (Table 4).
Table 4
Univariate and multivariate multinomial logistic regression analysis to find effects of variables on total IER-S score
Variable
|
Univariate
|
Multivariate
|
Low scores
|
Moderate scores
|
Low scores
|
Moderate scores
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
Total IES-R score
|
|
|
|
|
|
|
|
|
Age ≥ 66
|
0.981
(0.54–1.77)
|
0.949
|
1.07
(0.78–1.46)
|
0.669
|
-
|
-
|
-
|
-
|
Sex
(male vs. female)
|
0.986
(0.54–1.78)
|
0.962
|
0.888
(0.64–1.21)
|
0.460
|
-
|
-
|
-
|
-
|
Comorbidities
(yes vs. no)
|
0.937
(0.45–1.9)
|
0.860
|
1.028
(0.69–1.51)
|
0.890
|
-
|
-
|
-
|
-
|
APACHE II > 15
|
0.845
(0.45–1.5)
|
0.591
|
1.198
(0.85–1.68)
|
0.298
|
-
|
-
|
-
|
-
|
Physical restraint
(yes vs. no)
|
0.828
(0.45–1.51)
|
0.540
|
0.730
(0.52-1.00)
|
0.056
|
-
|
-
|
-
|
-
|
Sedation
(high vs. low and moderate)
|
0.291
(0.08–0.94)
|
0.040*
|
0.76
(0.49–1.16)
|
0.206
|
0.319
(0.09–1.04)
|
0.059
|
0.836
(0.54–1.28)
|
0.418
|
Nursing care
(high vs. low and moderate)
|
0.828
(0.41–1.65)
|
0.595
|
0.963
(0.65–1.41)
|
0.848
|
-
|
-
|
-
|
-
|
Sleep disturbance
(yes vs. no)
|
0.720
(0.39–1.31)
|
0.285
|
1.136
(0.81–1.59)
|
0.458
|
-
|
-
|
-
|
-
|
Alarm silence strategy
(high vs. low and moderate)
|
0.758
(0.33–1.73)
|
0.511
|
0.848
(0.53–1.33)
|
0.478
|
-
|
-
|
-
|
-
|
Groups
(ICU survivors with doulas vs. ICU survivors without doulas)
|
4.391
(1.92–10.01)
|
< 0.001*
|
2.493
(1.72–3.60)
|
< 0.001*
|
2.650
(1.45–4.81)
|
0.001*
|
2.703
(1.961–3.72)
|
< 0.001*
|
Groups
(ICU survivors with doulas vs. no ICU patients without doulas)
|
1.985
(1.02–3.85)
|
0.043*
|
3.059
(2.05–4.56)
|
< 0.001*
|
1.866
(0.95–3.63)
|
0.067
|
3.103
(2.07–4.64)
|
< 0.001*
|
Groups
(ICU survivors without doulas vs. no ICU patients without doulas)
|
0.425
(0.19–1.07)
|
0.071
|
1.227
(0.803–1.87)
|
0.344
|
-
|
-
|
-
|
-
|
* P-value < 0.05 considered as significant, for the total IES-R; low scores between 24–32 (represents of partial PSTD) and moderate scores between 33–38 (represents of PTSD) compared to high scores ≥ 39 (represents of high PTSD) as a reference group |
Multivariate analysis showed that the risk of lower score of HADS-depression over higher score of HADS-depression (normality vs. depression disorders) were significantly decreased in patients with comorbidities (OR: 0.676, 95% CI: 0.47–0.96, P = 0.029), sleep disturbance (OR: 0.591, 95% CI: 0.44–0.79, P < 0.001) and high alarm silence strategy (OR: 0.466, 95% CI: 0.31–0.69, P < 0.001). However, the risk of lower and moderate score of HADS-depression over higher score of HADS-depression (normality and suspected to depression disorders vs. definite diagnosis of depression disorders) were significantly increased in the ICU survivors with doulas compared with ICU survivors without doulas (OR: 2.731, 95% CI: 1.91–3.89, P < 0.001) and (OR: 3.722, 95% CI: 2.52–5.48, P < 0.001), and no ICU patients without doulas (OR: 2.512, 95% CI: 1.69–3.71, P < 0.001) and (OR: 3.570, 95% CI: 2.30–5.54, P < 0.001), respectively (Table 5)
Table 5
Univariate and multivariate multinomial logistic regression analysis to find effects of variables on HADS-depression
Variable
|
Univariate
|
Multivariate
|
Low scores
|
Moderate scores
|
Low scores
|
Moderate scores
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
HADS-depression score
|
|
|
|
|
|
|
|
|
Age ≥ 66
|
1.061
(0.81–1.38)
|
0.667
|
1.152
(0.84–1.57)
|
0.377
|
-
|
-
|
-
|
-
|
Sex
(male vs. female)
|
0.765
(0.58-1.00)
|
0.051
|
0.822
(0.60–1.12)
|
0.222
|
-
|
-
|
-
|
-
|
Comorbidities
(yes vs. no)
|
0.695
(0.49–0.97)
|
0.036*
|
0.960
(0.63–1.44)
|
0.843
|
0.676
(0.47–0.96)
|
0.029*
|
1.002
(0.65–1.52)
|
0.994
|
APACHE II > 15
|
0.985
(0.76–1.34)
|
0.934
|
0.878
(0.63–1.22)
|
0.443
|
-
|
-
|
-
|
-
|
Physical restraint
(yes vs. no)
|
0.876
(0.66–1.15)
|
0.342
|
0.916
(0.66–1.26)
|
0.590
|
-
|
-
|
-
|
-
|
Sedation
(high vs. low and moderate)
|
0.937
(0.66–1.33)
|
0.718
|
0.860
(0.57–1.27)
|
0.455
|
-
|
-
|
-
|
-
|
Nursing care
(high vs. low and moderate)
|
0.907
(0.65–1.25)
|
0.552
|
0.709
(0.47–1.05)
|
0.088
|
-
|
-
|
-
|
-
|
Sleep disturbance
(yes vs. no)
|
0.579
(0.43–0.76)
|
< 0.001*
|
0.960
(0.69–1.32)
|
0.803
|
0.591
(0.44–0.79)
|
< 0.001*
|
1.029
(0.73–1.43)
|
0.864
|
Alarm silence strategy
(high vs. low and moderate)
|
0.447
(0.30–0.66)
|
< 0.001*
|
0.66
(0.42–1.02)
|
0.065
|
0.466
(0.31–0.69)
|
< 0.001*
|
0.693
(0.44–1.08)
|
0.111
|
Groups
(ICU survivors with doulas vs. ICU survivors without doulas)
|
3.216
(2.20–4.70)
|
< 0.001*
|
3.992
(2.61–6.10)
|
< 0.001*
|
2.731
(1.91–3.89)
|
< 0.001*
|
3.722
(2.52–5.48)
|
< 0.001*
|
Groups
(ICU survivors with doulas vs. no ICU patients without doulas)
|
2.581
(1.75–3.80)
|
< 0.001*
|
3.485
(2.25–5.39)
|
< 0.001*
|
2.512
(1.69–3.71)
|
< 0.001*
|
3.570
(2.30–5.54)
|
< 0.001*
|
Groups
(ICU survivors without doulas vs. no ICU patients without doulas)
|
0.803
(0.592–1.08)
|
0.157
|
0.873
(0.601–1.26)
|
0.475
|
-
|
-
|
-
|
-
|
* P-value < 0.05 considered as significant, for the HADS depression score; low scores ≤ 7 (no depression) and moderate scores between 8–10 (suspected of depression) compared to high scores ≥ 11 (definitive diagnosis of depression) as a reference group |
Table 6
Univariate and multivariate multinomial logistic regression analysis to find effects of variables on HADS- anxiety score
Variable
|
Univariate
|
Multivariate
|
Low scores
|
Moderate scores
|
Low scores
|
Moderate scores
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
OR
(95% CI)
|
P-value
|
HADS-anxiety score
|
|
|
|
|
|
|
|
|
Age ≥ 66
|
0.917
(0.71–1.18)
|
0.503
|
0.853
(0.63–1.15)
|
0.298
|
-
|
-
|
-
|
-
|
Sex
(male vs. female)
|
0.932
(0.72–1.20)
|
0.590
|
0.741
(0.54-1.00)
|
0.053
|
-
|
-
|
-
|
-
|
Comorbidities
(yes vs. no)
|
0.666
(0.48–0.91)
|
0.010*
|
0.799
(0.55–1.16)
|
0.239
|
0.695
(0.50–0.95)
|
0.025*
|
0.810
(0.55–1.17)
|
0.270
|
APACHE II > 15
|
0.395
(0.13–1.44)
|
0.512
|
0.863
(0.63–1.18)
|
0.359
|
-
|
-
|
-
|
-
|
Physical restraint
(yes vs. no)
|
0.872
(0.67–1.12)
|
0.297
|
0.883
(0.65–1.19)
|
0.424
|
-
|
-
|
-
|
-
|
Sedation
(high vs. low and moderate)
|
0.966
(0.69–1.34)
|
0.839
|
0.906
(0.62–1.32)
|
0.611
|
-
|
-
|
-
|
-
|
Nursing care
(high vs. low and moderate)
|
0.847
(0.62–1.15)
|
0.291
|
0.993
(0.68–1.44)
|
0.970
|
-
|
-
|
-
|
-
|
Sleep disturbance
(yes vs. no)
|
0.969
(0.73–1.27)
|
0.822
|
0.883
(0.64–1.20)
|
0.437
|
-
|
-
|
-
|
-
|
Alarm silence strategy
(high vs. low and moderate)
|
0.864
(0.58–1.27)
|
0.459
|
0.791
(0.49–1.26)
|
0.328
|
-
|
-
|
-
|
-
|
Groups
(ICU survivors with doulas vs. ICU survivors without doulas)
|
3.093
(2.24–4.25)
|
< 0.001*
|
1.556
(1.06–2.28)
|
0.024*
|
3.050
(2.30–4.04)
|
< 0.001*
|
1.504
(1.06–2.12)
|
0.021*
|
Groups
(ICU survivors with doulas vs. no ICU patients without doulas)
|
3.073
(2.21–4.26)
|
< 0.001*
|
1.468
(0.99–2.16)
|
0.054
|
3.044
(2.19–4.23)
|
< 0.001*
|
1.461
(0.98–2.16)
|
0.057
|
Groups
(ICU survivors without doulas vs. no ICU patients without doulas)
|
0.993
(0.722–1.36)
|
0.968
|
0.944
(0.667–1.33)
|
0.744
|
-
|
-
|
-
|
-
|
* P-value < 0.05 considered as significant, for the HADS anxiety score; low scores ≤ 7 (no anxiety) and moderate scores between 8–10 (suspected of anxiety disorders) compared to high scores ≥ 11 (definitive diagnosis of anxiety) as a reference group |
In terms of HADS-anxiety score, multivariate analysis showed that the having comorbidities was associated with a decreased the risk of normality vs. anxiety disorders (OR: 0.695, 95% CI: 0.50–0.95, P = 0.025). However, in the ICU survivors with doulas the risk of normality vs. definite diagnosis of anxiety was increased compared with ICU survivors without doulas (OR: 3.050, 95% CI: 2.30–4.04, P < 0.001) and no ICU patients without doulas (OR: 3.044, 95% CI: 2.19–4.23, P < 0.001). In addition, the risk of suspected of anxiety disorders vs. definitive diagnosis of anxiety was significantly increased in the in the ICU survivors with doulas compared with ICU survivors without doulas (OR: 1.504, 95% CI: 1.06–2.12, P = 0.021).
Binary logistic regression finding
Univariate and multivariate binary logistic regression analysis to find effects of variables on outcomes (WHO-5 index, HADS-depression/anxiety scores and IES-R score) are presented in Figs. 1 to 4. In terms of WHO-5 index, multivariate analysis showed that the required more nursing care was associated with an increased the risk of mental distress (OR: 1.455, 95%CI: 1.089–1.944, P = 0.011). However, in the ICU survivor with doulas the risk of mental distress was significantly decreased compared with the ICU survivors without doulas (OR: 0.344, 95% CI: 0.265–0.446, P < 0.001) and no ICU patients without doulas (OR: 0.332, 95% CI: 0.246–0.448, P < 0.001) (Fig. 1B). As we can see in the Fig. 2B related to the HADS-depression score, comorbidities (OR: 1.481, 95% CI: 1.115–1.967, P = 0.007), sleep disturbance (OR: 1.721, 95% CI: 1.349–2.196, P < 0.001), and alarm silence strategy (OR: 1.821, 95% CI: 1.276-2.6, P < 0.001) were associated with an increased the risk of depression disorders. However, the risk of depression disorders was decreased in the ICU survivor with doulas than the ICU survivors without doulas (OR: 0.762, 95% CI: 0.593–0.979, P = 0.033). In terms of HADS-anxiety, multivariate analysis showed that the comorbidities were associated with an increased the risk of anxiety disorders (OR: 1.348, 95%CI: 1.004–1.812, P = 0.047). However, the risk of anxiety disorders was decreased in the ICU survivor with doulas than the ICU survivors without doulas (OR: 0.372, 95% CI: 0.287–0.481, P < 0.001) and no ICU patients without doulas (OR: 0.370, 95% CI: 0.273–0.501, P < 0.001) (Fig. 3B). According to the results of multivariate analysis related to the IES-R score, the OR of PTSD significantly decreased in the ICU survivor with doulas than the ICU survivors without doulas (OR: 0.367, 95%CI: 0.184–0.875, P < 0.001) (Fig. 4B).