在我的生存分析任务中,我使用 cox 比例模型来计算数据集不同组中的一致性指数 (c-index) 值。我想知道如何将 p 值添加到我的 c-index 图中以比较不同的组看起来像这个图?
这是我的代码:
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"))
提前致谢,