多项式的收敛

Question

我希望你不要认为这个问题没有问题，因为我一直在努力解决这个问题，而且我不明白为什么它不会收敛到同样的价值。

我写了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十进制之后。如果您想查看我的代码，请私下与您分享，因为我还没有准备好发布它。

谢谢

Answer 1

您的功能和R内置功能输出之间的差异是完全可以接受的。可以预期的是，这种类型的两个函数（编码方式不同）会带来不同的小错误，舍入等等。