我正试图在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时,我收到相同的错误消息。任何帮助将不胜感激!!
答案 0 :(得分:2)
您复制的示例显示了逻辑回归,但您正在拟合一个coxph模型,它们在处理方式上有很大不同。
如果您只想要一个风险比图,那么您的代码基本上可以正常工作(除非您添加到不存在的图,这可能是产生错误的图,尝试将add
更改为{ {1}})。
如果您想绘制生存曲线,然后使用FALSE
函数获取预测的生存信息并绘制出来。
答案 1 :(得分:0)
错误消息表明您没有打开设备,或者您尝试添加的图表可能还有其他问题?该代码使用我从coxph
帮助页面构建的玩具示例生成范围输入[0,1]上的绘图。也许你的'prov.yrs'的范围与现有的情节不同,或者没有设备打开?试试plot.new()
,绘制你要使用的其他内容,然后重新运行? (add = TRUE将禁止绘制框,轴和标签。)