ぎょーむ日誌 2004-02-26

2004 年 02 月 26 日 (木)

• 0750 起床． 朝飯． コーヒー． 0840 自宅発． 晴． 水たまりばかりの路上を歩いてるうちに， 今まで使いかたのわからなかったアカマツ部品データの使いかたを思いついた …… 確率分布の変化から， 成長速度 & 死亡速度を最尤推定する， ってのはどうだろう? 検討しつつ 0850 研究室着．
• で， ホントのところは昨日いろいろひねったぶんについて， とりあえず作図でもせんといかんところではあるんだが …… 「サイズ分布」(と仮によんでおこう) の相違から死亡速度と成長速度を憶測する， という計算が可能かどうかを調べてみる．
• いろいろ仮定が必要になるねえ ……
  • 現時点の age A の分布から A + 1 の分布が生成されるとする (ある種の「定常」のごとき状態)
  • A の分布は高次 logistic 分布に従う (ま， べつにパラメトリックな確率分布ぢゃなくてもいんだけど)
  • 死亡に関しては一次 logistic モデルで表現する (サイズが小さいと死にやすい，あるいはその逆)
  • 成長 (サイズの変化，射影) は何か確率分布に従う …… 例: ガンマ分布， 平均 = (定数) × (もとサイズ)， 分散 = (定数) × (平均)， これ以上に複雑なモデルの推定はかなり難しそう (死亡と区別がつかん)
  • 以上の確率分布を convolution したものが， 観測される A + 1 の確率分布である
  ということで， これまた複合確率分布の最尤推定問題に帰着されてしまった． どーせ数値的にしか解けんので， 成長の確率分布は上のようにやるのではなく， もっと conditional な状況を考えるべきかもしれん． 生き残ったものであれば， かならずサイズが増大する，というような …… あ， そうか． つまり「絶対」成長量が何か非負の確率分布にしたがう， というような状況ですな．
• てな考えをマトめてると時刻はすでに 1130 すぎ． うう
• ともかく， これまた計算がたいへんそうなので， とりあえず今までの中から多少は使えそうな部分をまとめる， という作業を先にやることにしてみる．
• ThinkPad X31 を修理に出すんで， 二ヵ月間ほど休眠させていた 240Z の戦力化作業 …… といっても rsync (あれこれデータの同期つまり差分バックアップ) と apt-get (システムまわりの最新化) すればいい． これだけなら 10 分もかからずに片づく． 自分の仕事環境を 10 分未満で別の計算機に完ペキに移せるようにするためには， 普段からのココロがけしかあるまい …… つまり 「いつでもネットワーク経由で差分バックアップが取れる」 ような整理整頓とシステム設定やってるかどーか， ということ．
• しかしながら 「ついでだから」 と Linux カーネルの作りなおし (2.4.20 → 2.4.25) とかやるから時間を食うんだよね …… えーと， IrDA の設定 はいぜんにどこかにメモしといたはず．
• 作業終了． X31 は止めて 240Z に移動． 外は風雪か． ということで北大構内走はあきらめて， ThinkPad の箱をかかえて北大生協へ． さーて， どうなることやら．
• 研究室にもどってくると， ドリルで壁に穴をあけたり鉄管をぎゃぃーんと切断する騒音が …… 有機溶剤のニオイもするぞ． なんだかよくわからん目的の修繕工事のようである． 学年度末における予算消化で忙しいんだから 「研究教育環境には配慮してるふり」 などなどの下らぬ「建前」なんぞはあっさり捨て去ってしまう， というのもまた北大事務の流儀ではあるんだが．
• 気にせず， アカマツ解析のつづきを， と思ったところで北大生協からお電話攻撃． ThinkPad まわりの部品いっさいがっさい持ってこいって? こちとら X31 を人質に出してしまったから， 何でも命令できると思ってるな (修理に時間かかって困るのは私なんで)． えーい． とりあえず， 小林さんに「仕事おくれてます」 いいわけメイルを書く．
• また風雪の中， 北大生協往復の旅． まあ， 歩いてる間もアタマの中では計算のやりかたとか考えられるし …… って， これは「しばらく放置しとく」はずの計算ではないか．
• とりあえず簡単にできる部分は簡単にやってしまって …… と単純な資源分配などの等分散正規分布なモデル選択をやる R プログラム書き． 作図する機能もつけてみるけど， まあなんともいんちきくさいかんぢが． とりあえずのメイル書き …… もう時間か．
• 今日はインドネシアから来てるユニさん (もうすぐ帰国) とカニを食べる会なのである． 1835 研究室発． けっこう雪がふったりやんだり． 1900 カニ料理屋着． カニばかりひたすら食べる …… というコース． 2100 同発． 2140 帰宅．
• [今日の食卓]
  • 朝 (0810): パンケイキ．
  • 昼 : 食ってない
  • 晩 (1910): すすきののカニ料理屋 えびかに合戦 でカニつくしコース． カニばかりだったけど， けっこうあきないもんだ． さすがわ専門店， というか．