Question

我正试图在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时，我收到相同的错误消息。任何帮助将不胜感激!!

Answer 1

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

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

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

Answer 2

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

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

2 个答案: