EPSON Endeavor のデスクトップを中古で購入しました。自宅用です。有線で LAN につなぎます。 安定していてよいですね。 数年前の i7 ですけれどこれで十分だ。 NONMEM, Phoenix, Monolix をインストールした。MATLAB も入れようかな。

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

これも SIR なのでやってみた。 単純なんだ。でも、実効性のある結果を出すためには大量の計算が必要になる。今のシステムでは取り扱いきれないんだな。

昨日の今日でなんですが、断念しました。 CPU パワーの問題ではなく、sample 発生のために必要となる「解析用データセット」が巨大になってしまう。例えば、 回シミュレーションしようとするならば、元のデータセットをその回数分積み重ねる必要がある。実際…

MCMC がまだ使えないので、SIR を使ってみた。SIR は初めて使ったのですが、あっさり計算できました。 個体間変動パラメータの事前分布はシンプルだから計算しやすいです。実に簡単にプログラムできた。元のサンプル数を 4000 個、resample を (≧20 がいいと…