Оценка параметров плотности Weibull (ошибка: «... начальное значение в« VMMIN »не конечно»)

Question

Я пытаюсь оценить форму и масштаб набора данных.
Я использовал два разных способа, и для обоих я получил сообщение об ошибке:

Во-первых, я пытался по моментам, используя survey упаковка:

survreg(Surv(all.ws)~1, dist="weibull")

Я получил сообщение об ошибке:

Неверное время выживания для этого распределения

Во-вторых, я пытался использовать fitdistr() Функция:

fitdistr(all.ws, densfun=dweibull, start=list(scale=1, shape=2))    

Я получил сообщение об ошибке:

Ошибка в OPTION (X = C (2.2, 2.1,1,9 ....:
Первоначальное значение в «VMMIN» не конечно

Что не так с чем я делаю?

Answer 1

Поиск Google "Fitdistr Weibull Error" показывает, что этот точный вопрос обсуждался год назад в списке рассылки R-HELT: http://r.789695.n4.nabilb.com/problems-with-fitdistr-td1334772.html.

Некоторые очки из этой ссылки:

• Zeros В ваших данных приведет к проблемам

• Используйте функцию Pelwei () из пакета LMOM

Answer 2

У меня была похожая проблема при использовании FITDDISTR () с бета-распределением. В этом случае оба одинаково и Zeros в данных дали эту ошибку.

Кроме того, я обнаружил, что когда предельный аргумент используется в вызове FITDDISTR (), производится другое ошибка: например, (где x - вектор образцов, содержащий 1,0 или 0):

fitdistr(x, "beta", list(shape1 = 1, shape2 = 0.2),lower=0.001)

   `Error in stats::optim: L-BFGS-B needs finite values of 'fn'`