我正在构建生存模型,我的代码看起来像这样:
library('survival')
my.model <- coxhp(Surv(time, event) ~ var_1+var_2+var_3+var_4+var_5, data =df)
问题是我有太多变量,变量集总是在变化,我想知道是否可以创建变量列表并将其输入模型。像这样:
my.var <- c(var_1+var_2+var_3+var_4+var_5)
my.model <- coxhp(Surv(time, event) ~ my.var, data =df)
我找到了线性模型linear model solution的类似帖子,但不知道如何使它适应coxph。
答案 0 :(得分:2)
你可以使用“。”表示“所有未使用的变量”。所以对你的例子来说,
my.model <- coxph(Surv(time, event) ~ ., data = df[,c("time", "event", my.var)])
应该有用。