0

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

4

2 回答 2

2

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

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

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

于 2013-06-13T17:06:59.260 回答
0

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

于 2013-06-13T17:27:06.610 回答