R: Simulation of Piecewise constant hazard models with two...

rcrisk {timereg} | R Documentation

## 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

    rcrisk(cumhaz1, cumhaz2, rr1, rr2, cens = NULL, rrc = NULL, ...)

### 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 , rate or cumumlative hazard
`rrc` | retlativ risk for censoring.
`...` | arguments for rchaz

### Author(s)

Thomas Scheike

### Examples

    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)

[Package timereg version 1.9.4 Index]