# ============================================================
# Cox Proportional Hazards
# Biostatistics in R — Rverse Analytics
# https://rverseanalytics.com/learn-biostatistics/cox
# ============================================================

# Cox proportional hazards
library(survival)
library(gtsummary)
# Fit a Cox model (lung cancer data)
cx <- coxph(Surv(time, status) ~ age + sex + ph.ecog, data = lung)
# Hazard ratios with 95% CI
tbl_regression(cx, exponentiate = TRUE)
# Check the proportional-hazards assumption
cox.zph(cx)