Сходится в многочлен

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

Я написал код Newtow Raphson R для оценки полиномиальной логистической модели, когда она сходится, и я получаю другое значимое значение из своего логарифма и встроенной функции R multinom. Это вывод из моего алгоритма

> t(new$coefficients)
        X0          X1        X2         X3         X4
1 1.055977  0.08972323 0.7070484 -0.4951407  0.1878713
2 1.845272 -0.26688254 0.6945098 -0.2749830 -2.0379031

а вот и вывод из встроенной функции

fit <- multinom(fm, data=df, Hess=TRUE)
> coef(fit)
        X0          X1        X2         X3        X4
1 1.055960  0.08972217 0.7070560 -0.4951362  0.187872
2 1.845253 -0.26688462 0.6945132 -0.2749767 -2.037911

Любой может помочь мне понять, почему это значимое после 4/5 десятичного числа. Пожалуйста, если вы хотите увидеть мой код, я поделюсь им с вами в частном порядке, потому что я не готов его опубликовать.

Спасибо

Источник

1 ответ

Решение

Расхождение между вашей функцией и выходами встроенной функции R вполне приемлемо. Этого можно ожидать, так как две функции этого типа, закодированные по-разному, будут нести разные незначительные ошибки, округления и т. Д.

Источник
Другие вопросы по тегам