# Logistic Regression
imp_dat <- df2milist(lbpmilr, impvar="Impnr")
res_stats <- with(data=imp_dat,
expr = cindex(glm(Chronic ~ Gender + Radiation,
family=binomial)))
res <- pool_cindex(res_stats)
res
# Cox regression
library(survival)
imp_dat <- df2milist(lbpmicox, impvar="Impnr")
res_stats <- with(data=imp_dat,
expr = cindex(coxph(Surv(Time, Status) ~ Pain + Radiation)))
res <- pool_cindex(res_stats)
res
Run the code above in your browser using DataLab