Как избежать проблемы вероятности журнала -Inf в функции mle из пакета stat4? - PullRequest
Новая
       41

Как избежать проблемы вероятности журнала -Inf в функции mle из пакета stat4?

0 голосов
/ 13 марта 2019

Я хочу максимизировать функцию логарифмического правдоподобия для соответствия некоторым данным, но функция mle останавливается с этой ошибкой, когда логарифмическое правдоподобие равно -Inf.

Ошибка в optim (start, f, method = method, hessian = TRUE, ...): L-BFGS-B требует конечных значений 'fn'

Код: 

  ll <- function(shape, scale) {
  x <-
    c(2,  12,  61,   9,  25, 110,   7,  27,  17,  15,  11,  17,  15,  25,   0)
  - sum(dgamma(x, shape, scale, log = TRUE))
  }

    m <- stats4::mle(ll, start = list(shape = 2.15, scale = 2), method = "L-BFGS-B", lower = c(0.01, 0.01))

Кто-нибудь знает, как избежать этой проблемы в чистом виде?

...