gistfile1.txt

pancreas2<-pancreas%>%dplyr::filter(!Patientnr%in%8 & !Patientnr%in%20) 
pancreas2$meds <- c(0,0,1,0,1,0,0,0,0,1,0,0,0,1,0,0,0,0,1,1,0,0,1,0,0,0,0,0)  
pancreas2$RT <- car::recode(pancreas2$Total_dose_spec, "1:2=1; 3=2")
table(pancreas2$RT)
class(pancreas2$RT)
fitSURV<-coxph(Surv(SurvTimeWeek, Event) ~ 
                 Gender+
                 NRS_WorstPain_at_consultation+
                 Disease_stage_at_consultation+
                 Age_at_consultation+
                 Tumorsite_in_pancreas+
                 Non_opioid_spec+
                 Opioid_spec+
                 Surgery_prior_to_consultation+
                 Total_dose_spec,
                 data = pancreas2, x = TRUE, 
                 na.action=na.exclude, 
                 cluster = Patientnr)
summary(fitSURV)