生存分析绘制校准图

Dr-GAO

想做校准图，但是程序一致报错，
cox1 <- cph(Surv(survival_time,status) ~ age + sex + race + smoking + radiation + pharmaceutical + stage_T + stage_N + stage_M + surgery,surv=T,x=T, y=T,time.inc = 13655,data=tcga)

cal <- calibrate(cox1, cmethod="KM", method="boot", u=13655, m= 180, B=1000)

pdf("calibrate5.pdf",12,8)
par(mar = c(10,5,3,2),cex = 1.0)
plot(cal,lwd=3,lty=2,errbar.col="black"
,xlim = c(0,1),ylim = c(0,1)
,xlab ="Nomogram Predicted Survival"
,ylab="Actual Survival",col="blue")
lines(cal,c('mean.predicted','KM')
,type = 'a',lwd = 3,col ="black" ,pch = 16)
mtext(" ")
box(lwd = 1)
abline(0,1,lty = 3,lwd = 3,col = "black")
dev.off()

Dr-GAO

这个是报错的内容，请大佬指点

lines(cal,c('mean.predicted','KM')

,type = 'a',lwd = 3,col ="black" ,pch = 16)
Error in plot.xy(xy.coords(x, y), type = type, ...) :
invalid plot type 'a'

fenguoerbian

你自己可以再看一下报错具体说的是什么啊：lines这个函数里报错了，提示是type参数的取值，字符串“a”不合适。那么去看一下文档里这个type参数到底支持哪些取值呢？从文档里看默认是type="l"，同时文档里还提示了可以在plot的文档里看看其它的支持：

"p" for points, "l" for lines, "b" for both points and lines, "c" for empty points joined by lines, "o" for overplotted points and lines, "s" and "S" for stair steps and "h" for histogram-like vertical lines. Finally, "n" does not produce any points or lines.

Dr-GAO

fenguoerbian 好的非常感谢

chuxinyuan

Dr-GAO 请你规范发帖。