timereg (version 1.9.3)

cause.pchazard.sim: Simulation of Piecewise constant hazard models with two causes (Cox).

Description

Simulates data from piecwise constant baseline hazard that can also be of Cox type. Censor data at highest value of the break points.

Usage

cause.pchazard.sim(cumhaz1, cumhaz2, rr1, rr2, cens = NULL, rrc = NULL,
  cens.cum.hazard = TRUE, ...)

Arguments

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

rrc

retlativ risk for censoring.

cens.cum.hazard

possible cumulative hazard for censoring.

...

arguments for pc.hazard

Examples

# NOT RUN {
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 <-  cause.pchazard.sim(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)

# }