ぎょーむ日誌 2005-08-07

2005 年 08 月 07 日 (日)

• 0800 起床． 朝飯． コーヒー． また怠業．
• お． 今日も曇ってきたな． 1110 自宅発北大構内走． 夏は代謝速度が高くなるのか …… 走ってる途中に急に空腹になることがある． いわゆるシャリばて． 1205 帰宅． 体重 74.4kg． 昼飯．
• まーた冷蔵庫の製氷部分に厚さ数センチの結氷． 中身とりだして電源おとしてドア開放して下に洗面器おいてしばらく放置． 一時間以上かかる． 小人閑居して不善を為す， ということで先日の biometry への質問投稿への返信． これまたちょうど一時間ぐらいかかった．
• 結氷はがしが終わったのに小人閑居もーどがつづく． 来月の 統計関連学会連合大会 の見物 (ホントに見物だけになってしまうわけだが) に行こうかと思ってんだけど …… 参加するとして宿泊どうしようかナヤむ． 大会会場がですね， 広島市内とはいえ 宇品のプリンスホテル ってへんぴなとこなんだよね． ぢゃあ， そこに泊まればいいぢゃん， とゆー案もあるんだけど …… 得失を勘案すると難しい． よさげなところは
  1. 元宇品公園横にあり， うんどう不足にならずにすみそう --- ジョギングコースいろいろありそう
  2. 無線 LAN は使える (hotspot)
  3. (ホテルが大会会場なんで) 移動しなくてよい
  しかしウラをかえせばヤなところになるわけで
  1. びぢねすほてるではなく観光地なんで値段高い (さらに飯もここで調達するとなると……)
  2. 有線 LAN ない
  3. 逃げ場 (?) ない
  うーむ ……
• ぢゃ， 会場から離れててもいーから， 良さげで LAN 使えるところ， と探してみるんだが …… 広島市って意外とないんだよねぇ．
• あれこれ検索． ウェルシティ広島 (広島厚生年金会館) ってどうかしらん． 一泊 8300 円ってべつに安くないぢゃんと思われるかもしれないけど， ぢつは「厚生年金被保険者」なる連中は 6900 円ですむんだよね．
• さーて， 問題のネット環境だけど …… 楽天トラベルの「モジュラージャック、インターネット接続 (LAN形式)」 なる記述はいまいち信用できん． 自サイト には何も記述ナシ． 不安． しょうがないので電話 (082-243-8881) で確認してしまった． なんとゆーか， 電話ぎらいのくせにそこまでして， とは思いつつ …… あ， なるほど， ホントに LAN の DHCP 接続でしたか (ケイブル貸し出しアリ)． まあ， 公園とかも近くにないわけではないし． 平和記念公園 とかだけど．
• 1550 自宅発． 曇． 1605 研究室着．
• 来月のすけぢゅーる (しかも「見物」だけ，という) を今ごろ確定してまっていいのだろうか， と疑いつつ …… まあ， 急にダメになったり 「広島なんか行ってる場合では」 と追いつめられたらもうあきらめましょう， ということで準備をすすめる．
• まずは 参加申し込み． 大会参加費 4500 円． 懇親会 5500 円， か． どうしよかな． ぢつはこの統計ゆにおんってほとんど知人いないんだよね． まあ， 根性で「参加」にするか．
• で， 支払いはクレジットカードでできる， と …… 生態学会関連がくれぢっとかーど決済になるのはいつのことやら． 何やら三井住友 VISA の V-payment なるサーヴィス使っているな． 「学会センター」って， まさかあの ……?

• 北大生協 で航空券手配． これはいつものごとく．
• 宿泊予約 …… ここに 「モジュラージャック付き LAN 設備」 とか書いてあったんだけど 「もぢゅらー」なる表現がちょっと不安だったんだよね． 厚生年金な施設のオンライン予約アカウントでログイン． 予約． どうか一ヶ月後にすけぢゅーる破綻にやられていませんように， と．
• 雑用というか， なんというか， ともかく一件落着．
• さて， (昨日とりあえずの一区切りとなった) アカマツに予想外にてこずってしまって， すけぢゅーりんぐの再考をせまられているところである．
• 今年前半の目的としては， 前年の釧路大会雑用 (cf. 一年前のぎょーむ日誌) でこうむった被害を圧縮し， これまでの仕事の「不良債権化」 (計算やったけど論文として出て行かない) を阻止する， という基本方針で．
• さて， この基本方針にもとづいて
  • 苫小牧 leafminer 寄生モデリング分担ぶんの再計算とまとめ
  • 苫小牧樹木直径成長モデルの解析まとめ (こちらも再計算必要?)
  を優先し， 「新しいこと」 にはできるだけ手をださないことにしたかったわけだが ……
• しかし今までのんびりしすぎてしまったので， どうも今月あたりはそうも言ってられなくて， 「新しいこと」 を先に始末しちまったほうがいいかも， というその場しのぎな状況に追い込まれつつある． 時間順で考えると
  1. 井田君の計算したうけ: 混合幾何分布の最尤推定など (3-4 日間ぐらい?)
  2. 8/18-8/23 (?) 小林さんがこちらに来られるので: 母子里葉群の三次元分布の MCMC 推定計算の基本わくぐみの構築 (10-20 日間ぐらい?)
  3. ツガザクラ下請けモデリング: ECOAS 「生物相互作用に着目した高山・亜高山生態系の 脆弱性評価システムの構築に関する研究」 なアレ (どれだけの労働時間を費せばいいのか， まったくわからない)
  あたりか． 1. は 進化学会大会 (仙台) がらみなんで， さっさとどうにかしたい． 2. は「モデル屋に仕事させる」方法論にますます熟達しつつある 小林さん監視下での作業なので逃れられない． いっぽうで， 3. は河田さんがらみなので逆に進化学会大会が終るまで (あるいはその後の あとしまつ・雑用ばて・もえつき期間 も含めてよい?) 放置しといてさしつかえなさそう． しかしながら， それ以上の サボタージュ はかえって危険な気がするのである (計算のみとおしがまったくたってないんで)．
• これらの 「新しいコト着手」 がアブないと見なされる最大の理由は， 論文書きとかやってるよりよほど 面白くてそれゆえに中毒性というか依存性が高くてぢゃんきー …… ということにつきる． ひどく危険だ．
• しばらく放置しても文句いわれそうにないのは …… 先日のアリコロニー戦争とかかなぁ． PipeTree の gcc-3.x ばぐとか早く修繕したいんだけどね． そして MCMC な推定計算の実用化研究， とかも． 今年度の修論下請けなんかはすごく不安なんだが．
• まずは混合幾何分布， と思ったんだけど空腹になってきたんで撤退することに． 1935 研究室発． 買いもの． 2005 帰宅． 晩飯．
• で， 混合幾何分布の推定問題について検討してみる．
• そうさなー， 世の中には beta geometric distribution とかヘンなこと考えてるやつとかってたくさんいそう …… と検索してみたら論文になってるのは少なかった． やはり使いにくいからだろう．
• で， 使いやすさで考えると …… gaussian logistic geometric distribution とか? さすがにこんなめちゃくちゃな名前を使うやつはいなかった． ぢゃ， mixed geometric distribution ならば? これは多少あるな．
• などとばかげたフレイズどもをアタマの中でひねくってるうちに …… いやー， けっこういい計算方法おもいついてしまいましたよ! 一言でいえば， 幾何分布 系のパラメーター推定も glm() 系の関数でやっちまえばいいんだ， ということ． 推定だけ ならば．
• R には幾何分布のパラメーター最尤推定用の関数はない (ま， そのうち誰かが CRAN に登録するだろうけど)． ましてや混合幾何分布をや． しかしながら， ここで発想を少しきりかえると glm(..., family = binomial) だの glmmML(..., family = binomial) だので幾何分布のパラメーター最尤推定計算がきちーんとできてしまう， とゆー理解に到達するわけで ……
• もうちょい詳しい説明は明日のぎょーむ日誌にでも． ともあれ， 混合幾何分布推定のため対数尤度を最大化する optim() よびだす関数をこしらえないといかんのか面倒面倒だなと思ってたんだけど， なんのこたァーない glm() とか glmmML() でケリがついちゃうとはね …… 「戦わずして勝つ」 「兵は詭道なり」 (孫子)， は大げさか．
• うん． これはきっと昨晩， 荒木飛呂彦の新刊 「スティール・ボール・ラン 5」 読んだ御利益にちがいない．
• [今日の運動]
  • 北大構内走 55 分間．
• [今日の食卓]
  • 朝 (0830): パン． ネギ・ニンニク茎・豆腐の炒めもの．
  • 昼 (1240): 蕎麦．
  • 晩 (2110): パン． タマネギ・ニンジン・ピーマン・エリンギ・豚ハツのカレー．