5

私の生存分析タスクでは、cox 比例モデルを使用して、データセットのさまざまなグループの一致指数 (c-index) 値を計算しました。c-index プロットに 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"))

前もって感謝します、

4

1 に答える 1