如何通过函数输入向R语言Cox模型中添加协变量

报错原因 #

R的公式语法~不会自动解析字符串变量的内容，你代码中写Surv(time, status) ~ covariates时，R会默认到数据集df中寻找名为covariates的列，找不到就会代入全局环境中存储的字符串对象，字符串长度和数据集行数不匹配，就会触发长度不一致的报错。
你之前使用的noquote()仅会改变字符串的打印效果，不会将字符串转换为公式可识别的表达式，因此无法解决问题。

解决方法 #

有两种常用的实现方式，都可以实现协变量的动态传入：

方法1：使用reformulate()直接生成公式（推荐） #

reformulate是R专门用来构建公式的内置函数，写法更简洁，不易出错：

# 封装为函数的示例
fit_cox <- function(covariates, df) {
  # 第一个参数传协变量向量，第二个参数传响应变量
  cox_formula <- reformulate(termlabels = covariates, response = "Surv(time, status)")
  cntrl_reduced <- coxph(cox_formula, data = df)
  return(cntrl_reduced)
}

# 调用测试
covariates = c('A', 'B', 'C')
model <- fit_cox(covariates, df)

方法2：手动拼接字符串后转公式 #

如果要沿用你之前拼接协变量的思路，只需要把拼接后的完整公式字符串转换为formula类型即可：

fit_cox <- function(covariates, df) {
  # 拼接协变量部分
  cov_str <- paste(covariates, collapse = ' + ')
  # 拼接完整公式字符串后转为formula类型
  cox_formula <- as.formula(paste0("Surv(time, status) ~ ", cov_str))
  cntrl_reduced <- coxph(cox_formula, data = df)
  return(cntrl_reduced)
}

内容的提问来源于stack exchange，提问作者JVDeasyas123