robin_surv: Covariate Adjusted and Stratified Survival Analysis

In RobinCar2: ROBust INference for Covariate Adjustment in Randomized Clinical Trials

Covariate Adjusted and Stratified Survival Analysis

Description

Calculate log-rank test as well as hazard ratio estimates for survival data, optionally adjusted for covariates and a stratification factor.

Usage

robin_surv(
  formula,
  data,
  treatment,
  comparisons,
  contrast = "hazardratio",
  test = "logrank",
  ...
)

Arguments

formula	(formula) A formula of analysis, of the form Surv(time, status) ~ covariates. (If no covariates should be adjusted for, use 1 instead on the right hand side.)
data	(data.frame) Input data frame.
treatment	(formula) A formula of treatment assignment or assignment by stratification.
comparisons	(list) An optional list of comparisons between treatment levels to be performed, see details. By default, all pairwise comparisons are performed automatically.
contrast	(character(1)) The contrast statistic to be used, currently only "hazardratio" is supported.
test	(character(1)) The test to be used, currently only "logrank" is supported.
...	Additional arguments passed to the survival analysis functions, in particular hr_se_plugin_adjusted (please see here for details).

Details

The user can optionally specify a list of comparisons between treatment levels to be performed. The list must have two elements:

• Treatment level indices of the treatment group.

• Treatment level indices of the control group.

So for example if you would like to compare level 3 with level 1, and also level 3 with level 2 (but not level 2 with level 1) then you can specify: comparisons = list(c(3, 3), c(1, 2))

Value

A surv_effect object containing the results of the survival analysis.

See Also

surv_effect_methods for S3 methods.

Examples

# Adjusted for covariates meal.cal and age and adjusted for stratification by sex:
robin_surv(
  formula = Surv(time, status) ~ meal.cal + age,
  data = surv_data,
  treatment = sex ~ strata
)

# Adjusted for stratification by strata but not for covariates:
robin_surv(
  formula = Surv(time, status) ~ 1,
  data = surv_data,
  treatment = sex ~ strata
)

# Unadjusted for covariates and stratification:
robin_surv(
  formula = Surv(time, status) ~ 1,
  data = surv_data,
  treatment = sex ~ 1
)

RobinCar2 documentation built on Sept. 9, 2025, 5:28 p.m.