対数尤度を計算できるようにしました. 次は, パラメータをサンプルしながら繰り返し計算させます