COX风险分析

# 构建公式，包括混杂变量
  formula <- as.formula(paste0(
    "Surv(follow_up_time, MI_state == 1) ~ ", protein,
    " + age + sex + ethnicity + TDI + BMI + smoking_status + alcohol_intake+season"
  ))

上述代码是r语言中进行多变量cox风险回归分析的基础代码。研究的是蛋白水平/代谢水平对MI发生的影响，回归模型中加入了多个混杂变量（confounders），目的是控制可能影响 MI 发生的其他因素。