survivalROC包比较两个AUC的P值
时间: 2023-07-15 08:12:00 浏览: 207
survivalROC包提供了两种方法来比较两个AUC的P值:
1. 根据DeLong方法,使用`pROC::roc.test()`函数。这个函数可以计算两个AUC之间的差异,并返回一个p-value。例如:
```
library(survivalROC)
library(survival)
data(veteran)
fit1 <- survfit(Surv(time, status) ~ trt, data = veteran)
fit2 <- survfit(Surv(time, status) ~ celltype, data = veteran)
auc1 <- survAUC(fit1, method = "simple")
auc2 <- survAUC(fit2, method = "simple")
p_value <- pROC::roc.test(auc1$AUC, auc2$AUC)$p.value
```
2. 根据比例风险模型,使用`survivalROC::coxtest()`函数。这个函数可以比较两个组之间的比例风险,并返回一个p-value。例如:
```
library(survivalROC)
library(survival)
data(veteran)
fit1 <- survfit(Surv(time, status) ~ trt, data = veteran)
fit2 <- survfit(Surv(time, status) ~ celltype, data = veteran)
p_value <- coxtest(fit1, fit2)$p.val
```
无论使用哪种方法,都需要先计算出两个组的AUC值,然后将它们传递给相应的函数进行比较。
阅读全文