24小时热门版块排行榜    

查看: 8298  |  回复: 3

lilylorry

金虫 (初入文坛)

[求助] R软件如何通过bootstrap方法得到ROC曲线下面积(AUC)的95%CI?

请教各位大侠如何通过bootstrap方法得到AUC的95%CI及其p值。
现找到一段程序,但是看不懂:
# Bootstrap 95% CI for R-Squared
library(boot)
# function to obtain R-Squared from the data
rsq <- function(formula, data, indices) {
d <- data[indices,] # allows boot to select sample
fit <- lm(formula, data=d)
return(summary(fit)$r.square)
}
# bootstrapping with 1000 replications
results <- boot(data=mtcars, statistic=rsq,
R=1000, formula=mpg~wt+disp)
我的程序中已作出AUC趋势图,但是得不到AUC的95%CI及p值,模仿上面编程出不来结果。
我的程序如下(无AUC的95%CI及p值):
fit<-coxph(Surv(survival.time,survival.status)~IPI+β2m,na.action=na.omit)
eta0<-fit$linear.predictor
AUC.CC=risksetAUC(Stime=survival.time,
status=survival.status, marker=eta, method="Cox", tmax=40)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冰苹果

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
这个问题可以去COS论坛里面问问。那里面有很多R语言的高手!希望对你有帮助

» 本帖已获得的红花(最新10朵)

2楼2012-02-27 10:28:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilylorry

金虫 (初入文坛)

送鲜花一朵
引用回帖:
: Originally posted by 冰苹果 at 2012-02-27 10:28:49:
这个问题可以去COS论坛里面问问。那里面有很多R语言的高手!希望对你有帮助

我已经在那个里面求助过,没有响应,不过还是谢谢你!
3楼2012-02-27 22:10:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zz97512468

铁虫 (初入文坛)

只会用别的软件计算额~
跟随那一缕阳光,就能找到最美的世界
4楼2012-05-14 15:25:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lilylorry 的主题更新
信息提示
请填处理意见