Я пытаюсь решить алгоритм Ньютона в K измерениях, не инвертируя матрицу гессиана, потому что она возвращает ошибку, когда K> 4 (функция convert и ginv возвращают матрицу 0 (гессиан не является обратимой)) Чтобы избежать этого, я пытаюсь решить линейную систему алгоритма Ньютона
GL=(hesse%*%b)-dL
#solving the linear system
b=lsolve.bicgstab(hesse,GL,maxiter=1000)$x
Я пытался использовать различные функции из библиотеки Rlinsolve, но всегда возвращалось сообщение об ошибке
сходимость не достигнута в пределах максимума.
Я использую алгоритм Ньютона для регрессии логистики c Я уверен, что могу пропустить эту ошибку, потому что при использовании функции glm я могу получить результат. когда я в K <4, я получаю тот же результат, что и glm (что означает, что остальная часть кода верна). Есть ли другой способ решения линейной системы в r? </p>