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

【问题标题】：plotting cox proportional hazard model in R在R中绘制cox比例风险模型
【发布时间】：2013-06-13 16:33:37
【问题描述】：

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

【问题讨论】：

为我运行示例... R 3.0.1 W7 64 位。

标签： r plot

【解决方案1】：

您复制的示例显示了逻辑回归，但您正在拟合 coxph 模型，它们的处理方式非常不同。

如果您只想绘制风险比图，那么您的代码基本上可以工作（除非您要添加到不存在的图，这可能是导致错误的原因，请尝试将 add 更改为 FALSE )。

如果您想绘制生存曲线，请使用survfit 函数获取预测的生存信息并绘制它。

【讨论】：

【解决方案2】：

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

【讨论】：