Я использую glmer
в R для запуска биномиальной модели для прогнозирования преждевременных родов с использованием лабораторного теста.
лабораторный тест разделен на две переменные: lab_order , который может принимать 0 (лабораторный тест не был заказан) или 1 (лабораторный тест был заказан) и lab_result , который является числовая переменная с результатом теста.
Я сделал это таким образом, чтобы отрегулировать тот факт, что lab_result можно заказать или не заказать.
Кроме того, лабораторный тест можно сдавать в разные недели беременности, поэтому я добавляю случайный эффект с помощью patient_index
для учета повторяющегося теста в разные недели для одного и того же пациента.
Модель, которую я запускаю, очень проста:
glmer(Term ~ lab_order + scale(lab_result) + (1|Patient_index),family=binomial)
но я получил следующую ошибку:
Error in pwrssUpdate(pp, resp, tol = tolPwrss, GQmat = GQmat, compDev = compDev, :
(maxstephalfit) PIRLS step-halvings failed to reduce deviance in pwrssUpdate
Я запустил его с и без масштабирования числовой переменной, потому что я прочитал, что эта ошибка может быть из-за проблемы масштабирования, но ошибка все равно не решена.
Есть идеи о том, что происходит?