Я хочу использовать функцию stats4 :: mle для оценки лучших параметров (2) распределения.
Я хотел бы убедиться, что мое использование правильное, и получить рекомендации, чтобы избежать ошибки
"Ошибка в optim (start, f, method = method, hessian = TRUE, ...): начальное значение в 'vmmin' не является конечным. Дополнительно: Предупреждающее сообщение: В журнале (mu): произведено NaN"
Функция, которую я хотел бы оценить, это exp (beta0 * a + beta1 * b), и я хотел бы оценить беты
Пример кода:
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, но я нене знаю, что я могу с этим сделать
Спасибо за вашу помощь.