Extracts adjusted hazard ratios (HR) for a categorical variable from a multivariable Cox proportional hazards model. Performs LRT for variables with >2 categories.

cox_multivar_fun(
  data_fun,
  res_time,
  res_event,
  var_model,
  var,
  multi_var_all = multi_var_all
)

Arguments

data_fun

Data frame containing survival data.

res_time

Character. Name of survival time variable.

res_event

Character. Name of event indicator (1 = event, 0 = censored).

var_model

Character vector of variables included in the multivariable model.

var

Character. Categorical predictor variable.

multi_var_all

Character vector of full model candidate variables. Used to determine whether var is included in var_model.

Value

A matrix of adjusted HRs and p-values.

Examples

if (FALSE) { # \dontrun{
cox_multivar_fun(mydata, "OS_time", "OS_event",
                 var_model = c("Stage","Age"),
                 var = "Stage")
} # }