私は現在、オプティマイザーを使用して対数尤度を計算するモデルの matlab セットfminunc
を R を使用して変換していますoptim 'BFGS'
。
すべての matlab モデルの初期値、最尤値、および最終パラメーターの結果があります。R 変換されたモデルのほとんどは、matlab と同じ初期パラメーターを使用して、optim、同じ対数尤度、および同じ終了パラメーター値を使用して見つけることができます。ただし、初期値として matlab の end パラメーター値を入力することで修正できる局所的な最適値に行き詰まるモデルもあり、これらのモデルは matlab の最尤値を見つけます。
Rのより強力な最適化はありますか、それはmatlabのものと同等ですか、それともRがローカルオプティマで行き詰まる可能性が高いため、最初のパラメータ値がより重要になり、立ち往生する最大のログ可能性を達成するために重要になりますローカル最適?
res<-optim(par=x,fn=BOTH4classnochange,hessian=TRUE,method='BFGS',control=list(maxit=MaxIter,abstol=TolFun,reltol=TolX))