ぎょーむ日誌 2005-10-21

2005 年 10 月 21 日 (金)

• 0620 起床． 朝飯． コーヒー． 0730 自宅発． 晴． 0745 研究室着．
• 論文セミナー予習．
• 1000 より 論文セミナー， 今日の担当は石橋君で BCI の小さい樹木の太さ成長 vs 節間長関係の解析． このあたりのあろめとらーな解析はいつ見てもアタマ痛い (…… がその改良は未だ「これからの問題」である)． なんでも等分散正規分布とか， log-log 回帰の平均値だしそこねとか， 「何でもいいからゆーい差さえあれば」主義とか， 推定値を使いまわした二重推定→検定だとか …… 要約すれば， 一見もっともらしいけどじつは無意味な数値こねくりまわし (データ解析に関して著者 King と甲山さんの間にいくらか発想の類似性がある， とわかったのは本日の発見ではあった)． しかし解説者はそもそもこの King 統計学がわかってないんで， King 統計学が正しいとおいた上での質問・確認に終始せざるをえなかった． 1140 ごろ終了．
• 昼飯．
• なんかアタマかきまわされた状態． これでは「データ解析しーづん」をのりきれない． えーと ……
  1. 苫小牧 leafminer 寄生問題: 現モデルの改善を考える …… そうそう， 今朝学校に来る途中に 「そーだ， Gaussian Random Field (GRF) を MCMC ではなく (spatial) pseudo likelihood で取り扱えばいいんだ」 という思いつきを得た． これはマトモなのかどうか， ちょっと検討してみる価値がある．
  2. 苫小牧樹木直径成長問題: モデルの変化に対応できるようにする …… 昨年 11 月にころに試験してた新しいモデル選択法を 実用化してみるべきか?
  3. 屋久島照葉樹シュート問題: シュート伸びる・伸びぬ問題はとりあえず牛原さんにおまかせ． いっぽうで， 欠測を確率論的にうめる葉寿命推定は未着手のままだな．
  4. ほかにもいろいろありそう: ツガザクラ問題とか放置しとくのは危険なんだが (いきなり「結果だせ」とか言われそうなんで) ……
• 1. の GRF まわり， 私の苦手な (とゆーか全然いままでまともにあつかったことのない) 多次元正規分布の計算が必要になる． ちょっと紙の上に尤度方程式とか書いてみたんだが …… これってさー， MCMC 法でやったほうがいっそラク， というたぐいの問題であるようなきがしてきた． 分散共分散行列の最尤推定なんかやってられるか …… とゆー気がするんだけど， 実際にはできるんだろうな． 各地点における「隠れ」変数のあつかいが難しいのである．
• と思ったら苫小牧の平尾君からメイルいただく． なんだかデータファイルの改訂の過程でややこしー事態が生じていたようで．
• 甲山さんの科研費申請ごたごた． 今年からネット上であるていどは作業ができるようになったのは 少しエラいと思う． システムは現時点では未熟であっても． しかし， このぶらうざー依存性ってのがちょっときつい． 甲山さんは昨日 Safari (MacOS X 上の web browser のひとつ) で作業してたんだけど， 結果的には作業途中の状態を「一時保存」できなかった． それでは， と MacOS X 用の MSIE5.2 で試してみると …… なぜかしらこれは動作しやがるんだよね．
• Gmail 音信不通 …… いや sinet あたりがおかしくなっていて， 国外につながらなくなってるみたいだな．
• 赤坂君の空間統計学こんさる． これってぜったいふつーの負の二項分布とかで説明できる現象ではない， とゆー気がするんだが． 自分でやっててど忘れのかなたに放置されてた library(spdep) まわりについて (ぎょーむ日誌今年の 三月下旬)， いろいろと思いだすことできた． 村上行列とかいろいろ工夫してたな ……
• 距離カーネルについてはいつもいつもわずらわしい状況になりがちなので， 七月ごろ に研究してた高速計算方法であるていど汎用性のあるライブラリとか作ったほうがよいかも．
• 竹中さん に生態学会さーばー雑用おしつけることになってしまった． ぢつは真にめんどうな部分は， だいたいにおいて竹中さんが処理されているのである …… が生態学会はわれわれをいつまで奴隷な技術者として使いつづけるつもりなのだろうか?
• 科研費まわりのめんどうは 粕谷さん に押しつけているし …… いやはや， どちらもおそれおおいことです． 粕谷さんから送られてきた承諾書のわーどファイル， OpenOffice.org 2.0beta でわりとうまく解毒できた．
• 1810 研究室発． 1825 帰宅． 晩飯．
• NHK らぢおのニュースによると， 本日の札幌の空を埋めつくすよーに飛び回っていたナゾの昆虫は ケヤキフシアブラムシ， 1995 年以来の 10 年ぶりに大発生だそーで (道新とかの報道もある)． そういや， 苫小牧 leafminer もしーづん終わりに「豊作」だったとか? どちらも気温との関連が指摘されてる． 昆虫個体数の時間変化とか， いつまでたってもつくづくわからん．
• 2240 自宅発北大構内走． 「一般化超幾何分布を事前分布とする推定……?」 といった正気さのいくぶん欠如したアイデアが夜道を徘徊している． かかる極論でアタマを刺激すると， 脳内りすとらがうまくいくことがある． 空間構造を介した相関の問題， われわれはかなり奇妙なとらえかたをしていた， というのがわかってきた． 惰性モデリングはまさに危険きわまりない． 2320 帰宅． 体重 74.4kg．
• [今日の運動]
  • 北大構内走 40 分間．
• [今日の食卓]
  • 朝 (0640): 米麦 0.7 合． ネギ・タマネギ・ニンジン・ブナシメジ・ハンペン・ サケの味噌汁． キャベツ酢のもの．
  • 昼 (1220): 研究室お茶部屋． 米麦 0.5 合． ネギ・タマネギ・ニンジン・ブナシメジ・ハンペン・ サケの味噌汁．
  • 晩 (1910): うどん． ネギ・タマネギ・ニンジン・ブナシメジ・ハンペン・ サケの味噌汁． キャベツ酢のもの．