data(TRACE)
dcut(TRACE) <- ~.
out1 <- phreg(Surv(time,status==9)~vf+chf+strata(wmicat.4),data=TRACE)
## robust standard errors default
summary(out1)
par(mfrow=c(1,2))
plot(out1)
## computing robust variance for baseline
rob1 <- robust.phreg(out1)
plot(rob1,se=TRUE,robust=TRUE)
## taking iid decomposition of regression parameters
betaiiid <- lava::iid(out1)
## making iid decomposition of baseline at a specific time-point
Aiiid <- mets::IIDbaseline.phreg(out1,time=30)
Run the code above in your browser using DataLab