# NOT RUN {
data(TRACE)
dcut(TRACE) <- ~.
mm <- model.matrix(~-1+factor(wmicat.4),data=TRACE)
m1 <- gofM.phreg(Surv(time,status==9)~vf+chf+wmi,data=TRACE,modelmatrix=mm)
summary(m1)
par(mfrow=c(2,2))
plot(m1)
m1 <- gofM.phreg(Surv(time,status==9)~strata(vf)+chf+wmi,data=TRACE,modelmatrix=mm)
summary(m1)
# }
Run the code above in your browser using DataLab