SIR のための尤度計算

多数の条件（パラメータ値）それぞれで尤度を計算して出力するために by 変数でループをまわす必要がある。これが大変だ。