mle не удалось оценить параметры Weibull-Modified Generalized - PullRequest
Новая
       58

mle не удалось оценить параметры Weibull-Modified Generalized

0 голосов
/ 11 октября 2018

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

"функция 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)

Изображения

enter image description here enter image description here

Образец: 

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)
...