simDC_weib {icenReg}

Simulate Doubly Censored Data

Description

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

Usage

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)