dxy logo
首页丁香园病例库全部版块
搜索
登录

【积分求助】关于R语言做calibration curve的相关问题

发布于 2015-05-13 · 浏览 3.2 万 · IP 山东山东
这个帖子发布于 10 年零 3 天前,其中的信息可能已发生改变或有所发展。
最近我想用R语言做诺谟图,参照版里的帖子已经做出了图,求出了C-index,但是做到校正曲线时总是出现问题,我参照的是肝胆医院发的一篇论文提到的一段校正曲线的代码:
cal < -calibrate(f, cmethod = ‘KM', method = “boot,” u = 3 or 5, m = 100, B = 1000)
par(mar = c,(8,5,3,2cex = 1.0)
plot(cal,lwd = 2,lty = 1,errbar.col = c(rgb(0,118,192,maxColorValue = 255)),xlim = c,(0,ylim = c,0,1xlab = ?,ylab = ?,col = c(rgb(192,98,83,maxColorValue = 255)))
lines(cal[,c('mean.predicted',‘KM')],type = ‘b',lwd = 2,col = c(rgb(192,98,83,maxColorValue = 255)),pch = 16)
mtext(“ ”)
box(lwd = 1)
abline(0,1,lty = 3,lwd = 2,col = c(rgb(0,118,192,maxColorValue = 255)))
我不是太明白u值指的是什么,如果是3年或5年生存率的话,因为我用的OS单位是月,所以,我改成了36和60,然后系统提示Using Cox survival estimates at 30 Days,为什么是30天,我并没有设置。而且图超丑,离45度线很远,我想知道是不是出了什么问题。
而且lines网上的代码似乎有点问题,我没有相关基础不知道该怎么修正啊,希望有人可以和我交流一下








最后编辑于 2022-10-09 · 浏览 3.2 万

66 82 7

全部讨论0

默认最新
avatar
66
分享帖子
share-weibo分享到微博
share-weibo分享到微信
认证
返回顶部