個別メモ

Latest update on 2012年3月5日 (月) at 10:54:46.

【第1957回】プログラミングする日曜（2011年8月14日）

久々に8:00まで惰眠を貪った。人口分析のテキストのドラフトを書き始めたり，死亡モデルのCのコードをRに移植し始めたり。当然，一日では終わらないが，死亡モデルの方は並行して英文も書き始めた。
Cに比べるとRのコーディングは非常に効率が良く，Cで10行必要なところがRだと1行で済んでしまう。実行速度が遅いのは速いコンピュータを使えば解決する問題だし，もうCには戻れないなあ。たぶん計算をベクトル化すれば，もっと速くなるはずだし。
Cでは乱数発生とか関数の最小値を求めるとかいった操作も１つ１つ関数定義する必要があったが，Rでそれをやっていると計算時間がかかって堪らないので，乱数発生はRNGkind("Mersenne-Twister")で，関数の最小値を求める部分はoptim()で置き換えることにした。たぶん，それだけで実用的な時間で結果が得られると思う。後はデータとのインターフェースで，Cでは予め加工してから読み込んで，結果出力は数値をファイルに書き出して，グラフ作成はgnuplotを使っていたが，Rでは全部内部でできるので，相当に効率化できるはず。

▼前【1956】（土曜（2011年8月13日））　▲次【1958】（月曜は普通に出勤（2011年8月15日））　●Top

△Read/Write COMMENTS