simDC_weib function

Simulate Doubly Censored Data

Simulates doubly censored data from a survival regression model with a Weibull baseline distribution.

simDC_weib(
  n = 100,
  b1 = 0.5,
  b2 = -0.5,
  model = "ph",
  shape = 2,
  scale = 2,
  lowerLimit = 0.75,
  upperLimit = 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
• lowerLimit: Lower censoring threshold
• upperLimit: Upper censoring threshold

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, all values less than lowerLimit are left censored and all values less than upperLimit are right censored.

Examples

simData <- simDC_weib()
fit <- ic_par(cbind(l, u) ~ x1 + x2, data = simData)