ぎょーむ日誌 2004-01-03

2004 年 01 月 03 日 (土)

• 0740 起床． 朝飯． コーヒー．
• う． またぎょーむ日誌生成まわりの Perl スクリプト手直しに 時間を使ってしまう．
• いつまでも怠業ばかりしてるわけにもいかんので 1100 自宅発． 晴． 1115 研究室着． 誰もいない． 飼い主たちが帰省してるのでお茶部屋ずまいになってる動植物の 水やり・エサやりとか． しかし仕事すすまんので， とりあえず昼飯．
• さーて， 仕事， と思ったらまた KuboLog スクリプトにとらっぷされるわけで． きりがないな．
• シウリザクラ推定問題， 昨年のやりのこし (森君問題のほう) に取り組もうとする． 簡単にできる問題は片づいてるわけだが， あとはなかなかややこしいな．
• かとーオフィス内， 奥手稲小屋なみに寒い …… とゆーのは大げさにしても， 室温 12°C はなかなかキビしい． とりあえずガスストーヴつける．
• 甲山さんがうろうろしてる．
• 森君のシウリザクラ推定問題に半日費して
  1. 「父親としての成功する要因」 の最尤推定 (GLM ではない) & 可能な全モデルの列挙 & その中から AIC によるモデル選択 …… 今回も R の optim() 版を作るよゆーがなかったんで 自作 Perl モジュールの組合せで解決
  2. その最良推定結果と連動して 「母親としての結実量を高める要因」 の一般化混合線形モデル (GLMM) による推定 …… これは大澤君のシウリザクラ問題用に開発した R の glmm() (in CRAN) & stepAIC.glmm() (自作) のコンビネイションわざで解決
• ありていに言って 1. で得られた推定値をもって 2. の推定計算やるのは少しイヤなんだが …… まあ， 修論だったらいいのかな．
• また 1. ではクレーンサイトと緑トンネルサイト を「プール」しちまっている， というか両調査地のパラメーターが同一と仮定して， 推定計算をやらせるようにした． なにしろクレーンサイトの調査個体分布ってのは ひどく偏っていてですね …… いまさらながら， 生態学における個体間「距離」を考慮した統計モデリングは 難しいもんだと思う．
• ともあれ， そのように灰色な作法を用いて得られた結果は ふつーなしろものだ． 父親としての成功， 幹直径にほぼ比例， 距離が近ければよい， 同一遺伝子型はダメ． 母親の結実， 同一遺伝子型樹木に包囲されるとがた落ち …… だけか． 例によって混合モデルのほうがあてはまりが断トツに良く， しかもそれは樹木単位ではなく花序単位でばらつく， と仮定したほうがよい． このデータセットとモデルの場合， 樹木ごとのランダム変量導入はほぼ意味ナシ， ということのよーで．
• ともあれ計算に関しては一件落着． しかし， 報告を書かんといかんのだろーな． 2220 研究室発． 北 12 生協で買い物． 生鮮食品の入荷が少なくておもしろい． 2240 帰宅． 体重 72.4kg． 晩飯．
• 昨日も今日もプログラミングばかりやってたよーな． 2500 消灯．
• [今日の素読]
  • なかなか再開できんなぁ
• [今日の運動]
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0800): 米麦 0.7 合． 納豆．
  • 昼 (1200): さっそく営業してた北 12 駅よこのパン屋 「北欧」で調達したあれこれ． ひじきパン 120 円． まずくはないしミネラル補給にもなるんだろーけど， やはり無理矢理感がつよい． カレーパン 90 円． 揚げてないやつ． こちらはふつー．
  • 晩 (2320): 米麦 0.7 合． ネギ・豆腐・カキ・コンブの味噌鍋．