coxph.detail(fit)
fit <- coxph(Surv(futime, fustat) ~ age + rx + ecog.ps, fleming, x=T) fitd <- coxph.details(fit) events <- fit$y[,2]==1 etime <- fit$y[events,1] # the event times -- may have duplicates indx <- match(etime, fitd$time) sresid <- fit$x[events,] - fitd$means[indx,]