rcrisk {timereg} | R Documentation |
Simulates data from piecwise constant baseline hazard that can also be of Cox type. Censor data at highest value of the break points.
rcrisk(cumhaz1, cumhaz2, rr1, rr2, cens = NULL, rrc = NULL, ...)
cumhaz1 |
cumulative hazard of cause 1 |
cumhaz2 |
cumulative hazard of cause 1 |
rr1 |
number of simulations or vector of relative risk for simuations. |
rr2 |
number of simulations or vector of relative risk for simuations. |
cens |
to censor further , rate or cumumlative hazard |
rrc |
retlativ risk for censoring. |
... |
arguments for rchaz |
Thomas Scheike
data(TRACE) cox1 <- cox.aalen(Surv(time,status==9)~prop(vf)+prop(chf)+prop(wmi), data=TRACE,robust=0) cox2 <- cox.aalen(Surv(time,status==0)~prop(vf)+prop(chf)+prop(wmi), data=TRACE,robust=0) X1 <- TRACE[,c("vf","chf","wmi")] n <- 1000 xid <- sample(1:nrow(X1),n,replace=TRUE) Z1 <- X1[xid,] Z2 <- X1[xid,] rr1 <- exp(as.matrix(Z1) %*% cox1$gamma) rr2 <- exp(as.matrix(Z2) %*% cox2$gamma) cumhaz1 <- cox1$cum cumhaz2 <- cox2$cum d <- rcrisk(cox1$cum,cox2$cum,rr1,rr2) dd <- cbind(d,Z1) sc1 <- cox.aalen(Surv(time,status==1)~prop(vf)+prop(chf)+prop(wmi), data=dd,robust=0) cbind(sc1$gamma, cox1$gamma) sc2 <- cox.aalen(Surv(time,status==2)~prop(vf)+prop(chf)+prop(wmi), data=dd,robust=0) cbind(sc2$gamma, cox2$gamma) par(mfrow=c(1,2)) plot(cox1); lines(sc1$cum,col=2) plot(cox2$cum,type="l"); lines(sc2$cum,col=2)