I decided to pick composite of CV death (not all cause death) and hospitalization, because pt ID 1 has good increase in SV and cancer death and this could interfere with our analysis. Here 2 pts with missing outcome are removed.
outcome_cvdeath | sv_0 | sv_1 | sv_2 | change_baseline_to_first | change_baseline_to_second | change_first_to_second |
---|---|---|---|---|---|---|
0 | 56.4 | 66.4 | 71.4 | 20.2 | 38.2 | 9.8 |
1 | 48.6 | 57.1 | 44.5 | 13.2 | 2.8 | -3.9 |
outcome_composite | sv_0 | sv_1 | sv_2 | change_baseline_to_first | change_baseline_to_second | change_first_to_second |
---|---|---|---|---|---|---|
0 | 56.8 | 70.1 | 76.9 | 25.6 | 46.5 | 10.7 |
1 | 53.8 | 58.9 | 54.8 | 11.4 | 12.3 | 5.4 |
Those in red box have outcome of CV death or hospitalization
id | sv_0 | sv_1 | sv_2 | change_baseline_to_first | change_baseline_to_second | change_first_to_second | mitra_clip_date | outcome_composite | outcome_cvdeath | outcome_date |
---|---|---|---|---|---|---|---|---|---|---|
1 | 80.7 | 97.1 | 67.9 | 20.3 | -15.9 | -30.1 | 2016-10-28 | 0 | 0 | 2018-12-14 |
3 | 54.0 | 80.0 | 114.9 | 48.1 | 112.8 | 43.6 | 2016-12-02 | 0 | 0 | 2019-08-06 |
4 | 75.1 | 71.4 | NA | -4.9 | NA | NA | 2016-12-16 | 1 | 0 | 2018-05-26 |
5 | 32.8 | 57.9 | 58.6 | 76.5 | 78.7 | 1.2 | 2017-01-20 | 1 | 0 | 2018-04-20 |
6 | 58.4 | 54.4 | 60.0 | -6.8 | 2.7 | 10.3 | 2017-02-17 | 1 | 1 | 2018-09-14 |
7 | 49.9 | 76.2 | 95.2 | 52.7 | 90.8 | 24.9 | 2017-03-17 | 0 | 0 | 2019-10-28 |
8 | 74.1 | 62.9 | 77.3 | -15.1 | 4.3 | 22.9 | 2017-04-21 | 1 | 0 | 2017-12-16 |
9 | 47.7 | 58.2 | 56.6 | 22.0 | 18.7 | -2.7 | 2017-05-05 | 0 | 0 | 2019-07-18 |
11 | 70.9 | 74.3 | 72.5 | 4.8 | 2.3 | -2.4 | 2017-06-09 | 0 | 0 | 2019-06-07 |
13 | 35.9 | NA | 47.5 | NA | 32.3 | NA | 2017-06-23 | 1 | 1 | 2017-12-24 |
14 | 48.1 | 69.7 | 106.5 | 44.9 | 121.4 | 52.8 | 2017-08-18 | 0 | 0 | 2019-08-19 |
16 | 34.3 | 30.8 | 63.7 | -10.2 | 85.7 | 106.8 | 2017-09-20 | 1 | 0 | 2018-04-20 |
17 | 68.7 | 60.5 | NA | -11.9 | NA | NA | 2017-09-20 | 1 | 1 | 2017-09-27 |
18 | 71.1 | 82.9 | 77.6 | 16.6 | 9.1 | -6.4 | 2017-10-18 | 0 | 0 | 2019-10-21 |
20 | 54.7 | 49.8 | NA | -9.0 | NA | NA | 2017-11-15 | 1 | 0 | 2020-02-12 |
21 | 80.1 | 64.7 | 66.9 | -19.2 | -16.5 | 3.4 | 2017-11-15 | 0 | 0 | 2018-01-26 |
22 | 49.4 | NA | 73.7 | NA | 49.2 | NA | 2017-12-20 | 0 | 0 | 2018-08-21 |
23 | 75.8 | 96.8 | 59.2 | 27.7 | -21.9 | -38.8 | 2018-01-17 | 0 | 0 | 2020-02-19 |
24 | 58.0 | 70.8 | 80.7 | 22.1 | 39.1 | 14.0 | 2018-02-21 | 0 | 0 | 2019-11-26 |
25 | 66.4 | 81.8 | 66.2 | 23.2 | -0.3 | -19.1 | 2018-03-21 | 1 | 0 | 2018-04-13 |
26 | 48.1 | 48.5 | 46.8 | 0.8 | -2.7 | -3.5 | 2018-05-16 | 1 | 0 | 2019-11-25 |
27 | 36.3 | 66.1 | 73.4 | 82.1 | 102.2 | 11.0 | 2018-05-16 | 0 | 0 | 2020-01-08 |
28 | 37.4 | 52.5 | 43.0 | 40.4 | 15.0 | -18.1 | 2018-05-30 | 1 | 1 | 2019-10-02 |
30 | 46.5 | 61.0 | NA | 31.2 | NA | NA | 2018-06-06 | 1 | 1 | 2019-09-06 |
31 | 63.2 | 66.0 | 59.4 | 4.4 | -6.0 | -10.0 | 2018-06-18 | 1 | 0 | 2018-08-17 |
32 | 44.5 | NA | 27.3 | NA | -38.7 | NA | 2018-09-12 | 1 | 1 | 2019-10-28 |
33 | 47.9 | 47.7 | NA | -0.4 | NA | NA | 2018-09-12 | 0 | 0 | NA |
35 | 35.2 | 48.7 | 52.5 | 38.4 | 49.1 | 7.8 | 2018-11-21 | 0 | 0 | 2020-03-13 |
36 | 52.8 | 68.5 | 39.9 | 29.7 | -24.4 | -41.8 | 2019-02-20 | 1 | 0 | 2019-12-07 |
37 | 40.2 | 58.5 | 98.7 | 45.5 | 145.5 | 68.7 | 2019-02-20 | 0 | 0 | 2020-01-23 |
38 | 71.4 | 76.7 | NA | 7.4 | NA | NA | 2019-06-17 | 0 | 0 | 2020-07-02 |
40 | 55.6 | 53.8 | 56.1 | -3.2 | 0.9 | 4.3 | 2019-10-16 | 0 | 0 | 2020-09-07 |
41 | 67.2 | NA | 67.9 | NA | 1.0 | NA | 2020-01-29 | 1 | NA | 2021-12-01 |
42 | 95.5 | NA | 87.2 | NA | -8.7 | NA | 2020-01-29 | NA | NA | 2020-09-08 |
44 | 49.9 | NA | 78.5 | NA | 57.3 | NA | 2020-05-20 | 0 | 0 | 2021-01-18 |
threshold <- 35.7
## Call: survfit(formula = surv_object ~ above_threshold_change, data = survival_dataset)
##
## 7 observations deleted due to missingness
## n events median 0.95LCL 0.95UCL
## above_threshold_change=0 17 10 558 411 NA
## above_threshold_change=1 11 2 NA NA NA
## Call:
## coxph(formula = surv_object ~ above_threshold_change, data = survival_dataset)
##
## coef exp(coef) se(coef) z p
## above_threshold_change1 -1.2345 0.2910 0.7757 -1.592 0.111
##
## Likelihood ratio test=3.23 on 1 df, p=0.07238
## n= 28, number of events= 12
## (7 observations deleted due to missingness)
id | baseline_to_first | baseline_to_second | first_to_second | outcome_com | outcome_cvdeath | above_threshold | surv_object |
---|---|---|---|---|---|---|---|
1 | 20.3 | -15.9 | -30.1 | 0 | 0 | 0 | 777+ |
3 | 48.1 | 112.8 | 43.6 | 0 | 0 | 1 | 977+ |
4 | -4.9 | NA | NA | 1 | 0 | NA | 526 |
5 | 76.5 | 78.7 | 1.2 | 1 | 0 | 1 | 455 |
6 | -6.8 | 2.7 | 10.3 | 1 | 1 | 0 | 574 |
7 | 52.7 | 90.8 | 24.9 | 0 | 0 | 1 | 955+ |
8 | -15.1 | 4.3 | 22.9 | 1 | 0 | 0 | 239 |
9 | 22.0 | 18.7 | -2.7 | 0 | 0 | 0 | 804+ |
11 | 4.8 | 2.3 | -2.4 | 0 | 0 | 0 | 728+ |
13 | NA | 32.3 | NA | 1 | 1 | 0 | 184 |
14 | 44.9 | 121.4 | 52.8 | 0 | 0 | 1 | 731+ |
16 | -10.2 | 85.7 | 106.8 | 1 | 0 | 1 | 212 |
17 | -11.9 | NA | NA | 1 | 1 | NA | 7 |
18 | 16.6 | 9.1 | -6.4 | 0 | 0 | 0 | 733+ |
20 | -9.0 | NA | NA | 1 | 0 | NA | 819 |
21 | -19.2 | -16.5 | 3.4 | 0 | 0 | 0 | 72+ |
22 | NA | 49.2 | NA | 0 | 0 | 1 | 244+ |
23 | 27.7 | -21.9 | -38.8 | 0 | 0 | 0 | 763+ |
24 | 22.1 | 39.1 | 14.0 | 0 | 0 | 1 | 643+ |
25 | 23.2 | -0.3 | -19.1 | 1 | 0 | 0 | 23 |
26 | 0.8 | -2.7 | -3.5 | 1 | 0 | 0 | 558 |
27 | 82.1 | 102.2 | 11.0 | 0 | 0 | 1 | 602+ |
28 | 40.4 | 15.0 | -18.1 | 1 | 1 | 0 | 490 |
30 | 31.2 | NA | NA | 1 | 1 | NA | 457 |
31 | 4.4 | -6.0 | -10.0 | 1 | 0 | 0 | 60 |
32 | NA | -38.7 | NA | 1 | 1 | 0 | 411 |
33 | -0.4 | NA | NA | 0 | 0 | NA | NA+ |
35 | 38.4 | 49.1 | 7.8 | 0 | 0 | 1 | 478+ |
36 | 29.7 | -24.4 | -41.8 | 1 | 0 | 0 | 290 |
37 | 45.5 | 145.5 | 68.7 | 0 | 0 | 1 | 337+ |
38 | 7.4 | NA | NA | 0 | 0 | NA | 381+ |
40 | -3.2 | 0.9 | 4.3 | 0 | 0 | 0 | 327+ |
41 | NA | 1.0 | NA | 1 | NA | 0 | 672 |
42 | NA | -8.7 | NA | NA | NA | 0 | 223? |
44 | NA | 57.3 | NA | 0 | 0 | 1 | 243+ |