MLE无法估计参数Weibull-Modified Generalized

Question

我正在尝试估计Weibull修改的广义分布参数，但是遇到以下错误：

“函数mle无法估计参数，错误代码为100”我该怎么办？

Weibull修正的广义分布

密度函数

   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
}

Harzard函数

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)