validate(f, method="boot", B=1000, dxy=T)
rcorrcens(Surv(time, status) ~ predict(f), data = mydata)
f3 <- cph(Surv(time, status) ~ AGE + ELN + GRADE + SIZE + SUR + SEX , x=T, y=T, surv=T, data=mydata, time.inc=365)
cal3 <- calibrate(f3, cmethod="KM", method="boot", u=365, m=250 , B=1000)
plot(cal3)


雷达卡



京公网安备 11010802022788号







