Не удалось оценить параметры Вейбулла-Модифицированный Обобщенный

Question

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

"Функция mle не смогла оценить параметры, с кодом ошибки 100" Что мне делать?

Модифицированное Вейбуллом обобщенное распределение

Функция плотности

   dgmw = function(x, beta,theta, gama,lambda)
{
  ((beta)*(theta)*x^(gama-1))*(gama+lambda*x)*exp((lambda*x)-(beta*x^(gama))*exp(lambda*x))*
    (1-exp(-beta*x^(gama)*exp(lambda*x)))^theta-1

}

Накопительная функция

pgmw = function(x, beta,theta, gama,lambda)
{
  (1-exp(-beta*x^(gama)*exp(lambda*x)))^theta
}

Функция Гарзарда

hgmw = function(x, beta,theta, gama,lambda)
{
  (((beta)*(theta)*x^(gama-1))*(gama+lambda*x)*exp((lambda*x)-(beta*x^(gama))*exp(lambda*x))*
     (1-exp(-beta*x^(gama)*exp(lambda*x)))^theta-1)/1-(1-exp(-beta*x^(gama)*exp(lambda*x)))^theta
}

Функция выживания

sgmw = function(x, beta,theta, gama,lambda)
{
  1-(1-exp(-beta*x^(gama)*exp(lambda*x)))^theta
}

Предварительный расчет

    paramMGW = fitdist(data = dadosp, distr = 'gmw', start = c(2,0.5,2,0.3),lower = c(0, 0))
summary(paramMGW)

Изображений

Образец:

dadosp = c(240.3,71.9,271.3, 186.3,241,253,287.4,138.3,206.9,176,270.4,73.3,118.9,203.1,139.7,31,269.6,140.2,205.1,133.2,107,354.6,277,27.6,186,260.9,350.4,242.6,292.5, 112.3,242.8,310.7,309.9,53.1,326.5,145.7,271.5, 117.5,264.7,243.9,182,136.7,103.8,188.3,236,419.8,338.6,357.7)