r - 如何在R中的一致性指数图上添加p值？

在我的生存分析任务中，我使用了Cox比例模型来计算数据集不同组中的一致性指标(c-index)值。
我想知道如何在我的c指数图中添加p值，以比较不同的组，使它看起来像这个数字吗？

这是我的代码:

 surv <- with(group, Surv(group$survival, group$time))

# calculate survival
 sum.surv_1 <-  with(group, summary(coxph(surv ~ group$1)))
 sum.surv.1_2 <-  with(group, summary(coxph(surv ~ group$1 + group$2,ties = T)))


c_index.1 <- sum.surv_1$concordance
c_index.1_2 <- sum.surv.1_2$concordance

Comb_cIndex = data.frame(rbind(c_index.1["concordance.concordant"],
                           c_index.1_2["concordance.concordant"]))

barplot(as.matrix(Comb_cIndex), beside=TRUE, axis.lty=1,
    ylab = "C Index", ylim = c(0, 0.8),
    col = c("green", "blue"))

提前致谢，

最佳答案

我设法通过为数据集计算有/无关系的一致性值来找到答案。
通过遵循以下示例:

测试两个忽略索引的c索引的p值

round(cindex.comp(c_index_no_ties1, c_index_no_ties2)$p.value,3)

测试两个c指数的领带的p值的函数
因变量的t检验用于显着性
输入变量是从第一个函数获得的对象

cindex.p.ties <- function(c_index_ties1, c_index_ties2, c_index_no_ties1, c_index_no_ties2) {
    eps <- 1E-15
    n <- c_index_no_ties1$n
    r <- cor(c_index_no_ties1$data$x, c_index_no_ties2$data$x, use="complete.obs", method="spearman")
    if ((1 - abs(r)) > eps) {
      t.stat <- (c_index_ties1$concordance - c_index_ties2$concordance) / sqrt(c_index_ties1$std.err^2 + c_index_ties2$std.err^2 - 2 * r * c_index_ties1$std.err * c_index_ties2$std.err)
      diff.ci.p <- pt(q=t.stat, df=n - 1, lower.tail=FALSE)
    } else { diff.ci.p <- 1 }
    return(list("p.value"=diff.ci.p))
  }
cindex.p.ties(c_index_ties1=c_index_ties1, c_index_ties2=c_index_ties2, c_index_no_ties1=c_index_no_ties1, c_index_no_ties2=c_index_no_ties2)

关于r - 如何在R中的一致性指数图上添加p值？，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/30279314/