Warning: file_get_contents(/data/phpspider/zhask/data//catemap/8/perl/11.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 单变量coxph函数在多个协变量中的应用_R_Cox Regression - Fatal编程技术网

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

R 单变量coxph函数在多个协变量中的应用,r,cox-regression,R,Cox Regression,对于我的数据集列,我有时间、状态、治疗、性别、BMI。我试图应用单变量coxph函数,并结合协变量治疗、性别和BMI 这是我的代码: univ_models <- lapply( univ_formulas, function(x){coxph(x, data = RemissionTimes)}) univ_results <- lapply(univ_models, function(x){

对于我的数据集列,我有时间、状态、治疗、性别、BMI。我试图应用单变量coxph函数,并结合协变量治疗、性别和BMI

这是我的代码:

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)

我不知道如何修理它。我看不到行的不同之处,也不明白为什么R找不到res。

对象res是否出现在前面的代码中?您还可以检查它们是否是表中的空白。我不知道这是否对您有帮助。

前面的代码中是否显示对象res?您还可以检查它们是否是表中的空白。我不知道这是否对您有帮助。

res以前没有出现在我的代码中,我的数据集中有0,但没有空格。更多的帮助我将不胜感激。在我的数据集中,我有一列a,B,C。这些都是因数,但这会在处理中增加更多的行。如何删除这一额外的行?对于您的第一个问题,我认为0仍然可以包含在生存数据中。你说res以前没有出现在代码中。试着用这种方式来拟合:>#单变量Cox分析的计算方法如下:#>res.Cox res.Cox然后用你发布的代码进行跟踪。至于列,我不确定我是否得到了clearlyres,它以前没有出现在我的代码中,在我的数据集中有0,但没有空格。更多的帮助我将不胜感激。在我的数据集中,我有一列a,B,C。这些都是因数,但这会在处理中增加更多的行。如何删除这一额外的行?对于您的第一个问题,我认为0仍然可以包含在生存数据中。你说res以前没有出现在代码中。试着用这种方式来拟合:>#单变量Cox分析的计算方法如下:#>res.Cox res.Cox然后用你发布的代码进行跟踪。至于专栏,我不知道我是否明白你的意思
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