Значення попередження про конвергенцію в glmer

Я використовую glmerфункцію з lme4пакету в R, і я використовую bobyqaоптимізатор (тобто за замовчуванням у моєму випадку). Я отримую попередження, і мені цікаво, що це означає.

Warning message:
In optwrap(optimizer, devfun, start, rho$lower, control = control,  :
  convergence code 3 from bobyqa: bobyqa -- a trust region step failed to reduce q

Я шукав "крок області довіри не зміг q". Знайшов деяку інформацію в пакеті minqa , де сказано "Зверніться до Пауелла для пояснення". Я зробив (ви теж можете, якщо хочете! Див. Посилання та посилання на них нижче), але я не розумію. Насправді мені нічого не вдалося знайти про зменшення q.

MJD Powell (2007) "Розробки NEWUOA для необмеженої мінімізації без похідних", Кембриджський університет, кафедра прикладної математики та теоретичної фізики, Група чисельних аналізів, звіт NA2007 / 05, http://www.damtp.cam.ac.uk/ user / na / NA_papers / NA2007_05.pdf .

MJD Powell (2009), "Алгоритм BOBYQA для обмеженої обмеженої оптимізації без похідних", Доповідь № DAMTP 2009 / NA06, Центр математичних наук, Кембриджський університет, Великобританія. http://www.damtp.cam.ac.uk/user/na/NA_papers/NA2009_06.pdf .

Ps Я знаю, що можу змінити оптимізатор, і буду бачити, чи зможу отримати вихід без попереджень або помилок. Я також перевірю градієнт і Гессіана, якщо зможу, згідно коментаря / відповіді Бен Болкер . Я використовую glmerв dredgeвід MuMInі я не впевнений , що якщо відповідь Бен буде працювати без будь - якого додаткового майструвати, але я буду працювати над цим , як тільки моїм комп'ютером обробка , що він робить, у всякому разі, я відволікся.

Оновлення

Відповідно до коментаря доктора Болкера нижче, я почав переглядати код FORTRAN ( ось код для всіх, хто зацікавлений шукати, але не завантажувати його ). "430" з'являється в частині bobyqb.f коду. Просто знайдіть "430" або "зменшити Q", щоб знайти відповідний код.

Це моя перша зустріч з кодом FORTRAN, але я думаю, що код говорить про те, що якщо будуть виконані наступні умови, виробляйте попередження: NTRITS> 0, VQUAD> = 0, IPRINT> 0. "Ціле число NTRITS встановлюється на число" довірчий регіон ". ітерації, що відбулися з часу останньої "альтернативної" ітерації ". VQUADз'являється кілька разів, і я ще не зрозумів його значення, оскільки його значення, здається, залежить від багатьох інших змінних, значення яких іноді залежать від інших змінних. Від bobyqa.f: "Значення IPRINT повинно бути встановлено на 0, 1, 2 або 3, що контролює кількість друку. Зокрема, немає виводу, якщо IPRINT = 0, і є вихід лише у зворотному випадку, якщо IPRINT = 1. ".

Отже, здається, завдання полягає в тому, щоб розібратися в значенні того, що VQUADбути> = 0, і, можливо, зрозуміти, як / коли IPRINTстало> 0. Мені доведеться повернутися до папери, щоб подивитися, але математику, або на принаймні його символічне вираження, є для мене трохи бар’єром. Якщо хтось не знає про алгоритм або має бажання дізнатися про нього, я думаю, що мені доведеться ітеративно посилити своє розуміння попередження, повертаючись туди-сюди між паперами, кодом та Інтернетом, поки я не зрозумію, що це засоби.

$f(x)$$x_k$$k$

• $\Delta_k$
• $f(x)$$x = x_k$$Q(x)$
• $s_k$$Q_k(x_k + s_k)$$||s_k|| \leq \Delta_k$
• $s_k$$x_{k+1} = x_k + s_k$
• В іншому випадку вдосконаліть модель і спробуйте ще раз

$s_k$

Щоб зрозуміти значення VQUAD, спершу треба зрозуміти кілька інших змінних. На щастя, є хороші коментарі прямо під декларацією від SUBROUTINE BOBYQB. Найважливішими змінними є:

• GOPT, градієнт моделі
• HQ, гессіана моделі
• D$s_k$

Почавши кілька рядків вище 410, ви побачите DO 410 J=1,N. З цього починається цикл for-циклу (і вкладений для-циклу), який оцінює зміни, передбачені моделлю, використовуючи пробний крок D. Це накопичує передбачувані зміни в VQUAD. Перша частина for-loop оцінює умови першого порядку, а вкладена for-loop оцінює умови другого порядку. Напевно, було б простіше читати, якби петлі були відступними, як ось такий:

    DO 410 J=1,N
        VQUAD=VQUAD+D(J)*GOPT(J)
        DO 410 I=1,J
            IH=IH+1
            TEMP=D(I)*D(J)
            IF (I .EQ. J) TEMP=HALF*TEMP
410         VQUAD=VQUAD+HQ(IH)*TEMP

Після цього є ще одна for-loop для включення в модель інших параметрів. Маю визнати, я не повністю розумію цього - я найкраще здогадуюсь, що це стосується того, як вони будують модель.

Зрештою, все це VQUADмає на увазі зміну цільової функції, передбачену моделлю. Тож якщо VQUADце негатив, це погано. Тепер цей конкретний вирішувач може використовувати альтернативне обчислення кроків (ймовірно, пошук рядків), куди NTRITSпотрібно заграти. Отже, логіка в 430 говорить: "Якщо в останній ітерації використовувались альтернативний крок обчислення І модель не передбачає зменшення AND IPRINT> 0, роздрукуйте попереджувальне повідомлення". Зауважте, що вирішувач закінчується незалежно від значення IPRINT.

Якщо говорити про IPRINT, це значення передається BOBYQAфункцією виклику. У цьому випадку ваша R-процедура є функцією виклику. Є verboseпараметр до glmer- я би став кумером до доларів, до якого передається те саме значення BOBYQA. Спробуйте встановити verbose0, і ви, ймовірно, не побачите попередження. Але це, звичайно, не змінить того, що відбувається під капотом.