I am currently doing cox regression models with the coxph() function. The formula of my model is:
coxph()
cox.death <- coxph(Surv(diff,event)~Delay+ag