library(survival)
# The simplest test data set from coxph function
test1 <- list(time=c(4,3,1,1,2,2,3),
status=c(1,1,1,0,1,1,0),
x=c(0,2,1,1,1,0,0),
sex=c(0,0,0,0,1,1,1))
AIC(coxph(Surv(time,status)~x + strata(sex), data=test1))
modelexpr <- "cox(time,status)/strata(sex)~exp(beta*x)"
AIC(epifit(modelexpr=modelexpr, data=test1))
Run the code above in your browser using DataLab