2012-02-10 Rでトービット・モデルによる打ち切りデータの推定 の「2. optim()関数で最尤推定」で ifelse を使っているのは，ダメダメ。 掛け算の方がまだ速い。これで，元のプログラムより 25% 速くなった（元のプログラムが 0.104 sec. で，このプログラム…

配列の拡張とパフォーマンス 差を強調したいがために，基本を忘れている。 > num = 50000> ### 配列を拡張していく ###> gc(); gc() used (Mb) gc trigger (Mb) max used (Mb)Ncells 340948 18.3 597831 32.0 597831 32.0Vcells 602119 4.6 1162592 8.9 1162…

2012-01-28 Rでマルコフ連鎖モンテカルロ法を試す の後半「3. MCMC（M-Hアルゴリズム）で推定してみる」で， for (n in 1:10000) を 100 倍の for (n in 1:1000000) にしたら，10 分経っても終わらない。Stop させてみるとやっと予定の半分くらいしかできて…