найти MLE модели смеси в R

Question

Я пытаюсь найти MLE для модели экспоненциальной смеси, используя R. Логарифмическая вероятность составляет

Из моего понимания, во-первых, яопределена логарифмическая функция правдоподобия, такая что

LL <- function (p, alpha) {
    R = log(p) + log(alpha) - alpha * dat$dat0[del == 1] + 
        log(1 - p + p * exp(-alpha * dat$dat0[del == 0])) - sum(log(R))}

, где dat - это кадр данных dat0 (значения x) и del (коэффициент индикатора)

Изатем я применил mle function:

mle(
    LL, 
    start = list(p = 0.2, alpha = 0.2), 
    method = "L-BFGS-B", 
    lower = c(-Inf, 0), upper = c(Inf, Inf))

Но возникают ошибки, говорящие о том, что

Ошибка в optim (start, f, method = method, hessian = TRUE, ...):
L-BFGS-B нужны конечные значения 'fn'
Кроме того: Предупреждающие сообщения:
1: В журнале (p) + журнал (альфа) - альфа * dat $ dat [del== 1] + log (1 - p +:
длинна объекта не кратна короткой длине объекта
2: В журнале (R): произведено NaN

Я думаю,это может быть результатом неправильного выбора начальных параметров.

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

Спасибо.