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

Question

我正在尝试在R中绘制一个cox比例风险模型（或logit模型）我使用了以下代码（我从https://sites.google.com/site/daishizuka/toolkits/plotting-logistic-复制了该代码回归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时，我收到相同的错误消息。 任何帮助将不胜感激！！

Answer 1

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

如果您只想要一个风险比图，那么您的代码基本上可以工作（除非您添加到不存在的图，这可能是产生错误的图，尝试将add更改为FALSE ）。

如果你想绘制生存曲线，那么使用survfit函数来获得预测的生存信息并绘制出来。

Answer 2

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