# rcrisk: Simulation of Piecewise constant hazard models with two... In timereg: Flexible Regression Models for Survival Data

 rcrisk 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 for either of the cumulatives.

### Usage

```rcrisk(cumhaz1, cumhaz2, rr1, rr2, n = NULL, 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. `n` number of simulation if rr not given `cens` to censor further , rate or cumumlative hazard `rrc` retlativ risk for censoring. `...` arguments for rchaz

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)

```

timereg documentation built on Nov. 10, 2022, 5:48 p.m.