powered by
mycph(formula, data, digits = 2)
require(survival) data(colon) attach(colon) colon$TS=Surv(time,status==1) out=mycph(TS~.,data=colon) out HRplot(out,type=2,show.CI=TRUE,main="Hazard ratios of all individual variables")
Run the code above in your browser using DataLab