使用各种协变量迭代coxph（）模型

Question

我对R来说还是一个新手，所以这可能是一个基本问题。

我正在寻找使用coxph（）的co-cox模型的风险评估。我必须使用各种变量组合迭代模型大约60次。由于模型的每次迭代都会有不同的协变量（和主要曝光），我想写一个函数来完成它。在经过年龄调整的模型中，我只是主要曝光，一切运行良好。我可以添加协变量，它运行...我只需要一种方法来编写单个函数，其中“covars”可以是我放入函数调用中的任何内容。

注意：这是一个简化版本，它运行得很好，我只想让它工作而不写出60个独特的迭代。

   subtype <- function(expo, covars){


    temp <- coxph(Surv(FAIL, OUTCOME) ~ joint[[expo]]*strata(EVENT2)+
                  covars+
                  cluster(ID)+strata(AGE_INT),
                  na.action=na.exclude,
                  data=joint)                  

    return(summary(temp))
    }
    results <- subtype("RACE", covars=...)
    results2 <- subtype("GENDER", covers=...

当我在SAS中进行宏编程时，很容易。

感谢您的帮助。