Displaying 1 result from an estimated 1 matches for "ddist1".
Did you mean:
ddist
2005 Jul 11
1
validation, calibration and Design
...ferent to an existing model
based on 60 month survival data generated by it (eg by McNemar's test).
Hence, somewhat shortened:
#using R 2.01 on Windows
library(Hmisc)
library(Design)
data1 #dataframe with predictor variables A and B, cens and time
columns (months)
ddist1 <- datadist(data1)
options(datadist='ddist1')
s1 <- Surv(data1$time, data1$cens)
cph.nomo <- cph(s1 ~ A+B, surv=T, x=T, y=T, time.inc=60)
survcph <- Survival(cph.nomo, x=T, y=T, time.inc=60, surv=T)
surv5 <- function(lp) survcph(60, lp)
nomogram(cph.nomo, lp=...