正确使用stats4 :: mle

Question

我想使用stats4 :: mle函数来估计分布的最佳参数（2）。

我想确保我的用法是正确的，并获得避免错误的指导

  

“ optim中的错误（开始，f，方法=方法，粗麻布= TRUE，...）：     'vmmin'中的初始值不是有限的   另外：警告消息：   在log（mu）中：产生的NaNs“

我想估计的函数是exp（beta0 * a + beta1 * b），我想估计betas

示例代码：

a <- mydata$a # first variable
b <- mydata$b # second variable
y <- mydata$y # observed result

nll <- function(beta0, beta1) {
mu = y - exp(beta0 * a + beta1 * b)
- sum(log(mu))
}

est <- stats4::mle(minuslog = nll, start = list(beta0 = 0.0001, beta1 = 0.0001))
est

所以：

• 这是正确的做事方式吗？
• 对于错误，我理解这是由于mu的值变为0，但我不知道该怎么办

感谢您的帮助。