Ошибка MLE в R: ошибка в optim (start, f, method = method, hessian = TRUE, ...): не конечное разностное значение [1]

Question

Я работаю над подбором модели, используя MLE в R. Я оценил параметры в excel, но когда я попытался воспроизвести тот же результат в R., появляется ошибка.

loglikegomp <- function(b,c)
{
  mu <- b*c^gompdata$age
  u = dpois(gompdata$deaths,gompdata$exposure*mu,log=TRUE)
  -sum(u)
}
loglikegomp()

fittedgomp <- mle(loglikegomp, start = list(b=0.00004,c=1.06),method="BFGS")

Ошибка было

Error in optim(start, f, method = method, hessian = TRUE, ...) : 
non-finite finite-difference value [1]
In addition: Warning message:
In dpois(gompdata$deaths, gompdata$exposure * mu, log = TRUE) :
NaNs produced

Я пытался изменить начальные значения, но это не помогло.