coxph(formula = Surv(tstart, tstop, event) ~ Treatment + cluster(ID), data=data))