Follow up manuscript FINAL FEB 25

Author

KS, MT

Published

March 29, 2025

Metadata and analysis of the file: CPXE_26.2.23.follow-up stary + nowy.xlsx
Outcomes are derived from file: 25.2.24.CPXE_3_26_follow_up dla MT_12.03_DATA_ANKIETY.xlsx

Each row represents one patient.

Outcome:
combined_end_point_card_neo_death

Beginning of follow up:
First MI date, “mi_date_hosp” variable

Number of observations:
100

Number of variables:
676

RER threshold:
0

RER variable:
rer_maxload_cpxe_1

EF threshold:
40

EF variable:
ef_rest_prec_cpxe_1

Excluding patient with ID 7 (hyperthyroid):
TRUE

Outcome

# A tibble: 3 × 2
  combined_end_point_card_neo_death     n
                              <dbl> <int>
1                                 0    74
2                                 1    16
3                                NA    10

Excluding patient with hyperthyroidism

 [1] 1.14 1.44 1.05 1.07 1.22 1.22 1.20 1.23 1.33 1.21 1.06 1.05 1.20 1.22 1.31
[16] 1.15 1.22 1.07 1.16 1.05 1.08 1.07   NA 1.11 1.11 1.12 1.28 1.25 1.00 1.12
[31] 1.22 1.06 1.12 1.14 1.05 0.92 0.93 1.20 1.20 0.92 1.05 1.23 1.26 1.07 1.12
[46] 1.05 0.99 1.08 0.98 1.21 1.05 1.08 1.11 1.05 1.00 1.12 1.20 0.83 1.05 1.16
[61] 1.05 0.97 1.25 1.15 1.05 1.05 1.05 1.10 1.05 1.18 0.99 0.95 1.17 0.96 1.12
[76] 1.18 1.12 1.15 1.19 0.99 1.07 1.21 1.37 1.39 1.30 1.05 1.30 1.30 0.94 1.30
[91] 1.11 1.32 1.12 1.20 1.15 1.05 1.05 1.09 1.14
 [1] 1.14 1.44 1.05 1.07 1.22 1.22 1.20 1.23 1.33 1.21 1.06 1.05 1.20 1.22 1.31
[16] 1.15 1.22 1.07 1.16 1.05 1.08 1.07   NA 1.11 1.11 1.12 1.28 1.25 1.00 1.12
[31] 1.22 1.06 1.12 1.14 1.05 0.92 0.93 1.20 1.20 0.92 1.05 1.23 1.26 1.07 1.12
[46] 1.05 1.00 1.08 1.00 1.21 1.05 1.08 1.11 1.05 1.00 1.12 1.20 0.83 1.05 1.16
[61] 1.05 1.00 1.25 1.15 1.05 1.05 1.05 1.10 1.05 1.18 1.00 1.00 1.17 1.00 1.12
[76] 1.18 1.12 1.15 1.19 1.00 1.07 1.21 1.37 1.39 1.30 1.05 1.30 1.30 0.94 1.30
[91] 1.11 1.32 1.12 1.20 1.15 1.05 1.05 1.09 1.14

Number of observations after removal of patients with hyperthyroid:
99

All our RER initialy

 [1] 1.14 1.44 1.05 1.07 1.22 1.22 1.20 1.23 1.33 1.21 1.06 1.05 1.20 1.22 1.31
[16] 1.15 1.22 1.07 1.16 1.05 1.08 1.07   NA 1.11 1.11 1.12 1.28 1.25 1.00 1.12
[31] 1.22 1.06 1.12 1.14 1.05 0.92 0.93 1.20 1.20 0.92 1.05 1.23 1.26 1.07 1.12
[46] 1.05 1.00 1.08 1.00 1.21 1.05 1.08 1.11 1.05 1.00 1.12 1.20 0.83 1.05 1.16
[61] 1.05 1.00 1.25 1.15 1.05 1.05 1.05 1.10 1.05 1.18 1.00 1.00 1.17 1.00 1.12
[76] 1.18 1.12 1.15 1.19 1.00 1.07 1.21 1.37 1.39 1.30 1.05 1.30 1.30 0.94 1.30
[91] 1.11 1.32 1.12 1.20 1.15 1.05 1.05 1.09 1.14

RER below threshold patients

# A tibble: 0 × 3
# ℹ 3 variables: id <dbl>, rer_maxload_cpxe_1 <dbl>,
#   combined_end_point_card_neo_death <dbl>

All out EF initially

 [1] 54 38 59 56 35 51 52 55 48 68 50 74 50 46 59 58 60 47 57 64 57 48 NA 57 53
[26] 69 58 58 47 41 54 48 51 70 54 56 46 55 61 51 59 71 57 56 54 55 57 49 68 48
[51] 45 57 67 66 60 56 59 61 65 56 57 64 59 59 53 50 48 67 51 51 61 59 66 42 53
[76] 32 47 27 59 73 55 55 54 64 79 62 60 53 42 73 58 61 66 59 62 68 66 63 73

EF below threshold

# A tibble: 4 × 3
     id ef_rest_prec_cpxe_1 combined_end_point_card_neo_death
  <dbl>               <dbl>                             <dbl>
1     2                  38                                NA
2     5                  35                                 0
3    77                  32                                 0
4    79                  27                                 0

All patients removed in this analysis (RER and EF below threshold)

# A tibble: 4 × 4
     id ef_rest_prec_cpxe_1 rer_maxload_cpxe_1 combined_end_point_card_neo_death
  <dbl>               <dbl>              <dbl>                             <dbl>
1     2                  38               1.44                                NA
2     5                  35               1.22                                 0
3    77                  32               1.18                                 0
4    79                  27               1.15                                 0

Number of observations after removal of below EF and RER threshold + removal of hyperthyroid
95

Patients with missing outcome

# A tibble: 9 × 2
     id combined_end_point_card_neo_death
  <dbl>                             <dbl>
1    24                                NA
2    25                                NA
3    28                                NA
4    29                                NA
5    39                                NA
6    42                                NA
7    58                                NA
8    81                                NA
9    99                                NA

Number of observations after removal of patients with missing outcome + low RER (if any) + low EF + hyperthyroid:
86

ROC curve

  PANEL group     AUC
1     1    -1 0.72723
  threshold sensitivity specificity
1  66.86627   0.8666667   0.5492958
Sensitivity: 0.8666667 
Specificity: 0.5492958 
PPV: 0.2888889 
NPV: 0.9512195 

Survival analysis

Survival analysis table

id como_end_point_card_neo_death days_difference mi_date_hosp date_endpoint_or_survey
19 1 104 2016-05-08 2016-08-20
3 1 131 2015-11-14 2016-03-24
10 1 144 2016-03-18 2016-08-09
91 0 193 2018-07-05 2019-01-14
92 1 203 2018-06-04 2018-12-24
87 0 220 2018-05-28 2019-01-03
85 1 229 2018-04-16 2018-12-01
89 0 231 2018-05-18 2019-01-04
62 1 270 2017-05-07 2018-02-01
23 1 346 2016-06-10 2017-05-22
76 0 368 2017-12-31 2019-01-03
71 1 393 2017-05-31 2018-06-28
70 1 471 2017-09-30 2019-01-14
48 1 531 2017-01-04 2018-06-19
65 0 603 2017-05-09 2019-01-02
59 0 626 2017-04-16 2019-01-02
61 0 628 2017-04-14 2019-01-02
56 0 636 2017-04-06 2019-01-02
54 0 644 2017-03-29 2019-01-02
49 0 647 2017-03-22 2018-12-29
95 0 668 2018-09-30 2020-07-29
47 0 729 2016-12-30 2018-12-29
46 0 739 2017-01-05 2019-01-14
94 0 772 2018-11-05 2020-12-16
15 0 773 2016-04-19 2018-06-01
37 0 809 2016-10-17 2019-01-04
50 1 814 2017-04-03 2019-06-26
31 0 849 2016-08-31 2018-12-28
30 0 852 2016-08-27 2018-12-27
27 0 905 2016-07-05 2018-12-27
14 0 932 2016-04-11 2018-10-30
21 0 959 2016-05-20 2019-01-04
8 1 1020 2016-03-17 2019-01-01
11 0 1047 2016-03-31 2019-02-11
6 0 1053 2015-12-12 2018-10-30
67 0 1097 2017-07-28 2020-07-29
80 0 1187 2017-12-14 2021-03-15
43 0 1343 2016-11-24 2020-07-29
38 0 1381 2016-10-17 2020-07-29
100 0 1422 2018-07-24 2022-06-15
96 0 1433 2018-11-07 2022-10-10
98 0 1446 2018-10-25 2022-10-10
74 0 1450 2018-11-03 2022-10-23
97 0 1450 2018-10-21 2022-10-10
4 1 1458 2016-01-04 2020-01-01
90 0 1465 2018-05-28 2022-06-01
88 0 1481 2018-05-12 2022-06-01
93 0 1520 2018-04-17 2022-06-15
86 0 1641 2018-04-26 2022-10-23
75 0 1651 2017-11-23 2022-06-01
83 0 1684 2017-10-21 2022-06-01
84 0 1698 2018-02-28 2022-10-23
78 0 1746 2017-12-29 2022-10-10
73 0 1827 2017-10-09 2022-10-10
82 0 1832 2017-10-04 2022-10-10
72 0 1844 2017-10-05 2022-10-23
69 0 1859 2017-09-07 2022-10-10
68 0 1905 2017-08-05 2022-10-23
66 0 1984 2017-05-05 2022-10-10
55 0 2014 2017-04-05 2022-10-10
64 0 2015 2017-05-24 2022-11-29
57 0 2019 2017-03-31 2022-10-10
63 0 2024 2017-05-15 2022-11-29
36 0 2038 2016-10-01 2022-05-01
51 0 2043 2017-03-20 2022-10-23
35 0 2047 2016-09-22 2022-05-01
60 0 2050 2017-04-21 2022-12-01
52 0 2060 2017-03-03 2022-10-23
53 0 2063 2017-02-15 2022-10-10
17 1 2097 2016-05-06 2022-02-01
1 0 2118 2015-08-23 2021-06-10
44 0 2130 2016-12-10 2022-10-10
22 0 2151 2016-06-10 2022-05-01
41 0 2163 2016-11-07 2022-10-10
40 0 2178 2016-10-23 2022-10-10
18 1 2219 2016-05-04 2022-06-01
33 0 2223 2016-09-08 2022-10-10
32 0 2228 2016-09-03 2022-10-10
34 0 2230 2016-09-14 2022-10-23
26 0 2288 2016-07-05 2022-10-10
45 0 2331 2016-05-23 2022-10-10
20 0 2340 2016-05-14 2022-10-10
12 0 2383 2016-04-01 2022-10-10
13 0 2383 2016-04-01 2022-10-10
9 0 2395 2016-03-20 2022-10-10
16 0 2415 2016-04-19 2022-11-29

Models

Logistic regression model

Characteristic

OR

1

95% CI

1

p-value

hr_proc_pred_peak_cpxe_1 0.89 0.79, 0.98 0.034
lvedv_index_rest_ml_m2_cpxe_1 1.06 0.99, 1.13 0.091
sex_hosp


    0
    1 0.84 0.12, 6.05 0.9
age_mi_hosp 0.99 0.91, 1.08 0.8
physical_activity_before_mi_small_hosp


    0
    1 4.34 0.83, 24.4 0.081
smoking_past_hosp


    0
    1 5.27 1.06, 27.8 0.042
bisoprolol_equ_dose 1.03 0.78, 1.34 0.8
creatinine_discharge_hosp 0.06 0.00, 17.2 0.4
1

OR = Odds Ratio, CI = Confidence Interval

COX regression model

Characteristic

HR

1

95% CI

1

p-value

hr_proc_pred_peak_cpxe_1 0.87 0.80, 0.95 0.003
lvedv_index_rest_ml_m2_cpxe_1 1.04 0.98, 1.09 0.2
sex_hosp


    0
    1 1.28 0.24, 6.89 0.8
age_mi_hosp 0.98 0.91, 1.04 0.5
physical_activity_before_mi_small_hosp


    0
    1 2.43 0.72, 8.22 0.2
smoking_past_hosp


    0
    1 5.92 1.59, 22.0 0.008
bisoprolol_equ_dose 1.00 0.81, 1.22 >0.9
creatinine_discharge_hosp 0.06 0.00, 9.53 0.3
1

HR = Hazard Ratio, CI = Confidence Interval

COX model parameters and VIF

Call:
coxph(formula = Surv(days_difference, outcome) ~ ., data = data_to_cox_regression_model)

                                            coef exp(coef) se(coef)      z
hr_proc_pred_peak_cpxe_1                -0.13682   0.87213  0.04603 -2.972
lvedv_index_rest_ml_m2_cpxe_1            0.03489   1.03550  0.02711  1.287
sex_hosp1                                0.24940   1.28326  0.85713  0.291
age_mi_hosp                             -0.02464   0.97566  0.03352 -0.735
physical_activity_before_mi_small_hosp1  0.88852   2.43153  0.62172  1.429
smoking_past_hosp1                       1.77812   5.91870  0.66995  2.654
bisoprolol_equ_dose                     -0.00377   0.99624  0.10402 -0.036
creatinine_discharge_hosp               -2.80835   0.06030  2.58300 -1.087
                                              p
hr_proc_pred_peak_cpxe_1                0.00296
lvedv_index_rest_ml_m2_cpxe_1           0.19814
sex_hosp1                               0.77107
age_mi_hosp                             0.46231
physical_activity_before_mi_small_hosp1 0.15296
smoking_past_hosp1                      0.00795
bisoprolol_equ_dose                     0.97109
creatinine_discharge_hosp               0.27693

Likelihood ratio test=29.39  on 8 df, p=0.0002707
n= 86, number of events= 15 
              hr_proc_pred_peak_cpxe_1          lvedv_index_rest_ml_m2_cpxe_1 
                              1.311033                               2.150919 
                              sex_hosp                            age_mi_hosp 
                              2.343571                               1.659278 
physical_activity_before_mi_small_hosp                      smoking_past_hosp 
                              1.368142                               1.415072 
                   bisoprolol_equ_dose              creatinine_discharge_hosp 
                              1.165885                               1.982311 

Table

# A tibble: 5 × 4
# Groups:   combined_end_point_card_neo_death [2]
  combined_end_point_card_neo_death rer_maxload_group count percentage
                              <dbl> <chr>             <int>      <dbl>
1                                 0 high                 60      84.5 
2                                 0 low                   5       7.04
3                                 0 med                   6       8.45
4                                 1 high                 13      86.7 
5                                 1 med                   2      13.3 

Characteristic

N

Overall, N = 65

1

0, N = 52

1

1, N = 13

1

p-value

2
rehab_reh 65 65 (100%) 52 (100%) 13 (100%)
cr_k 65 45 (69%) 35 (67%) 10 (77%) 0.7
cr_ik 65 20 (31%) 17 (33%) 3 (23%) 0.7
cr 65 65 (100%) 52 (100%) 13 (100%)
karta_informacyjna_1_jest_0_brak 65 54 (83%) 44 (85%) 10 (77%) 0.7
model_rehab_poczatkowy 65


0.3
    A
5 (7.7%) 5 (9.6%) 0 (0%)
    B
41 (63%) 34 (65%) 7 (54%)
    C
19 (29%) 13 (25%) 6 (46%)
model_rehab_koncowy 65


0.6
    A
15 (23%) 13 (25%) 2 (15%)
    B
44 (68%) 35 (67%) 9 (69%)
    C
6 (9.2%) 4 (7.7%) 2 (15%)
lliczba_treningow_wykonanych 65 21.00 (21.00, 24.00) 21.00 (21.00, 24.00) 21.00 (21.00, 21.00) 0.8
training_sessions_planned 65 21.00 (21.00, 24.00) 21.00 (21.00, 24.00) 21.00 (21.00, 21.00) 0.5
c_training_hr_limit_40_50_perc 6 6 (100%) 4 (100%) 2 (100%)
    Unknown
59 48 11
b_training_hr_limit_50_60_perc 44 44 (100%) 35 (100%) 9 (100%)
    Unknown
21 17 4
a_training_hr_limit_60_80_perc 15 15 (100%) 13 (100%) 2 (100%)
    Unknown
50 39 11
1

n (%); Median (IQR)

2

Fisher’s exact test; Wilcoxon rank sum test

Mean and SD for group with outcome

[1] 1.129333
[1] 0.1157872

Mean and SD for group without outcome

[1] 1.118732
[1] 0.1079805

Mean and SD for whole

[1] 1.120581
[1] 0.1087507