Skip to contents

Performs cure model analysis for interval-censored survival data using the ICGOR (Interval-Censored Generalized Odds Rate) methodology. This approach is specifically designed for situations where the exact event time is not known, but falls within an interval, and a fraction of the population may be cured (i.e., will never experience the event).

Usage

intervalcensorcure(
  data,
  leftTime,
  rightTime,
  covariates,
  survivalDistribution = "weibull",
  cureModel = "mixture",
  includeUncured = TRUE,
  confidenceLevel = 0.95,
  maxIterations = 1000,
  tolerance = 1e-06,
  bootstrapSamples = 500,
  performBootstrap = FALSE,
  modelComparison = TRUE,
  plotSurvival = TRUE,
  plotCure = TRUE
)

Arguments

data

the data as a data frame

leftTime

the left endpoint of the interval for interval-censored observations

rightTime

the right endpoint of the interval for interval-censored observations

covariates

covariates to include in the cure model

survivalDistribution

the parametric survival distribution for the uncured fraction

cureModel

the type of cure model to fit

includeUncured

whether to include detailed analysis of the uncured fraction

confidenceLevel

the confidence level for parameter estimates

maxIterations

maximum number of iterations for model convergence

tolerance

convergence tolerance for parameter estimation

bootstrapSamples

number of bootstrap samples for confidence intervals

performBootstrap

whether to perform bootstrap analysis for confidence intervals

modelComparison

whether to perform model comparison across different distributions

plotSurvival

whether to plot estimated survival functions

plotCure

whether to plot cure fraction estimates

Value

A results object containing:

results$instructionsa html
results$modelSummarya table
results$cureResultsa table
results$survivalResultsa table
results$modelFita table
results$modelComparisona table
results$bootstrapResultsa table
results$survivalPlotan image
results$curePlotan image
results$modelSummaryTexta html

Tables can be converted to data frames with asDF or as.data.frame. For example:

results$modelSummary$asDF

as.data.frame(results$modelSummary)

Details

Key features:

  • Mixture and non-mixture cure models

  • Various survival distributions (Weibull, exponential, log-normal, log-logistic)

  • Covariate effects on both cure fraction and survival parameters

  • Model comparison and goodness-of-fit assessment

  • Bootstrap confidence intervals

  • Parametric and non-parametric approaches

Examples

# Interval-censored cure model with Weibull distribution
intervalcensorcure(
    data = data,
    leftTime = "left_time",
    rightTime = "right_time",
    covariates = c("age", "treatment"),
    survivalDistribution = "weibull",
    cureModel = "mixture",
    includeUncured = TRUE
)