## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(cox.rvph)
library(survival)

## ----step-example-------------------------------------------------------------
if (requireNamespace("KMsurv", quietly = TRUE)) {

  data("psych", package = "KMsurv")
  psych$sex <- factor(psych$sex)

  fit_step <- cox.rvph(
    data = psych,
    time = "time",
    event = "death",
    covariate = "age",
    adjust_vars = "sex",
    method = "step",
    verbose = FALSE
  )

  print(fit_step)
  summary(fit_step)
}

## ----timev-example------------------------------------------------------------
data("pbc", package = "survival")

pbc$status2 <- ifelse(pbc$status == 2, 1, 0)
pbc$ascites <- factor(pbc$ascites)

fit_timev <- cox.rvph(
  data = pbc,
  time = "time",
  event = "status2",
  covariate = "bili",
  adjust_vars = c("ascites", "edema", "protime"),
  method = "timev",
  verbose = FALSE
)

print(fit_timev)
summary(fit_timev)