R生存分析：使用Weibull的survreg錯誤

Question

我正在嘗試使用以前在SAS中生成的Weibull分布來復制一些生存分析-我現在使用的是無執照的計算機，因此使用R（均來自Windows）。 我的（經過審查）輸入數據如下所示：

> head(mydata)
  ID         Key  Time  Score    Event    Censor
1 1231231    ZXC   28   182.34   0      1
2 4564564    ASD   28   320.04   0      1
3 7897897    QWE   28   306.32   0      1
4 9879879    QWE   28   211.92   0      1
5 6546546    ASD   28   276.14   0      1
6 3213213    ZXC   28   331.50   0      1

“事件”和“檢查器”是二進制文件，“分數”在150到450之間變化，“時間”在1到28之間變化。輸入數據集中大約有30,000行。

當我嘗試：

mydatasr <- survreg(Surv(Time, Censor) ~ Score, dist = "w")

我收到警告消息：

在survreg.fit（X，Y，weights，offset，init = init，controlvals = control，中：跑出迭代並且沒有收斂，

而且沒有輸出。

我已經在線（並且通過本網站）搜索了此味精，但沒有找到任何可以指出問題所在的信息。 我沒有收斂問題，可以通過SAS中的（過程邏輯和）lifereg來放入相同的數據。

Answer 1

在沒有數據的情況下很難知道。 您可以將迭代次數加倍（或如下圖所示，加倍），默認為30：

(mydatasr <- survreg(Surv(Time, Censor) ~ Score , dist = "w", control = list(maxiter=90) )

有關更多選項，請參見？survreg.control。 我還猜想您可能已經錯過了Surv對象在公式之前有一個右括號的~

Answer 2

嘗試這個：

   survreg(Surv(Time, Censor) ~ Score, data=mydata, dist = "w", scale=1)