simCS_weib function

Simulate Current Status Data

Simulates current status data from a survival regression model with a Weibull baseline distribution.

simCS_weib(n = 100, b1 = 0.5, b2 = -0.5, model = "ph", shape = 2, scale = 2)

Arguments

• n: Number of observations
• b1: Regression coefficient 1
• b2: Regression coefficient 2
• model: Regression model to use. Choices are "ph", "po" or "aft"
• shape: Baseline shape parameter
• scale: Baseline scale parameter

Details

Exact event times are simulated according to the given survival regression model. Two covariates are used; x1 = rnorm(n), x2 = 1 - 2 * rbinom(n, 1, .5). After event times are simulated, current status inspection times are simulated following the exact same conditional distribution as event time (so each event time necessarily has probability 0.5 of being right censored).

Returns data in current status format, i.e. inspection time and event indicator. Use cs2ic to convert to interval censored format (see example).

Examples

simData <- simCS_weib()
fit <- ic_par(cs2ic(time, event) ~ x1 + x2, data = simData)