Я пытаюсь воспроизвести некоторый анализ выживания, используя дистрибутив Weibull, который я ранее создал в SAS - теперь я работаю с нелицензионной машины, поэтому использую 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
где Event и Censor являются двоичными, Score варьируется от 150 до 450, а Time — от 1 до 28. Во входном наборе данных около 30 000 строк.
Когда я пытаюсь:
mydatasr <- survreg(Surv(Time, Censor) ~ Score, dist = "w")
Я получаю предупреждающее сообщение:
В survreg.fit(X, Y, weights, offset, init = init, controlvals = control, : исчерпаны итерации и не сходятся,
И никакого выхода.
Я искал это сообщение в Интернете (и на этом сайте), но мне не удалось найти ничего, что указывало бы на то, в чем может быть проблема. У меня не было проблем с конвергенцией при передаче одних и тех же данных через (процесс логистики и) lifereg в SAS.
print.survreg
. Присваивание не обязательно приводит к выполнениюprint
-методов. - person IRTFM   schedule 25.08.2015