stata - Stata は最尤推定で私の開始値を無視しているようです

Question

最尤モデルを推定しようとしていますが、Stata で収束の問題が発生しています。実際のモデルは非常に複雑ですが、適切な開始値が与えられた場合、R で問題なく収束します。しかし、私が提供する開始値をStataに受け入れさせることができないようです。

ポアソン分布の平均を推定する簡単な例を以下に示します。これは私が推定しようとしている実際のモデルではありませんが、私の問題を示しています。トレース変数を設定すると、Stata が尤度面を検索するときにパラメーターを確認できます。

init を使用して開始値を 0.5 に設定しましたが、最初の反復では、Stata が係数 4 を試行していることを示しています。

どうしてこれなの？推定手順で開始値を使用するように強制するにはどうすればよいですか?

ありがとう！

generate y = rpoisson(4)

capture program drop mypoisson
program define mypoisson

args lnf mu

quietly replace `lnf' = $ML_y1*ln(`mu') - `mu' - lnfactorial($ML_y1) 
end

ml model lf mypoisson (mean:y=)
ml init 0.5, copy
ml maximize, iterations(2) trace

出力：

   Iteration 0:
   Parameter vector:
        mean:
       _cons
   r1         4