单变量coxph函数在多变量中的应用

Question

对于我的数据集列，我有时间、状态、治疗、性别、BMI。我正在尝试应用单变量coxph函数与协变量治疗，性别和体重指数。

这是我的代码：

univ_models <- lapply( univ_formulas, function(x){coxph(x, data = RemissionTimes)})
univ_results <- lapply(univ_models,
                       function(x){ 
                           x <- summary(x)
                           p.value<-signif(x$wald["pvalue"], digits=2)
                           wald.test<-signif(x$wald["test"], digits=2)
                           beta<-signif(x$coef[1], digits=2);#coeficient beta
                           HR <-signif(x$coef[2], digits=2);#exp(beta)
                           HR.confint.lower <- signif(x$conf.int[,"lower .95"], 2)
                           HR.confint.upper <- signif(x$conf.int[,"upper .95"],2)
                           HR <- paste0(HR, " (", 
                                        HR.confint.lower, "-", HR.confint.upper, ")")
                           res<-c(beta, HR, wald.test, p.value)
                           names(res)<-c("beta", "HR (95% CI for HR)", "wald.test", 
                                         "p.value")
                           return(res)
                           #return(exp(cbind(coef(x),confint(x))))
                       })
res <- t(as.data.frame(univ_results, check.names = FALSE))
as.data.frame(res)

但是，当我运行它时会得到这样的结果：

Error in (function (..., row.names = NULL, check.rows = FALSE, check.names = TRUE,  : 
  arguments imply differing number of rows: 5, 4
> as.data.frame(res)
Error in as.data.frame(res) : object 'res' not found

我不知道该怎么解决它。我看不出行的不同之处，我也不明白为什么R找不到res。

Answer 1

对象res是否出现在前面的代码中？还可以检查表中的空格是否为空格。我不知道这是否会对你有所帮助。