在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將禁止繪制框，軸和標簽。）