我首先要声明的是，我无法确定这是 BFGS 问题还是 R 包问题。

我一直在使用 R 包“mlogit”进行一些混合 logit 回归。我总共有 32 个随机参数和一个要估计的固定参数。随机参数之间允许相关。使用 BFGS 方法，代码在 13 次迭代后停止。从卡方值来看，远未收敛。然而，使用 BHHH 方法，它通常会收敛（尽管非常缓慢）。这是包装问题吗？还是在我的上下文中它比 BFGS 更适合 BHHH 方法？

我怀疑它是否与步长有关，但这只是我的怀疑。我使用的代码和数据在这里。