# NOT RUN {
library(rms)
# logistic
data(lung)
lung$status=lung$status-1
f1=lrm(status~age+sex,lung)
f2=lrm(status~age+sex+ph.ecog,lung)
idi(f1,f2)
# survival
head(lung)
range(lung$time)
lung=lung[complete.cases(lung),]
f1=cph(Surv(time,status)~age+sex,lung)
f2=cph(Surv(time,status)~age+sex+ph.ecog,lung)
idi(f1,f2)
# }
Run the code above in your browser using DataLab