最尤推定プログラムが動かない原因の1つが

尤度関数を2つに分けて最尤推定値を計算するようにしたことで、パラメータベクトルの分割がうまくいっていなかったことが判明。

次週の目標： プログラムを修正して最尤推定値っぽい値を計算する。