Survival Omar

Author

Dr Omar Bin Nazmi

Published

January 21, 2024

Show the Code
library(haven)
Warning: package 'haven' was built under R version 4.4.2
Show the Code
library(here)
Warning: package 'here' was built under R version 4.4.2
here() starts at C:/Users/H P/OneDrive - The Goose and Duck/Documents
Show the Code
library(survminer)
Warning: package 'survminer' was built under R version 4.4.2
Loading required package: ggplot2
Loading required package: ggpubr
Warning: package 'ggpubr' was built under R version 4.4.2
Show the Code
library(lubridate)

Attaching package: 'lubridate'
The following objects are masked from 'package:base':

    date, intersect, setdiff, union
Show the Code
library(survival)

Attaching package: 'survival'
The following object is masked from 'package:survminer':

    myeloma
Show the Code
library(readxl)
Warning: package 'readxl' was built under R version 4.4.2
Show the Code
library(gtsummary)
Warning: package 'gtsummary' was built under R version 4.4.2
Show the Code
datasurv <- read_excel("C:/Users/H P/OneDrive - The Goose and Duck/Desktop/Data Survival/Data Survival NHL.xlsx")
str(datasurv)
tibble [448 × 41] (S3: tbl_df/tbl/data.frame)
 $ Age                                             : num [1:448] 38 19 21 60 81 41 69 49 55 48 ...
 $ <60=0>60=1                                      : num [1:448] 0 0 0 1 1 0 1 0 0 0 ...
 $ 15-30, 31-59, >60                               : num [1:448] 1 0 0 2 2 1 2 1 1 1 ...
 $ AgeGroup 1544:1,4554:2,5564 :3,6574:4.>75:5)    : num [1:448] 1 1 1 3 5 1 4 2 3 2 ...
 $ Agegroupbaru 15400, 41601,61-752,753            : num [1:448] 0 0 0 1 3 1 2 1 1 1 ...
 $ Topography                                      : num [1:448] 779 779 779 779 779 779 779 779 779 779 ...
 $ Indolent0aggressive1veryaggressiive2            : num [1:448] 1 0 1 1 0 1 1 1 1 1 ...
 $ very aggressive0,agrresive1,indolent2           : num [1:448] 1 2 1 1 2 1 1 1 1 1 ...
 $ Indolent0,Aggressive1very aggressive2latest     : num [1:448] 1 0 1 1 0 1 1 1 1 1 ...
 $ Indolent :0,Agressive:1,NOS:2                   : num [1:448] 2 0 2 2 0 2 2 2 1 2 ...
 $ 0 OTHER, 1 B CELL                               : num [1:448] 1 1 1 1 1 1 1 1 0 1 ...
 $ Bcell0, othertcellnkcell1unclassifiable 2       : num [1:448] 1 0 0 0 0 0 0 0 2 0 ...
 $ Bcell1othercell0unclassifiable2                 : num [1:448] 0 1 1 1 1 1 1 1 2 1 ...
 $ B cell0,tcell1,precursor,2,unclasifiable3       : num [1:448] 1 0 0 0 0 0 0 0 3 0 ...
 $ 1=B cell,0=Tcell,2 precursor, Unclassifiable    : num [1:448] 1 1 1 1 1 1 1 1 3 1 ...
 $ Aggressive= 1, indolent = 2                     : num [1:448] 1 2 1 1 2 1 1 1 1 1 ...
 $ Morphology                                      : num [1:448] 9591 9691 9591 9591 9695 ...
 $ stage 10,stage21,stage32,satge43                : num [1:448] 2 0 4 2 4 2 4 4 2 4 ...
 $ Stage10,stage231,stage42unnkown 3               : num [1:448] 1 0 3 1 3 1 3 3 1 3 ...
 $ stage120,stage31satge42unkown3                  : num [1:448] 1 0 3 1 3 1 3 3 1 3 ...
 $ stage 0=stage 1,2, 1 stage 3,4, 2= unknown group: num [1:448] 1 0 2 1 2 1 2 2 1 2 ...
 $ Stage                                           : num [1:448] 3 1 9 3 9 3 9 9 3 9 ...
 $ ICD-10                                          : chr [1:448] "C859" "C821" "C859" "C859" ...
 $ Chemotherapyyes0mo1unkown2                      : num [1:448] 0 1 1 0 1 0 1 0 0 1 ...
 $ Chemotherapy 0= NO, Yes= 1, 2= unknown          : num [1:448] 1 0 0 1 0 1 0 1 1 0 ...
 $ Radiotherapy0yes,no1,unknown2                   : num [1:448] 1 1 1 1 1 1 1 1 1 1 ...
 $ Radiotherapy 0=No, 1= yes, 2= unkown            : num [1:448] 0 0 0 0 0 0 0 0 0 0 ...
 $ Combine0yes1no                                  : num [1:448] 1 1 1 1 1 1 1 1 1 1 ...
 $ Combine rad No, 0, Yes= 1                       : num [1:448] 0 0 0 0 0 0 0 0 0 0 ...
 $ Sex male : 1, female 0                          : num [1:448] 1 0 1 1 1 1 0 1 0 1 ...
 $ Ethnic group 1: malay, 2, chinese, 3 : Indian   : num [1:448] 3 1 1 2 1 1 1 1 1 1 ...
 $ ethnic malay =0. non malay = 1                  : num [1:448] 1 0 0 1 0 0 0 0 0 0 ...
 $ Duration in days                                : num [1:448] 5478 355 1777 286 1801 ...
 $ Duration in years                               : num [1:448] 15.008 0.973 4.868 0.784 4.934 ...
 $ Incidence date                                  : POSIXct[1:448], format: "2008-01-01" "2008-01-01" ...
 $ duration in month Guna (Outcome)                : num [1:448] 180.1 11.7 58.4 9.4 59.2 ...
 $ Updated latest                                  : chr [1:448] "ACTIVE" "DIED" "DIED" "DIED" ...
 $ Alive : 0, Death : 1                            : num [1:448] 0 1 1 1 1 1 1 0 0 1 ...
 $ NoPresent Status, Alive : 1, Died : 2           : num [1:448] 1 2 2 2 2 2 2 1 1 2 ...
 $ Date still alive                                : POSIXct[1:448], format: "2022-12-31" NA ...
 $ Date of died                                    : POSIXct[1:448], format: NA "2008-12-21" ...
Show the Code
library(janitor)
Warning: package 'janitor' was built under R version 4.4.2

Attaching package: 'janitor'
The following objects are masked from 'package:stats':

    chisq.test, fisher.test
Show the Code
surv_data <- datasurv %>% clean_names()
str(surv_data)
tibble [448 × 41] (S3: tbl_df/tbl/data.frame)
 $ age                                          : num [1:448] 38 19 21 60 81 41 69 49 55 48 ...
 $ x60_0_60_1                                   : num [1:448] 0 0 0 1 1 0 1 0 0 0 ...
 $ x15_30_31_59_60                              : num [1:448] 1 0 0 2 2 1 2 1 1 1 ...
 $ age_group_1544_1_4554_2_5564_3_6574_4_75_5   : num [1:448] 1 1 1 3 5 1 4 2 3 2 ...
 $ agegroupbaru_15400_41601_61_752_753          : num [1:448] 0 0 0 1 3 1 2 1 1 1 ...
 $ topography                                   : num [1:448] 779 779 779 779 779 779 779 779 779 779 ...
 $ indolent0aggressive1veryaggressiive2         : num [1:448] 1 0 1 1 0 1 1 1 1 1 ...
 $ very_aggressive0_agrresive1_indolent2        : num [1:448] 1 2 1 1 2 1 1 1 1 1 ...
 $ indolent0_aggressive1very_aggressive2latest  : num [1:448] 1 0 1 1 0 1 1 1 1 1 ...
 $ indolent_0_agressive_1_nos_2                 : num [1:448] 2 0 2 2 0 2 2 2 1 2 ...
 $ x0_other_1_b_cell                            : num [1:448] 1 1 1 1 1 1 1 1 0 1 ...
 $ bcell0_othertcellnkcell1unclassifiable_2     : num [1:448] 1 0 0 0 0 0 0 0 2 0 ...
 $ bcell1othercell0unclassifiable2              : num [1:448] 0 1 1 1 1 1 1 1 2 1 ...
 $ b_cell0_tcell1_precursor_2_unclasifiable3    : num [1:448] 1 0 0 0 0 0 0 0 3 0 ...
 $ x1_b_cell_0_tcell_2_precursor_unclassifiable : num [1:448] 1 1 1 1 1 1 1 1 3 1 ...
 $ aggressive_1_indolent_2                      : num [1:448] 1 2 1 1 2 1 1 1 1 1 ...
 $ morphology                                   : num [1:448] 9591 9691 9591 9591 9695 ...
 $ stage_10_stage21_stage32_satge43             : num [1:448] 2 0 4 2 4 2 4 4 2 4 ...
 $ stage10_stage231_stage42unnkown_3            : num [1:448] 1 0 3 1 3 1 3 3 1 3 ...
 $ stage120_stage31satge42unkown3               : num [1:448] 1 0 3 1 3 1 3 3 1 3 ...
 $ stage_0_stage_1_2_1_stage_3_4_2_unknown_group: num [1:448] 1 0 2 1 2 1 2 2 1 2 ...
 $ stage                                        : num [1:448] 3 1 9 3 9 3 9 9 3 9 ...
 $ icd_10                                       : chr [1:448] "C859" "C821" "C859" "C859" ...
 $ chemotherapyyes0mo1unkown2                   : num [1:448] 0 1 1 0 1 0 1 0 0 1 ...
 $ chemotherapy_0_no_yes_1_2_unknown            : num [1:448] 1 0 0 1 0 1 0 1 1 0 ...
 $ radiotherapy0yes_no1_unknown2                : num [1:448] 1 1 1 1 1 1 1 1 1 1 ...
 $ radiotherapy_0_no_1_yes_2_unkown             : num [1:448] 0 0 0 0 0 0 0 0 0 0 ...
 $ combine0yes1no                               : num [1:448] 1 1 1 1 1 1 1 1 1 1 ...
 $ combine_rad_no_0_yes_1                       : num [1:448] 0 0 0 0 0 0 0 0 0 0 ...
 $ sex_male_1_female_0                          : num [1:448] 1 0 1 1 1 1 0 1 0 1 ...
 $ ethnic_group_1_malay_2_chinese_3_indian      : num [1:448] 3 1 1 2 1 1 1 1 1 1 ...
 $ ethnic_malay_0_non_malay_1                   : num [1:448] 1 0 0 1 0 0 0 0 0 0 ...
 $ duration_in_days                             : num [1:448] 5478 355 1777 286 1801 ...
 $ duration_in_years                            : num [1:448] 15.008 0.973 4.868 0.784 4.934 ...
 $ incidence_date                               : POSIXct[1:448], format: "2008-01-01" "2008-01-01" ...
 $ duration_in_month_guna_outcome               : num [1:448] 180.1 11.7 58.4 9.4 59.2 ...
 $ updated_latest                               : chr [1:448] "ACTIVE" "DIED" "DIED" "DIED" ...
 $ alive_0_death_1                              : num [1:448] 0 1 1 1 1 1 1 0 0 1 ...
 $ no_present_status_alive_1_died_2             : num [1:448] 1 2 2 2 2 2 2 1 1 2 ...
 $ date_still_alive                             : POSIXct[1:448], format: "2022-12-31" NA ...
 $ date_of_died                                 : POSIXct[1:448], format: NA "2008-12-21" ...
Show the Code
surv_data <- surv_data %>% mutate(dur = surv_data$date_of_died %--% surv_data$date_still_alive) %>%
  mutate(dur = as.duration(dur))

surv_data$chemotherapy_0_no_yes_1_2_unknown <- as.factor(surv_data$chemotherapy_0_no_yes_1_2_unknown)
surv_data$x60_0_60_1 <- as.factor(surv_data$x60_0_60_1)
surv_data$stage120_stage31satge42unkown3 <- as.factor(surv_data$stage120_stage31satge42unkown3)
surv_data$alive_0_death_1 <- as.factor(surv_data$alive_0_death_1)
str(surv_data)
tibble [448 × 42] (S3: tbl_df/tbl/data.frame)
 $ age                                          : num [1:448] 38 19 21 60 81 41 69 49 55 48 ...
 $ x60_0_60_1                                   : Factor w/ 2 levels "0","1": 1 1 1 2 2 1 2 1 1 1 ...
 $ x15_30_31_59_60                              : num [1:448] 1 0 0 2 2 1 2 1 1 1 ...
 $ age_group_1544_1_4554_2_5564_3_6574_4_75_5   : num [1:448] 1 1 1 3 5 1 4 2 3 2 ...
 $ agegroupbaru_15400_41601_61_752_753          : num [1:448] 0 0 0 1 3 1 2 1 1 1 ...
 $ topography                                   : num [1:448] 779 779 779 779 779 779 779 779 779 779 ...
 $ indolent0aggressive1veryaggressiive2         : num [1:448] 1 0 1 1 0 1 1 1 1 1 ...
 $ very_aggressive0_agrresive1_indolent2        : num [1:448] 1 2 1 1 2 1 1 1 1 1 ...
 $ indolent0_aggressive1very_aggressive2latest  : num [1:448] 1 0 1 1 0 1 1 1 1 1 ...
 $ indolent_0_agressive_1_nos_2                 : num [1:448] 2 0 2 2 0 2 2 2 1 2 ...
 $ x0_other_1_b_cell                            : num [1:448] 1 1 1 1 1 1 1 1 0 1 ...
 $ bcell0_othertcellnkcell1unclassifiable_2     : num [1:448] 1 0 0 0 0 0 0 0 2 0 ...
 $ bcell1othercell0unclassifiable2              : num [1:448] 0 1 1 1 1 1 1 1 2 1 ...
 $ b_cell0_tcell1_precursor_2_unclasifiable3    : num [1:448] 1 0 0 0 0 0 0 0 3 0 ...
 $ x1_b_cell_0_tcell_2_precursor_unclassifiable : num [1:448] 1 1 1 1 1 1 1 1 3 1 ...
 $ aggressive_1_indolent_2                      : num [1:448] 1 2 1 1 2 1 1 1 1 1 ...
 $ morphology                                   : num [1:448] 9591 9691 9591 9591 9695 ...
 $ stage_10_stage21_stage32_satge43             : num [1:448] 2 0 4 2 4 2 4 4 2 4 ...
 $ stage10_stage231_stage42unnkown_3            : num [1:448] 1 0 3 1 3 1 3 3 1 3 ...
 $ stage120_stage31satge42unkown3               : Factor w/ 4 levels "0","1","2","3": 2 1 4 2 4 2 4 4 2 4 ...
 $ stage_0_stage_1_2_1_stage_3_4_2_unknown_group: num [1:448] 1 0 2 1 2 1 2 2 1 2 ...
 $ stage                                        : num [1:448] 3 1 9 3 9 3 9 9 3 9 ...
 $ icd_10                                       : chr [1:448] "C859" "C821" "C859" "C859" ...
 $ chemotherapyyes0mo1unkown2                   : num [1:448] 0 1 1 0 1 0 1 0 0 1 ...
 $ chemotherapy_0_no_yes_1_2_unknown            : Factor w/ 3 levels "0","1","2": 2 1 1 2 1 2 1 2 2 1 ...
 $ radiotherapy0yes_no1_unknown2                : num [1:448] 1 1 1 1 1 1 1 1 1 1 ...
 $ radiotherapy_0_no_1_yes_2_unkown             : num [1:448] 0 0 0 0 0 0 0 0 0 0 ...
 $ combine0yes1no                               : num [1:448] 1 1 1 1 1 1 1 1 1 1 ...
 $ combine_rad_no_0_yes_1                       : num [1:448] 0 0 0 0 0 0 0 0 0 0 ...
 $ sex_male_1_female_0                          : num [1:448] 1 0 1 1 1 1 0 1 0 1 ...
 $ ethnic_group_1_malay_2_chinese_3_indian      : num [1:448] 3 1 1 2 1 1 1 1 1 1 ...
 $ ethnic_malay_0_non_malay_1                   : num [1:448] 1 0 0 1 0 0 0 0 0 0 ...
 $ duration_in_days                             : num [1:448] 5478 355 1777 286 1801 ...
 $ duration_in_years                            : num [1:448] 15.008 0.973 4.868 0.784 4.934 ...
 $ incidence_date                               : POSIXct[1:448], format: "2008-01-01" "2008-01-01" ...
 $ duration_in_month_guna_outcome               : num [1:448] 180.1 11.7 58.4 9.4 59.2 ...
 $ updated_latest                               : chr [1:448] "ACTIVE" "DIED" "DIED" "DIED" ...
 $ alive_0_death_1                              : Factor w/ 2 levels "0","1": 1 2 2 2 2 2 2 1 1 2 ...
 $ no_present_status_alive_1_died_2             : num [1:448] 1 2 2 2 2 2 2 1 1 2 ...
 $ date_still_alive                             : POSIXct[1:448], format: "2022-12-31" NA ...
 $ date_of_died                                 : POSIXct[1:448], format: NA "2008-12-21" ...
 $ dur                                          :Formal class 'Duration' [package "lubridate"] with 1 slot
  .. ..@ .Data: num [1:448] NA NA NA NA NA NA NA NA NA NA ...
Show the Code
KM1 <- survfit(Surv(time = surv_data$duration_in_month_guna_outcome, event = alive_0_death_1 == '1') ~ 1,
               type = "kaplan-meier", data = surv_data)
summary(KM1)
Call: survfit(formula = Surv(time = surv_data$duration_in_month_guna_outcome, 
    event = alive_0_death_1 == "1") ~ 1, data = surv_data, type = "kaplan-meier")

     time n.risk n.event survival std.err lower 95% CI upper 95% CI
   0.0000    448      10    0.978 0.00698        0.964        0.991
   0.0329    438       2    0.973 0.00763        0.958        0.988
   0.0658    436       3    0.967 0.00850        0.950        0.983
   0.0986    433       7    0.951 0.01021        0.931        0.971
   0.1315    426       3    0.944 0.01084        0.923        0.966
   0.1644    423       4    0.935 0.01162        0.913        0.958
   0.1973    419       4    0.926 0.01234        0.902        0.951
   0.2301    415       3    0.920 0.01284        0.895        0.945
   0.2630    412       5    0.908 0.01362        0.882        0.936
   0.2959    407       3    0.902 0.01406        0.875        0.930
   0.3288    404       3    0.895 0.01448        0.867        0.924
   0.3616    401       4    0.886 0.01501        0.857        0.916
   0.4274    397       4    0.877 0.01550        0.847        0.908
   0.4603    393       3    0.871 0.01586        0.840        0.902
   0.4932    390       2    0.866 0.01609        0.835        0.898
   0.5260    388       4    0.857 0.01653        0.825        0.890
   0.5589    384       7    0.842 0.01725        0.808        0.876
   0.5918    377       3    0.835 0.01754        0.801        0.870
   0.6247    374       1    0.833 0.01764        0.799        0.868
   0.6575    373       2    0.828 0.01782        0.794        0.864
   0.6904    371       3    0.821 0.01809        0.787        0.858
   0.7562    368       2    0.817 0.01827        0.782        0.854
   0.7890    366       6    0.804 0.01877        0.768        0.841
   0.8219    360       2    0.799 0.01893        0.763        0.837
   0.8548    358       2    0.795 0.01909        0.758        0.833
   0.8877    356       1    0.792 0.01916        0.756        0.831
   0.9205    355       3    0.786 0.01939        0.749        0.825
   0.9534    352       1    0.783 0.01946        0.746        0.823
   0.9863    351       2    0.779 0.01960        0.742        0.818
   1.0192    349       2    0.775 0.01974        0.737        0.814
   1.0521    347       2    0.770 0.01988        0.732        0.810
   1.0849    345       2    0.766 0.02001        0.727        0.806
   1.1836    343       1    0.763 0.02008        0.725        0.804
   1.2493    342       2    0.759 0.02021        0.720        0.800
   1.2822    340       1    0.757 0.02027        0.718        0.797
   1.3151    339       1    0.754 0.02033        0.716        0.795
   1.4137    338       3    0.748 0.02052        0.709        0.789
   1.4466    335       2    0.743 0.02064        0.704        0.785
   1.4795    333       1    0.741 0.02070        0.702        0.783
   1.5123    332       1    0.739 0.02075        0.699        0.781
   1.6767    331       1    0.737 0.02081        0.697        0.779
   1.7096    330       1    0.734 0.02087        0.695        0.776
   1.7753    329       2    0.730 0.02098        0.690        0.772
   1.8082    327       1    0.728 0.02103        0.688        0.770
   1.8740    326       1    0.725 0.02109        0.685        0.768
   1.9068    325       1    0.723 0.02114        0.683        0.766
   1.9397    324       3    0.717 0.02129        0.676        0.759
   1.9726    321       1    0.714 0.02134        0.674        0.757
   2.0384    320       1    0.712 0.02139        0.671        0.755
   2.1041    319       3    0.705 0.02154        0.664        0.749
   2.1370    316       2    0.701 0.02163        0.660        0.745
   2.2027    314       2    0.696 0.02172        0.655        0.740
   2.2356    312       1    0.694 0.02177        0.653        0.738
   2.2685    311       1    0.692 0.02181        0.651        0.736
   2.3014    310       1    0.690 0.02186        0.648        0.734
   2.3342    309       1    0.687 0.02190        0.646        0.732
   2.4000    308       1    0.685 0.02194        0.644        0.730
   2.4658    307       1    0.683 0.02198        0.641        0.728
   2.5315    306       1    0.681 0.02202        0.639        0.725
   2.5644    305       1    0.679 0.02206        0.637        0.723
   2.6301    304       2    0.674 0.02214        0.632        0.719
   2.7288    302       1    0.672 0.02218        0.630        0.717
   2.7945    301       1    0.670 0.02222        0.627        0.715
   2.8932    300       2    0.665 0.02230        0.623        0.710
   2.9918    298       3    0.658 0.02240        0.616        0.704
   3.0247    295       1    0.656 0.02244        0.614        0.702
   3.1233    294       1    0.654 0.02247        0.611        0.700
   3.2219    293       1    0.652 0.02251        0.609        0.697
   3.2877    292       2    0.647 0.02257        0.605        0.693
   3.3534    290       1    0.645 0.02261        0.602        0.691
   3.4521    289       1    0.643 0.02264        0.600        0.689
   3.4849    288       2    0.638 0.02270        0.595        0.684
   3.5178    286       2    0.634 0.02276        0.591        0.680
   3.6164    284       1    0.632 0.02279        0.589        0.678
   3.6493    283       1    0.629 0.02282        0.586        0.676
   3.8137    282       3    0.623 0.02290        0.579        0.669
   3.8466    279       1    0.621 0.02293        0.577        0.667
   3.8795    278       2    0.616 0.02298        0.573        0.663
   3.9123    276       2    0.612 0.02303        0.568        0.658
   3.9781    274       2    0.607 0.02307        0.564        0.654
   4.0767    272       1    0.605 0.02310        0.561        0.652
   4.1096    271       1    0.603 0.02312        0.559        0.650
   4.1425    270       1    0.600 0.02314        0.557        0.648
   4.1753    269       2    0.596 0.02318        0.552        0.643
   4.2411    267       3    0.589 0.02324        0.545        0.637
   4.4712    264       1    0.587 0.02326        0.543        0.634
   4.5699    263       1    0.585 0.02328        0.541        0.632
   4.6027    262       2    0.580 0.02332        0.536        0.628
   4.6356    260       1    0.578 0.02333        0.534        0.626
   4.7342    259       1    0.576 0.02335        0.532        0.624
   4.9315    258       1    0.574 0.02337        0.530        0.621
   5.0959    257       1    0.571 0.02338        0.527        0.619
   5.1616    256       1    0.569 0.02340        0.525        0.617
   5.2603    255       1    0.567 0.02341        0.523        0.615
   5.4575    254       1    0.565 0.02342        0.521        0.613
   5.5562    253       1    0.562 0.02344        0.518        0.610
   5.6219    252       1    0.560 0.02345        0.516        0.608
   6.0493    251       2    0.556 0.02348        0.512        0.604
   6.0822    249       1    0.554 0.02349        0.509        0.602
   6.1479    248       1    0.551 0.02350        0.507        0.599
   6.2466    247       1    0.549 0.02351        0.505        0.597
   6.3123    246       1    0.547 0.02352        0.503        0.595
   6.8384    245       1    0.545 0.02353        0.500        0.593
   6.9041    244       1    0.542 0.02354        0.498        0.591
   7.2986    243       1    0.540 0.02355        0.496        0.588
   7.5945    242       1    0.538 0.02355        0.494        0.586
   7.6274    241       1    0.536 0.02356        0.491        0.584
   7.6932    240       1    0.533 0.02357        0.489        0.582
   7.7260    239       1    0.531 0.02358        0.487        0.580
   7.8575    238       1    0.529 0.02358        0.485        0.577
   8.0548    237       1    0.527 0.02359        0.483        0.575
   8.0877    236       1    0.525 0.02359        0.480        0.573
   8.1205    235       2    0.520 0.02360        0.476        0.568
   8.3178    233       1    0.518 0.02361        0.474        0.566
   8.5479    232       1    0.516 0.02361        0.471        0.564
   8.6137    231       1    0.513 0.02361        0.469        0.562
   8.6466    230       1    0.511 0.02362        0.467        0.560
   8.6795    229       1    0.509 0.02362        0.465        0.557
   8.8438    228       2    0.504 0.02362        0.460        0.553
   8.9425    226       1    0.502 0.02362        0.458        0.551
   9.0082    225       1    0.500 0.02362        0.456        0.549
   9.1397    224       1    0.498 0.02362        0.454        0.546
   9.4027    223       1    0.496 0.02362        0.451        0.544
   9.5342    222       1    0.493 0.02362        0.449        0.542
   9.9616    221       1    0.491 0.02362        0.447        0.540
  10.0932    220       1    0.489 0.02362        0.445        0.537
  10.9151    219       1    0.487 0.02361        0.442        0.535
  11.0137    218       1    0.484 0.02361        0.440        0.533
  11.0795    217       1    0.482 0.02361        0.438        0.531
  11.1452    216       1    0.480 0.02360        0.436        0.528
  11.1781    215       2    0.475 0.02359        0.431        0.524
  11.2110    213       1    0.473 0.02359        0.429        0.522
  11.6712    212       3    0.467 0.02357        0.423        0.515
  12.0000    209       2    0.462 0.02355        0.418        0.511
  12.0329    207       1    0.460 0.02355        0.416        0.508
  12.0986    206       1    0.458 0.02354        0.414        0.506
  12.2301    205       1    0.455 0.02353        0.412        0.504
  12.4274    204       1    0.453 0.02352        0.409        0.502
  12.8219    203       1    0.451 0.02351        0.407        0.499
  12.9205    202       1    0.449 0.02350        0.405        0.497
  13.4466    201       1    0.446 0.02349        0.403        0.495
  13.4795    200       1    0.444 0.02348        0.400        0.493
  14.5644    199       2    0.440 0.02345        0.396        0.488
  15.1562    197       1    0.438 0.02344        0.394        0.486
  15.3205    196       2    0.433 0.02341        0.390        0.481
  15.5507    194       1    0.431 0.02340        0.387        0.479
  15.8137    193       2    0.426 0.02337        0.383        0.475
  16.4055    191       1    0.424 0.02335        0.381        0.472
  16.7014    190       1    0.422 0.02333        0.379        0.470
  16.7342    189       1    0.420 0.02332        0.376        0.468
  16.8000    188       1    0.417 0.02330        0.374        0.466
  17.0959    187       1    0.415 0.02328        0.372        0.463
  17.2932    186       1    0.413 0.02326        0.370        0.461
  17.4575    185       1    0.411 0.02324        0.368        0.459
  17.6219    184       1    0.408 0.02322        0.365        0.457
  18.0493    183       1    0.406 0.02320        0.363        0.454
  19.8904    182       1    0.404 0.02318        0.361        0.452
  20.4164    181       1    0.402 0.02316        0.359        0.450
  20.5151    180       1    0.400 0.02314        0.357        0.448
  21.0082    179       1    0.397 0.02312        0.354        0.445
  21.8630    178       1    0.395 0.02310        0.352        0.443
  21.9945    177       1    0.393 0.02307        0.350        0.441
  22.3233    176       1    0.391 0.02305        0.348        0.439
  22.5534    175       2    0.386 0.02300        0.344        0.434
  22.9479    173       1    0.384 0.02298        0.341        0.432
  24.5260    172       1    0.382 0.02295        0.339        0.429
  24.5918    171       1    0.379 0.02293        0.337        0.427
  24.8219    170       1    0.377 0.02290        0.335        0.425
  24.8548    169       1    0.375 0.02287        0.333        0.423
  24.9863    168       1    0.373 0.02285        0.331        0.420
  25.2164    167       1    0.371 0.02282        0.328        0.418
  25.6767    166       1    0.368 0.02279        0.326        0.416
  26.1041    165       1    0.366 0.02276        0.324        0.414
  26.1699    164       1    0.364 0.02273        0.322        0.411
  26.4329    163       1    0.362 0.02270        0.320        0.409
  27.1562    162       1    0.359 0.02267        0.318        0.407
  27.5507    161       1    0.357 0.02264        0.315        0.404
  28.1425    160       1    0.355 0.02261        0.313        0.402
  29.9507    159       1    0.353 0.02257        0.311        0.400
  30.8055    158       1    0.350 0.02254        0.309        0.398
  31.3973    157       1    0.348 0.02251        0.307        0.395
  31.5616    156       1    0.346 0.02247        0.305        0.393
  31.7918    155       1    0.344 0.02244        0.302        0.391
  33.1397    154       1    0.342 0.02240        0.300        0.388
  33.9945    153       1    0.339 0.02237        0.298        0.386
  34.7178    152       1    0.337 0.02233        0.296        0.384
  35.2767    151       1    0.335 0.02230        0.294        0.382
  36.5918    150       1    0.333 0.02226        0.292        0.379
  36.6904    149       1    0.330 0.02222        0.290        0.377
  38.0712    148       1    0.328 0.02218        0.287        0.375
  40.6356    147       1    0.326 0.02214        0.285        0.372
  41.5233    146       1    0.324 0.02210        0.283        0.370
  42.8384    145       1    0.321 0.02206        0.281        0.368
  43.3644    144       1    0.319 0.02202        0.279        0.365
  43.4301    143       1    0.317 0.02198        0.277        0.363
  48.0658    142       1    0.315 0.02194        0.275        0.361
  50.7945    141       1    0.313 0.02190        0.272        0.359
  51.7479    140       1    0.310 0.02186        0.270        0.356
  55.2000    139       1    0.308 0.02181        0.268        0.354
  55.2658    138       1    0.306 0.02177        0.266        0.352
  55.2986    137       1    0.304 0.02172        0.264        0.349
  57.0740    136       1    0.301 0.02168        0.262        0.347
  57.8301    135       2    0.297 0.02159        0.257        0.342
  58.4219    133       2    0.292 0.02149        0.253        0.338
  58.5534    131       1    0.290 0.02144        0.251        0.335
  59.2110    130       1    0.288 0.02139        0.249        0.333
  59.8027    129       1    0.286 0.02134        0.247        0.331
  59.8356    128       1    0.283 0.02129        0.245        0.328
  60.8877    126       1    0.281 0.02124        0.243        0.326
  61.6767    125       1    0.279 0.02119        0.240        0.324
  62.4329    123       1    0.277 0.02114        0.238        0.321
  66.4110    113       1    0.274 0.02109        0.236        0.319
  66.5753    112       1    0.272 0.02105        0.234        0.316
  66.6082    111       1    0.269 0.02100        0.231        0.314
  68.1205    108       1    0.267 0.02095        0.229        0.311
  73.2493    102       1    0.264 0.02091        0.226        0.309
  77.7534     97       1    0.262 0.02087        0.224        0.306
  79.0027     95       1    0.259 0.02083        0.221        0.303
  80.5479     91       1    0.256 0.02080        0.218        0.300
  81.0740     89       1    0.253 0.02076        0.215        0.297
  82.1260     84       1    0.250 0.02073        0.213        0.294
  83.4082     80       1    0.247 0.02071        0.209        0.291
  94.8164     62       1    0.243 0.02075        0.205        0.287
  97.6438     57       1    0.239 0.02082        0.201        0.283
  98.2027     56       1    0.234 0.02088        0.197        0.279
 103.5945     50       1    0.230 0.02098        0.192        0.275
 104.0548     48       1    0.225 0.02108        0.187        0.270
 112.6027     46       1    0.220 0.02119        0.182        0.266
 115.9562     43       1    0.215 0.02130        0.177        0.261
 118.9808     41       1    0.210 0.02142        0.172        0.256
 121.8411     39       1    0.204 0.02153        0.166        0.251
 125.6877     35       1    0.198 0.02169        0.160        0.246
 167.7041     12       1    0.182 0.02542        0.138        0.239
 175.0356      7       1    0.156 0.03246        0.104        0.235
Show the Code
ggsurvplot(KM1, data = surv_data, risk.table = TRUE)

Show the Code
KM1.age <- survfit(Surv(time = surv_data$duration_in_month_guna_outcome, event = alive_0_death_1 == '1') ~ x60_0_60_1,
               type = "kaplan-meier", data = surv_data)
summary(KM1.age)
Call: survfit(formula = Surv(time = surv_data$duration_in_month_guna_outcome, 
    event = alive_0_death_1 == "1") ~ x60_0_60_1, data = surv_data, 
    type = "kaplan-meier")

                x60_0_60_1=0 
     time n.risk n.event survival std.err lower 95% CI upper 95% CI
   0.0000    288       6    0.979 0.00842        0.963        0.996
   0.0658    282       2    0.972 0.00968        0.953        0.991
   0.0986    280       6    0.951 0.01267        0.927        0.977
   0.1315    274       2    0.944 0.01350        0.918        0.971
   0.1644    272       2    0.938 0.01426        0.910        0.966
   0.1973    270       3    0.927 0.01532        0.898        0.958
   0.2301    267       2    0.920 0.01597        0.889        0.952
   0.2630    265       2    0.913 0.01659        0.881        0.946
   0.2959    263       1    0.910 0.01689        0.877        0.943
   0.3288    262       1    0.906 0.01718        0.873        0.941
   0.3616    261       2    0.899 0.01773        0.865        0.935
   0.4274    259       3    0.889 0.01852        0.853        0.926
   0.4603    256       1    0.885 0.01877        0.849        0.923
   0.4932    255       2    0.878 0.01925        0.842        0.917
   0.5260    253       3    0.868 0.01994        0.830        0.908
   0.5589    250       6    0.847 0.02120        0.807        0.890
   0.5918    244       1    0.844 0.02140        0.803        0.887
   0.6247    243       1    0.840 0.02159        0.799        0.884
   0.6575    242       2    0.833 0.02196        0.791        0.878
   0.6904    240       2    0.826 0.02232        0.784        0.871
   0.7562    238       1    0.823 0.02249        0.780        0.868
   0.7890    237       2    0.816 0.02283        0.772        0.862
   0.8219    235       2    0.809 0.02316        0.765        0.856
   0.9205    233       1    0.806 0.02332        0.761        0.853
   0.9863    232       2    0.799 0.02363        0.754        0.846
   1.0192    230       2    0.792 0.02393        0.746        0.840
   1.0521    228       1    0.788 0.02408        0.742        0.837
   1.0849    227       2    0.781 0.02436        0.735        0.830
   1.2493    225       1    0.778 0.02450        0.731        0.827
   1.4137    224       2    0.771 0.02477        0.724        0.821
   1.5123    222       1    0.767 0.02490        0.720        0.818
   1.6767    221       1    0.764 0.02503        0.716        0.815
   1.7096    220       1    0.760 0.02515        0.713        0.811
   1.7753    219       1    0.757 0.02527        0.709        0.808
   1.9397    218       3    0.747 0.02563        0.698        0.798
   2.0384    215       1    0.743 0.02575        0.694        0.795
   2.1041    214       1    0.740 0.02586        0.691        0.792
   2.1370    213       1    0.736 0.02597        0.687        0.789
   2.2027    212       2    0.729 0.02619        0.680        0.782
   2.2356    210       1    0.726 0.02629        0.676        0.779
   2.2685    209       1    0.722 0.02639        0.672        0.776
   2.3342    208       1    0.719 0.02649        0.669        0.773
   2.4658    207       1    0.715 0.02659        0.665        0.769
   2.5644    206       1    0.712 0.02669        0.661        0.766
   2.6301    205       1    0.708 0.02678        0.658        0.763
   2.7288    204       1    0.705 0.02688        0.654        0.760
   2.8932    203       1    0.701 0.02697        0.650        0.756
   3.1233    202       1    0.698 0.02706        0.647        0.753
   3.2877    201       1    0.694 0.02714        0.643        0.750
   3.4849    200       2    0.688 0.02731        0.636        0.743
   3.5178    198       1    0.684 0.02739        0.632        0.740
   3.6164    197       1    0.681 0.02747        0.629        0.737
   3.6493    196       1    0.677 0.02755        0.625        0.733
   3.8466    195       1    0.674 0.02763        0.622        0.730
   3.8795    194       1    0.670 0.02770        0.618        0.727
   3.9123    193       1    0.667 0.02778        0.614        0.723
   3.9781    192       1    0.663 0.02785        0.611        0.720
   4.0767    191       1    0.660 0.02792        0.607        0.717
   4.1096    190       1    0.656 0.02799        0.604        0.713
   4.2411    189       1    0.653 0.02805        0.600        0.710
   4.4712    188       1    0.649 0.02812        0.596        0.707
   4.6027    187       1    0.646 0.02818        0.593        0.703
   4.6356    186       1    0.642 0.02824        0.589        0.700
   4.7342    185       1    0.639 0.02830        0.586        0.697
   5.0959    184       1    0.635 0.02836        0.582        0.694
   5.1616    183       1    0.632 0.02842        0.579        0.690
   5.2603    182       1    0.628 0.02847        0.575        0.687
   5.4575    181       1    0.625 0.02853        0.572        0.683
   5.5562    180       1    0.622 0.02858        0.568        0.680
   5.6219    179       1    0.618 0.02863        0.564        0.677
   6.0493    178       1    0.615 0.02868        0.561        0.673
   6.1479    177       1    0.611 0.02873        0.557        0.670
   6.2466    176       1    0.608 0.02877        0.554        0.667
   6.8384    175       1    0.604 0.02882        0.550        0.663
   6.9041    174       1    0.601 0.02886        0.547        0.660
   7.2986    173       1    0.597 0.02890        0.543        0.657
   7.5945    172       1    0.594 0.02894        0.540        0.653
   7.6274    171       1    0.590 0.02898        0.536        0.650
   7.6932    170       1    0.587 0.02902        0.533        0.647
   7.7260    169       1    0.583 0.02905        0.529        0.643
   7.8575    168       1    0.580 0.02908        0.526        0.640
   8.0548    167       1    0.576 0.02912        0.522        0.636
   8.0877    166       1    0.573 0.02915        0.519        0.633
   8.1205    165       1    0.569 0.02918        0.515        0.630
   8.3178    164       1    0.566 0.02921        0.512        0.626
   8.6137    163       1    0.563 0.02923        0.508        0.623
   8.6466    162       1    0.559 0.02926        0.505        0.619
   8.8438    161       1    0.556 0.02928        0.501        0.616
   9.0082    160       1    0.552 0.02930        0.498        0.613
   9.1397    159       1    0.549 0.02932        0.494        0.609
   9.5342    158       1    0.545 0.02934        0.491        0.606
   9.9616    157       1    0.542 0.02936        0.487        0.602
  10.0932    156       1    0.538 0.02938        0.484        0.599
  11.0795    155       1    0.535 0.02939        0.480        0.596
  11.1452    154       1    0.531 0.02941        0.477        0.592
  11.1781    153       1    0.528 0.02942        0.473        0.589
  11.6712    152       2    0.521 0.02944        0.466        0.582
  12.0000    150       2    0.514 0.02945        0.459        0.575
  12.0986    148       1    0.510 0.02946        0.456        0.572
  12.2301    147       1    0.507 0.02946        0.452        0.568
  12.4274    146       1    0.503 0.02946        0.449        0.565
  13.4466    145       1    0.500 0.02946        0.445        0.561
  14.5644    144       1    0.497 0.02946        0.442        0.558
  15.5507    143       1    0.493 0.02946        0.439        0.554
  16.4055    142       1    0.490 0.02946        0.435        0.551
  16.8000    141       1    0.486 0.02945        0.432        0.547
  17.0959    140       1    0.483 0.02945        0.428        0.544
  17.2932    139       1    0.479 0.02944        0.425        0.540
  17.4575    138       1    0.476 0.02943        0.421        0.537
  17.6219    137       1    0.472 0.02942        0.418        0.534
  18.0493    136       1    0.469 0.02941        0.415        0.530
  20.5151    135       1    0.465 0.02939        0.411        0.527
  21.0082    134       1    0.462 0.02938        0.408        0.523
  21.8630    133       1    0.458 0.02936        0.404        0.520
  22.3233    132       1    0.455 0.02934        0.401        0.516
  22.5534    131       2    0.448 0.02930        0.394        0.509
  22.9479    129       1    0.444 0.02928        0.391        0.506
  24.5260    128       1    0.441 0.02926        0.387        0.502
  24.5918    127       1    0.438 0.02923        0.384        0.499
  24.8548    126       1    0.434 0.02921        0.380        0.495
  24.9863    125       1    0.431 0.02918        0.377        0.492
  25.2164    124       1    0.427 0.02915        0.374        0.488
  25.6767    123       1    0.424 0.02912        0.370        0.485
  26.1041    122       1    0.420 0.02908        0.367        0.481
  26.4329    121       1    0.417 0.02905        0.363        0.478
  27.1562    120       1    0.413 0.02902        0.360        0.474
  27.5507    119       1    0.410 0.02898        0.357        0.471
  30.8055    118       1    0.406 0.02894        0.353        0.467
  31.3973    117       1    0.403 0.02890        0.350        0.464
  31.7918    116       1    0.399 0.02886        0.347        0.460
  33.9945    115       1    0.396 0.02882        0.343        0.457
  34.7178    114       1    0.392 0.02877        0.340        0.453
  35.2767    113       1    0.389 0.02873        0.336        0.449
  36.6904    112       1    0.385 0.02868        0.333        0.446
  38.0712    111       1    0.382 0.02863        0.330        0.442
  40.6356    110       1    0.378 0.02858        0.326        0.439
  41.5233    109       1    0.375 0.02853        0.323        0.435
  42.8384    108       1    0.372 0.02847        0.320        0.432
  43.4301    107       1    0.368 0.02842        0.316        0.428
  48.0658    106       1    0.365 0.02836        0.313        0.425
  50.7945    105       1    0.361 0.02830        0.310        0.421
  55.2658    104       1    0.358 0.02824        0.306        0.418
  55.2986    103       1    0.354 0.02818        0.303        0.414
  57.0740    102       1    0.351 0.02812        0.300        0.410
  58.4219    101       1    0.347 0.02805        0.296        0.407
  59.8027    100       1    0.344 0.02799        0.293        0.403
  59.8356     99       1    0.340 0.02792        0.290        0.400
  66.4110     91       1    0.337 0.02786        0.286        0.396
  66.5753     90       1    0.333 0.02780        0.283        0.392
  73.2493     83       1    0.329 0.02775        0.279        0.388
  77.7534     80       1    0.325 0.02771        0.275        0.384
  80.5479     76       1    0.320 0.02767        0.271        0.380
  81.0740     74       1    0.316 0.02764        0.266        0.375
  82.1260     69       1    0.311 0.02761        0.262        0.371
  83.4082     66       1    0.307 0.02759        0.257        0.366
  94.8164     50       1    0.301 0.02772        0.251        0.360
  98.2027     47       1    0.294 0.02785        0.244        0.354
 103.5945     42       1    0.287 0.02806        0.237        0.348
 104.0548     40       1    0.280 0.02826        0.230        0.341
 115.9562     36       1    0.272 0.02853        0.222        0.334
 121.8411     33       1    0.264 0.02883        0.213        0.327
 167.7041     10       1    0.238 0.03606        0.176        0.320
 175.0356      6       1    0.198 0.04701        0.124        0.315

                x60_0_60_1=1 
     time n.risk n.event survival std.err lower 95% CI upper 95% CI
   0.0000    160       4   0.9750  0.0123       0.9511        0.999
   0.0329    156       2   0.9625  0.0150       0.9335        0.992
   0.0658    154       1   0.9563  0.0162       0.9251        0.988
   0.0986    153       1   0.9500  0.0172       0.9168        0.984
   0.1315    152       1   0.9438  0.0182       0.9087        0.980
   0.1644    151       2   0.9313  0.0200       0.8929        0.971
   0.1973    149       1   0.9250  0.0208       0.8851        0.967
   0.2301    148       1   0.9188  0.0216       0.8774        0.962
   0.2630    147       3   0.9000  0.0237       0.8547        0.948
   0.2959    144       2   0.8875  0.0250       0.8399        0.938
   0.3288    142       2   0.8750  0.0261       0.8252        0.928
   0.3616    140       2   0.8625  0.0272       0.8108        0.918
   0.4274    138       1   0.8563  0.0277       0.8036        0.912
   0.4603    137       2   0.8438  0.0287       0.7893        0.902
   0.5260    135       1   0.8375  0.0292       0.7822        0.897
   0.5589    134       1   0.8313  0.0296       0.7752        0.891
   0.5918    133       2   0.8188  0.0305       0.7612        0.881
   0.6904    131       1   0.8125  0.0309       0.7542        0.875
   0.7562    130       1   0.8063  0.0312       0.7473        0.870
   0.7890    129       4   0.7813  0.0327       0.7198        0.848
   0.8548    125       2   0.7688  0.0333       0.7061        0.837
   0.8877    123       1   0.7625  0.0336       0.6993        0.831
   0.9205    122       2   0.7500  0.0342       0.6858        0.820
   0.9534    120       1   0.7438  0.0345       0.6791        0.815
   1.0521    119       1   0.7375  0.0348       0.6724        0.809
   1.1836    118       1   0.7313  0.0350       0.6657        0.803
   1.2493    117       1   0.7250  0.0353       0.6590        0.798
   1.2822    116       1   0.7188  0.0355       0.6524        0.792
   1.3151    115       1   0.7125  0.0358       0.6457        0.786
   1.4137    114       1   0.7063  0.0360       0.6391        0.780
   1.4466    113       2   0.6938  0.0364       0.6259        0.769
   1.4795    111       1   0.6875  0.0366       0.6193        0.763
   1.7753    110       1   0.6813  0.0368       0.6127        0.757
   1.8082    109       1   0.6750  0.0370       0.6062        0.752
   1.8740    108       1   0.6688  0.0372       0.5997        0.746
   1.9068    107       1   0.6625  0.0374       0.5931        0.740
   1.9726    106       1   0.6563  0.0375       0.5866        0.734
   2.1041    105       2   0.6438  0.0379       0.5737        0.722
   2.1370    103       1   0.6375  0.0380       0.5672        0.717
   2.3014    102       1   0.6313  0.0381       0.5607        0.711
   2.4000    101       1   0.6250  0.0383       0.5543        0.705
   2.5315    100       1   0.6188  0.0384       0.5479        0.699
   2.6301     99       1   0.6125  0.0385       0.5415        0.693
   2.7945     98       1   0.6063  0.0386       0.5351        0.687
   2.8932     97       1   0.6000  0.0387       0.5287        0.681
   2.9918     96       3   0.5813  0.0390       0.5096        0.663
   3.0247     93       1   0.5750  0.0391       0.5033        0.657
   3.2219     92       1   0.5688  0.0392       0.4970        0.651
   3.2877     91       1   0.5625  0.0392       0.4907        0.645
   3.3534     90       1   0.5563  0.0393       0.4844        0.639
   3.4521     89       1   0.5500  0.0393       0.4781        0.633
   3.5178     88       1   0.5438  0.0394       0.4718        0.627
   3.8137     87       3   0.5250  0.0395       0.4531        0.608
   3.8795     84       1   0.5188  0.0395       0.4468        0.602
   3.9123     83       1   0.5125  0.0395       0.4406        0.596
   3.9781     82       1   0.5063  0.0395       0.4344        0.590
   4.1425     81       1   0.5000  0.0395       0.4282        0.584
   4.1753     80       2   0.4875  0.0395       0.4159        0.571
   4.2411     78       2   0.4750  0.0395       0.4036        0.559
   4.5699     76       1   0.4688  0.0395       0.3975        0.553
   4.6027     75       1   0.4625  0.0394       0.3914        0.547
   4.9315     74       1   0.4563  0.0394       0.3852        0.540
   6.0493     73       1   0.4500  0.0393       0.3792        0.534
   6.0822     72       1   0.4438  0.0393       0.3731        0.528
   6.3123     71       1   0.4375  0.0392       0.3670        0.522
   8.1205     70       1   0.4313  0.0392       0.3610        0.515
   8.5479     69       1   0.4250  0.0391       0.3549        0.509
   8.6795     68       1   0.4188  0.0390       0.3489        0.503
   8.8438     67       1   0.4125  0.0389       0.3429        0.496
   8.9425     66       1   0.4063  0.0388       0.3369        0.490
   9.4027     65       1   0.4000  0.0387       0.3309        0.484
  10.9151     64       1   0.3938  0.0386       0.3249        0.477
  11.0137     63       1   0.3875  0.0385       0.3189        0.471
  11.1781     62       1   0.3813  0.0384       0.3130        0.464
  11.2110     61       1   0.3750  0.0383       0.3070        0.458
  11.6712     60       1   0.3688  0.0381       0.3011        0.452
  12.0329     59       1   0.3625  0.0380       0.2952        0.445
  12.8219     58       1   0.3563  0.0379       0.2893        0.439
  12.9205     57       1   0.3500  0.0377       0.2834        0.432
  13.4795     56       1   0.3438  0.0375       0.2775        0.426
  14.5644     55       1   0.3375  0.0374       0.2716        0.419
  15.1562     54       1   0.3313  0.0372       0.2658        0.413
  15.3205     53       2   0.3188  0.0368       0.2541        0.400
  15.8137     51       2   0.3063  0.0364       0.2425        0.387
  16.7014     49       1   0.3000  0.0362       0.2368        0.380
  16.7342     48       1   0.2938  0.0360       0.2310        0.374
  19.8904     47       1   0.2875  0.0358       0.2253        0.367
  20.4164     46       1   0.2813  0.0355       0.2195        0.360
  21.9945     45       1   0.2750  0.0353       0.2138        0.354
  24.8219     44       1   0.2688  0.0350       0.2081        0.347
  26.1699     43       1   0.2625  0.0348       0.2025        0.340
  28.1425     42       1   0.2563  0.0345       0.1968        0.334
  29.9507     41       1   0.2500  0.0342       0.1912        0.327
  31.5616     40       1   0.2438  0.0339       0.1855        0.320
  33.1397     39       1   0.2375  0.0336       0.1799        0.314
  36.5918     38       1   0.2313  0.0333       0.1743        0.307
  43.3644     37       1   0.2250  0.0330       0.1688        0.300
  51.7479     36       1   0.2188  0.0327       0.1632        0.293
  55.2000     35       1   0.2125  0.0323       0.1577        0.286
  57.8301     34       2   0.2000  0.0316       0.1467        0.273
  58.4219     32       1   0.1938  0.0312       0.1412        0.266
  58.5534     31       1   0.1875  0.0309       0.1358        0.259
  59.2110     30       1   0.1813  0.0305       0.1304        0.252
  60.8877     28       1   0.1748  0.0300       0.1248        0.245
  61.6767     27       1   0.1683  0.0296       0.1192        0.238
  62.4329     26       1   0.1618  0.0292       0.1136        0.230
  66.6082     22       1   0.1545  0.0288       0.1072        0.223
  68.1205     21       1   0.1471  0.0283       0.1009        0.215
  79.0027     17       1   0.1385  0.0279       0.0932        0.206
  97.6438     10       1   0.1246  0.0284       0.0798        0.195
 112.6027      8       1   0.1090  0.0288       0.0650        0.183
 118.9808      7       1   0.0935  0.0286       0.0513        0.170
 125.6877      6       1   0.0779  0.0277       0.0388        0.157
Show the Code
ggsurvplot(KM1.age, data = surv_data, risk.table = TRUE,
           linetype = c(3,6), pval = TRUE)

Show the Code
Log_rank_age <- survdiff(Surv(time = surv_data$duration_in_month_guna_outcome, event = alive_0_death_1 == '1') ~ x60_0_60_1, data = surv_data, rho = 1)
Log_rank_age
Call:
survdiff(formula = Surv(time = surv_data$duration_in_month_guna_outcome, 
    event = alive_0_death_1 == "1") ~ x60_0_60_1, data = surv_data, 
    rho = 1)

               N Observed Expected (O-E)^2/E (O-E)^2/V
x60_0_60_1=0 288    124.7    145.7      3.03      13.6
x60_0_60_1=1 160     88.6     67.6      6.52      13.6

 Chisq= 13.6  on 1 degrees of freedom, p= 2e-04 
Show the Code
coxmodel1 <- coxph(Surv(time = surv_data$duration_in_month_guna_outcome, event = alive_0_death_1 == '1') ~ x60_0_60_1 + stage120_stage31satge42unkown3 + chemotherapy_0_no_yes_1_2_unknown, data = surv_data)

coxmodel1
Call:
coxph(formula = Surv(time = surv_data$duration_in_month_guna_outcome, 
    event = alive_0_death_1 == "1") ~ x60_0_60_1 + stage120_stage31satge42unkown3 + 
    chemotherapy_0_no_yes_1_2_unknown, data = surv_data)

                                       coef exp(coef) se(coef)      z        p
x60_0_60_11                         0.44014   1.55293  0.11391  3.864 0.000112
stage120_stage31satge42unkown31     0.04104   1.04189  0.28304  0.145 0.884725
stage120_stage31satge42unkown32     0.56182   1.75385  0.24336  2.309 0.020966
stage120_stage31satge42unkown33     0.10394   1.10954  0.22567  0.461 0.645085
chemotherapy_0_no_yes_1_2_unknown1 -0.39696   0.67236  0.12601 -3.150 0.001631
chemotherapy_0_no_yes_1_2_unknown2 -0.51271   0.59887  0.19811 -2.588 0.009653

Likelihood ratio test=40.18  on 6 df, p=4.193e-07
n= 448, number of events= 347 
Show the Code
coxmodel1.ia <- coxph(Surv(time = surv_data$duration_in_month_guna_outcome, event = alive_0_death_1 == '1') ~ x60_0_60_1 + stage120_stage31satge42unkown3 + chemotherapy_0_no_yes_1_2_unknown + x60_0_60_1:stage120_stage31satge42unkown3, data = surv_data)

coxmodel1.ia
Call:
coxph(formula = Surv(time = surv_data$duration_in_month_guna_outcome, 
    event = alive_0_death_1 == "1") ~ x60_0_60_1 + stage120_stage31satge42unkown3 + 
    chemotherapy_0_no_yes_1_2_unknown + x60_0_60_1:stage120_stage31satge42unkown3, 
    data = surv_data)

                                                coef exp(coef) se(coef)      z
x60_0_60_11                                 -0.06304   0.93891  0.50625 -0.125
stage120_stage31satge42unkown31              0.17113   1.18664  0.33407  0.512
stage120_stage31satge42unkown32              0.48621   1.62614  0.27928  1.741
stage120_stage31satge42unkown33             -0.10614   0.89929  0.25854 -0.411
chemotherapy_0_no_yes_1_2_unknown1          -0.41834   0.65814  0.12787 -3.272
chemotherapy_0_no_yes_1_2_unknown2          -0.51638   0.59668  0.19816 -2.606
x60_0_60_11:stage120_stage31satge42unkown31 -0.11426   0.89202  0.64194 -0.178
x60_0_60_11:stage120_stage31satge42unkown32  0.32571   1.38501  0.57378  0.568
x60_0_60_11:stage120_stage31satge42unkown33  0.65627   1.92759  0.52287  1.255
                                                  p
x60_0_60_11                                 0.90091
stage120_stage31satge42unkown31             0.60848
stage120_stage31satge42unkown32             0.08170
stage120_stage31satge42unkown33             0.68140
chemotherapy_0_no_yes_1_2_unknown1          0.00107
chemotherapy_0_no_yes_1_2_unknown2          0.00917
x60_0_60_11:stage120_stage31satge42unkown31 0.85873
x60_0_60_11:stage120_stage31satge42unkown32 0.57027
x60_0_60_11:stage120_stage31satge42unkown33 0.20943

Likelihood ratio test=45.55  on 9 df, p=7.306e-07
n= 448, number of events= 347 
Show the Code
prop.h <- cox.zph(coxmodel1, 
                  transform = 'km', 
                  global = TRUE)
plot(prop.h)

Show the Code
coxmodel1 %>% 
  tbl_regression(
    exponentiate = TRUE)
Characteristic HR1 95% CI1 p-value
x60_0_60_1


    0
    1 1.55 1.24, 1.94 <0.001
stage120_stage31satge42unkown3


    0
    1 1.04 0.60, 1.81 0.9
    2 1.75 1.09, 2.83 0.021
    3 1.11 0.71, 1.73 0.6
chemotherapy_0_no_yes_1_2_unknown


    0
    1 0.67 0.53, 0.86 0.002
    2 0.60 0.41, 0.88 0.010
1 HR = Hazard Ratio, CI = Confidence Interval