基于R的MLE定制分布的参数估计。

Question

我想使用mle（stats4）估算自定义分发的参数。 这个R代码是

    <div class="layout">
       <app-sidebar></app-sidebar>
         <!-- Page Container START -->
       <div class="page-container">
          <app-header></app-header>
          <router-outlet></router-outlet>
          <app-footer></app-footer>
       </div>
    <!-- Page Container END -->
   </div>

然后输出

library(stats4)
data=rweibull(1000, shape=5, scale = 1.5)
rn=data
n=length(rn)
r=max(rn)+0.05
ll= function(b,d){-(sum(log(r-rn)+n*log(d)-n*log(gamma(1+1/d)))-2*n*log(b)
         + (d-1)*sum(r-rn)-n*(d-1)*log(b)-sum((r-rn)^d)/(b^d))
    }
fit <- mle(ll, start = list(b=2,d=1.5), nobs = NROW(rn), method = "L-BFGS-B", lower = 1, upper = 100)
fit

现在的问题是，b的估计值每次都会随着数据（随机数）的变化而变化，但是d的估计值保持不变，并且给出的值I设置为L-BFGS-B的上限。 请为此提供帮助