SurvTrunc (version 0.2.0)

coxDTdep: Fit Cox Proportional Hazards Regression Model Under Dependent Double Truncation

Description

Fits a Cox proportional hazards regression model under dependent double truncation. That is, when the survival time is subject to left truncation and/or right truncation and the survival times are dependent on at least one of the truncation times.

Usage

coxDTdep(
  formula,
  L,
  R,
  data,
  error = 1e-04,
  n.iter = 1000,
  n.boot = 100,
  CI.level = 0.95
)

Value

Displays the estimate, standard error, lower and upper bounds of confidence interval, Wald test statistic and p-value for each regression coefficient

Arguments

formula

a formula object, with the response on the left of a ~ operator, and the terms on the right. The response must be a survival object as returned by the Surv. function. NOTE: coxDTdep does not handle censoring.

L

vector of left truncation times. If only right truncation is present, set L = -infinity.

R

vector right truncation times. If only left truncation is present, set R=infinity.

data

mandatory data.frame matrix needed to interpret variables named in the formula

error

tolerance for convergence, default is 1e-4.

n.iter

maximun number of iterations for EM algorithm, default is 1000.

n.boot

number of bootstraps for computing standard errors, default is 100.

CI.level

a numeric value between 0.5 and 1 representing the confidence level for two-sided confidence intervals, default is 0.95

Details

Fits a Cox proportional hazards model in the presence of left, right, or double truncation when the survival times are dependent on at least one of the truncation times. An EM algorithm is employed to obtain point estimates for the regression coefficients. The standard errors are calculated using the bootstrap method. This method assumes no censoring is present in the data. Note: If only left truncation is present, set R=infinity. If only right truncation is present, set L = -infinity.

References

Rennert L and Xie SX (2022). Cox regression model under dependent truncation. Biometrics. 78(2), 460-473. https://doi.org/10.1111/biom.13451

Examples

###### Example: AIDS data set #####
if (FALSE) coxDTdep(Surv(Induction.time)~Adult, L=AIDS$L.time, R=AIDS$R.time, data=AIDS, n.boot=2)

# WARNING: To save computation time, number of bootstrap resamples for standard error set to 2.
# Note: The minimum recommendation is 100, which is the default setting.