# 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)
## cumulative sums in covariates, via design matrix mm
mm <- cumContr(TRACE$wmi,breaks=10,equi=TRUE)
m1 <- gofM.phreg(Surv(time,status==9)~strata(vf)+chf+wmi,data=TRACE,
modelmatrix=mm,silent=0)
summary(m1)
## cumulative sums in covariates, via design matrix mm
m1 <- gofZ.phreg(Surv(time,status==9)~strata(vf)+chf+wmi,data=TRACE,vars="wmi")
summary(m1)
# }
Run the code above in your browser using DataLab