在R中绘制cox比例风险模型

时间:2013-06-13 16:33:38

标签: r plot

我正试图在R中绘制一个cox比例风险模型(或logit模型) 我使用了以下代码(我从https://sites.google.com/site/daishizuka/toolkits/plotting-logistic-regression-in-r复制了)

c<-coxph(formula=Surv(year, promo)~prov.yrs, data=cul)
curve(predict(c, data.frame(prov.yrs=x), type="risk"), add=TRUE)

我收到错误消息

Error in plot.xy(xy.coords(x, y), type = type, ...) : 
  invalid graphics state

我认为绘制这个有什么问题,所以我想知道是否有办法来绘制这个。当我使用glm时,我收到相同的错误消息。任何帮助将不胜感激!!

2 个答案:

答案 0 :(得分:2)

您复制的示例显示了逻辑回归,但您正在拟合一个coxph模型,它们在处理方式上有很大不同。

如果您只想要一个风险比图,那么您的代码基本上可以正常工作(除非您添加到不存在的图,这可能是产生错误的图,尝试将add更改为{ {1}})。

如果您想绘制生存曲线,然后使用FALSE函数获取预测的生存信息并绘制出来。

答案 1 :(得分:0)

错误消息表明您没有打开设备,或者您尝试添加的图表可能还有其他问题?该代码使用我从coxph帮助页面构建的玩具示例生成范围输入[0,1]上的绘图。也许你的'prov.yrs'的范围与现有的情节不同,或者没有设备打开?试试plot.new(),绘制你要使用的其他内容,然后重新运行? (add = TRUE将禁止绘制框,轴和标签。)