Top

個別メモ

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

【第1964回】 偶然（2011年8月21日）

• 朝，祭りの準備で用具チェックをする予定だったので7:00に集合したのだが，天気が悪くて1週間延期になった。
• Warningが出る原因は，パラメータの定義域として正の場合しか考えていないのに，どれかが負になると返り値がNAになってしまう場合があることらしいとわかったのでコード改良の方針が立った。本来ならoptim()の引数としてパラメータの定義域に制約をつけられればいいと思うが，その方法がわからないので，値を最小化するパラメータを見つけたい関数の中で番兵をつけることにした。つまり，パラメータのどれかが負であったら計算をせず，確実に最適ではない値を返すようにした。これでWarningはでなくなったし，無駄な計算を省けたので，さらに高速化されたはず。optim()を確実に収束させるためには，control=list(maxiter=10000))のmaxiterを増やすといいのだが，ここまでもってくれば増やせるかもしれない。
• あとは計算させて2系統のグラフを書いて，もっともらしい解釈ができれば短報にはできるはずだから，英文も平行して書こう。面白い結果ならLancetに出してみたいところだが，面白くなくてもAm. J. Epidemiol.とかに出したい。しかし，最悪の場合でも，新しい方法を含むRのパッケージは作れそうだから，R Journalには載せられるだろう（書き方は変わってしまうが）。
• とりあえずデスクトップでプログラムを実行させ始めたので，続きは草取り後にしよう。
• ちなみに今日の見出しの「偶然」とは，今日が自分の誕生日なのだが，このメモの連番を西暦年と見なせば誕生年だったという，他愛もない話。偶然思いついたアルゴリズムが素晴らしい結果を出したとかだったら嬉しいのだが，そんなことはなくて，研究は地道にしか進まない。研究成果についてのセレンディピティ(serendipity)は属人的なもので，今のところ，ぼくには縁が無いようだ。
• 気温は低いのだが，蚊はしぶとく活動していて，長袖ウィンドブレーカを着て草取りをしている間に，顔を３ヶ所刺された。やはりTropical StrengthのInsect Repellentを顔に塗ってからやるべきだった。

▼前【1963】（別解（2011年8月20日） ）　▲次【1965】（打ち上げ（2011年8月22日） ）　●Top

△Read/Write COMMENTS

Notice to cite or link here | [TOP PAGE]