ぎょーむ日誌 2004-04-(01-10)

2004 年 04 月 01 日 (木)

• 0820 起床． 生活周期が修復しない． 朝飯． コーヒー． 0920 自宅発． 晴かつ雪． 0930 研究室着．
• 昨日しめきりだった公募シンポジウム・自由集会もうしこみを 事務局作業用 wiki はりつけ作業のこり片づける …… なんつーか， 生態学会ってすごいイキおいで土建屋じみてきてるんぢゃないのか? そういうことをやらざるを得ないのか， それとも， そもそもそういうことをやりたかったのか．
• いやはや …… 苫小牧直径成長モデリング， ちょいヘマをした． 不幸中のさいわいというやつで， 推定結果を左右するものではなく， 単にまぬけなむだ計算量を増やしていた， ですんだんだけど． 2003 年は光量子密度 (PPFD) の観測値が 6 月ぐらいまでしかないんだった． で， このへんを回避するために気象値は 気温 × 降水量だけで計算してたんだよね．
• さて， ここで「前年の光合成量だけが重要」 という極端な仮定を置いたとしよう． いやいや， 苫小牧データみるかぎり， 当年の気象状態も重要かもという気にさせられるんだが …… それはまた別の種類の気象フィルターが必要らしい， という気がする． ということで， 限定された時間でよい推定のため前年だけ， と限定すると 2003 年の PPFD 観測データが足りないのは 問題なくなる (1997-2002 年の期間はだいたいそろってるため)．
• というような状況の計算やらせてみる． もはやこのあたり R のプログラムたちを連動させるしくみが ほぼ自動化されてしまってるんで， 計算やりなおしなんかが簡単にできるようになってる．
• 「ついにエコカップが日本フットサル発祥の地、北海道にやってきました」 ということで エコカップ2004 サイトが永光さんによって準備 & 公開された． 大会ペイジからリンクする． ときどきサーヴァーがオチてるようで …… 全国のサッカー好き生態学者たちがいっせいにアクセスしてるためか?
• 2003 年の PPFD も送っていただく …… これで解析に必要なデータは， だいたいでそろったはずだ．
• 計算走らせつつ北大構内走に出る． もどって昼飯．
• なぜか PPFD をくみこむと計算の収束がひどく悪くなる． うーむ……?

[光量子密度]

って， そもそもこんなにちゃかちゃか揺れるものなんだろーか． 天気悪いと日累積がゼロ近くにまで低下する? この落ちつきの無さが， 推定計算に悪影響およぼしてるよーな気がしてならん …… とゆーのは責任転嫁だろうか．

• えーい， また釧路大会関係の …… しかも， 甲山さんに撹乱されてしまって わけわからん質問をメイリングリストに投稿してしまった． すぐにおわび書き． 皆さん， 参加費のしはらいは なるべく和文誌とじこみの払込票を使ってください．
• あらら …… 生態学会事務局のペイジから釧路大会サイトへリンクすると， frame の中で開いてしまってじつにマヌケだ． というか， 実際のところすごく不便で不都合だ． さっそくペイジを管理してるらしい 広島大にそんなことやめてくれメイルを送る． いやはや． 意味もなく frame なんか使うなよな． ちなみに釧路大会サイトではまったく使ってない．
• その他ネット雑用を終ってみれば時刻はすでに 1630． ウラで計算はさせてたけどね． 収束がよろしくない， というのが確認できただけかもしれんけど．
• お茶部屋 iMac (MacOS X 10.3.3) がふちょー …… かとーさんに診てもらうと， command + option + P + R 起動の「PRAM クリア」わざ (まだあったのか!) で解決された． なんでも OS X になってもアプリケイションソフトウェアが PRAM 値 (これはいわば firmware = BIOS の設定値ごときものか?) を書き換えやがることがある， ということらしい． なるほど …… って， そういう腐れたソフトウェアがあるわけね． command + S 起動で single user mode， というわざも教えていただいた．
• 本日から M2 になった平林さんの R こんさるてぃんぐ． お， apply() わざを修得しつつあるようで．
• 直径成長モデリング， 苦闘の再開． PPFD もくみこんでなんとかまっとうそうな答えをだしたい …… が， 数時間ついやして， こいつはどうも一般化線形混合モデルでは無理なんでは， という気がしてきた． 少なくとも R の glmmML() ではうまく解けんよーな． かのおそるべき start.sigma 依存性を確認してしまったからだ． データが足りんちうこと? ううう ……
• しょうがないんで glm() + stepAIC() にもどす． これは計算むちゃくちゃ速い． しかも確率 1 で着実に収束する． 推定結果の内容も， まぁ悪くないのかなあ …… 尤度がすごくヘボいってことを除けば． すなわち， やっぱ overdispersion ありありまくりの混合モデル使わせやがれコノやろ状況ということですよ， これは．
• まあ， こういうこともあろーかと問題をあらかじめ離散化しちまってるから， 一般化線形「負の二項分布」モデル glm.nb() (in MASS library) 利用という悪辣なる策も残されてはいますな …… 待てよ， これって stepAIC() はウケつけないんだっけ? 実験してみる． ダメぢゃん …… Ripley 御大は何を考えておられるのか? stepAIC.glmmML() を自作したのに続いて stepAIC.glm.nb() も作れって試練なわけ?
• といいますか， 便利なのか不便なのかよーわからん glmmML() やめて自分で計算関数かけば混合 Poisson 分布でもイケると思うんだけど．
• しかしながら， とりあえずは， 新アイデアたる気象フィルターの性質を把握することが先決だ． いかな「粗い」結果とやらが発注内容とはいえ頑健ぢゃないと 精神衛生によろしくないので， stepAIC(glm()) を使って時間をかけてひととーり調べるよう命じておいて 1940 研究室発． 2000 帰宅． 体重 73.8kg． 晩飯． なんか雑用 & 計算ばてばて．
• [今日の運動]
  • 北大構内走 1310-1345． ストレッチング．
• [今日の食卓]
  • 朝 (0840): 米麦 0.8 合． 卵焼き．
  • 昼 (1410): 弁当． 研究室お茶部屋． 米麦 0.7 合． ハクサイ・ニンジン・タマネギのカレー．
  • 晩 (2200): 米麦 1.2 合． ハクサイ・ニンジン・タマネギのカレー．

2004 年 04 月 02 日 (金)

• 0750 起床． 朝飯． コーヒー． 0900 自宅発． 曇． 0910 研究室着．
• さーて， 324 とーりの気象フィルターについて調べさせていた昨晩の計算結果だが …… まあ， こんなもんか， という推定結果． ついでに記録しとくと， 324 とーりの気象値セットを生成するのに 50 分ぐらい， stepAIC(glm()) な計算は 1 分ほど使っている．
• でてきた推定結果を視覚化する作図を工夫してみる． ついでに個体ごとの挙動を重ねあわせてみたところ， つくづく混合モデル的な状況だというのがよくわかった (今回は random effect ヌキの推定計算なんだが)． また苫小牧にご注進メイル．

[太さと太さ成長の関係]

…… が年によって異なり， それは前年 & 当年の気象に影響されてるんだろう， という仮定のもとで推定されたモデル (とその「もと」になった観測データ)． たしかに， 樹木のサイズによって， 年変動する気象への応答のしかたが異なるように見える (そしてそのように推定された)．

• 昨日依頼しておいた 生態学会事務局ペイジから釧路大会サイトへのリンク， すばやく変更していただいた． 助かります．
• 午前の残りは， Matt の熱帯林のろわれモデリングにまきこまれる． 限りなくややこしくなりそうな問題に歯止めをかけるべく， 問題の原点にたちかえって検討しなおす． すると ``global'' と称してる分布の計算のやりかたに， もうひとつの分岐が存在することに気づいた． こちらのほうがやや人工的ではあるけれど， 同時に現存種の性質にあまり左右されない， という特性があることに気づいた．
• と朝からクタびれてしまって， 北大構内走にでも出ようかとしたんだけど …… 真冬のごとき豪雪． ということで， 昼飯．
• 午後は甲山さんのアヤしげモデルにとらっぷされる …… いや， かなり単純化されているのでいいかんぢにはなってきている． しかし 5 変数あるんで， 見とおしはよくない． 数値的に結果を示すほかないんではないかな． で， それをやりたくないんだったら， わりと面白い問題だとは思うんだけど， それはそこまででしょう． 私としては， 授業で院生に 「簡単な森林動態モデル」 みたいなことを教える機会なければ使うことなさそうな．
• でも， 教育につかえる数理モデルは， それはそれで貴重ではある． 少なくとも， ものを考える出発点にはなりうるからね．
• ふーむ …… 生物学的にはあまり意味のない変数であっても， 曲線のあてはめには貢献しうることがあるなぁ， と気づいた．
• で， 日没ごろまでまたまた熱帯林問題にまきこまれたり …… 無作為化した量をうまくつかえる局面が見つかったので， 少しばかり呪われ度が減少した． このあたり Matt の計算は素早い． ばぐをかなり出すけど．
• 熱帯林ばてしてるところで苫小牧からの追加発注来信 …… イタヤカエデ の次は アサダ ですか ……
• データ変換やら何やらで 2 時間ほどしてアサダの推定結果得られた …… おお， イタヤカエデとまるっきり違うぢゃん．
• アタマが働かなくなってきたので， 計算機力まかせな演算を命じておいて 2040 研究室発． 2100 帰宅． 晩飯． ばて．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0810): 米麦 0.8 合． ハクサイ・ニンジン・タマネギ・豆腐のカレー．
  • 昼 (1250): 弁当． 研究室お茶部屋． 米麦 0.8 合． ハクサイ・ニンジン・タマネギ・豆腐のカレー．
  • 晩 (2140): 米麦 1.0 合． ハクサイ・ニンジン・タマネギ・豆腐のカレー．

2004 年 04 月 03 日 (土)

• 0900 起床． 朝飯． コーヒー． 0950 自宅発． 曇． 1000 研究室着．
• また苫小牧直径成長モデリングにとりくむ． この「粗い」推定方法もそろそろ限界が見えてきたなぁ …… というのは確かなんだけど， ひとつ見落としてたポイントに気づいた． ただしこれは後になってから． 気象値なるものは (フィルターされた気温) × (フィルターされた降水量) × (フィルターされた PPFD) なる構成になってるわけだ． この 3 要因のうちひとつかふたつ落したい場合は， フィルターのパラメーターを， すごく極端にすればいいぢゃないか． ふーむ
• 1300 から 1700 まで釧路大会雑用会議． 公募シンポジウム 15 (!)， 自由集会 23 (!!) すべて採択されましたよ． おそるべし． というか， 生態学会大会なる集まりというのが， いろいろな方面において あれこれ破綻しつつあるような予感．
• で， 会議の席上， 私がこなさねばならぬネット雑用が， やはりけっこう多量にあると確認されてしまった …… どうなるんだろうね． あ， ポスター賞とかまたやるらしい． それだけならいいんだけど， これまたネット雑用増大に貢献してくれるんだよねえ．
• なんというか， ひどくバテ感の高い 4 時間でした． で， 終ったらそちら関係を整理すべきなんだが …… そういうぎょーむなんぞは放棄して， イタヤカエデとアサダのアヤしげなる世界に没頭 → ますます疲労度を高める． うう
• ふらふらと 2140 研究室発． 2200 帰宅． 釧路ばて & 苫小牧ばてで体重減か …… といってもまだ 73.2kg もあるわけだが． 晩飯．
• [今日の運動]
  • うんどう休養日になってしまった
• [今日の食卓]
  • 朝 (0920): パンケイキ．
  • 昼 (1230): 「北欧」のベイグルサンドイッチ．
  • 晩 (2310): モヤシ・ニラ・マイタケの あんかけやきそば．

2004 年 04 月 04 日 (日)

• 0900 起床． 朝飯． コーヒー．
• でれでれとする週末はよいものですな …… とでれでれしてる場合なのだろうか．
• 1210 自宅発北大構内走． 1300 帰宅． 体重 73.2kg． 減らないなぁ． 昼飯．
• 1410 自宅発． 晴． 1420 研究室着．
• 昨日から， 当年の気象要因として累積気温ってのを組み込もうとしている． 累積気温なるものはあれこれと趣旨不明瞭なところがあるんで， あまり使いたくないんだが …… とりあえず氷点下の気温を累積しても意味なさそうなんで， ゼロ℃以上を足していくことに．
• で， これを使うとアサダの成長モデルは多少マシになるような気がする． しかしながらイタヤカエデのほうはずたずたになる． いろいろと作図してみて影響のキキかたを確認する． ならば次の一手は …… パラメーター数をこれ以上は増やしたくないんだけど， その条件つき累積気温なるものを展葉前と展葉後に分割して， 線形予測子を構成してみる． よーし． これでイタヤカエデのほうはだいたいもとどーりになった．
• アサダのほうは …… うーむ， またまた逆気象フィルタリングになってるな． しかし， これによって現象はよく再現できてはいる． どうしたもんかのう．
• とりあえず昨日おもいついた 「極端なパラメーター値で気象要因を落す」 技法を適用して， 気象値テイブルをつくりなおしてみる． これでまた 2-3 時間かかりそうだな． 1608 計算開始．
• あ， こういうふうに極端値をいれて要因を落すのはよいとして， これは「使ってるパラメーター数」に反映されないぞ． えーい． このあたり修正するのに 1 時間以上かかってしまった． やはりこういう姑息なごまかしはダメなのかしらん? 1725 計算再開． えーと， 忘れんうちに書いとくと， 気温フィルターは降水・PPFD のそれとはカタチが違うんで， 「落す」 のは簡単ではない …… いや， 簡単か． しまった． まぁ， 気温の効果はいつでもアリってことでいいや．
• 計算まち時間は釧路大会ネット雑用など．
• PPFD 以外の気象値は「落す」とわけわからなくなる， ということのよーで． うーむ
• 計算はあまりはかどらずに 2120 研究室発． 2140 帰宅． 晩飯．
• やはり気になるんで， また 2700 ごろまで推定計算をひねってみる． しかし， いまいちですなぁ．
• [今日の運動]
  • 北大構内走 1210-1300． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0920): リンゴ．
  • 昼 (1340): モヤシ・ニラ・マイタケの あんかけスパゲッティー．
  • 晩 (2230): 米麦 0.8 合． キャベツ・タマネギ・ショウガ・エノキダケの炒めもの．

2004 年 04 月 05 日 (月)

• 0840 起床． 朝飯． コーヒー． 0920 自宅発． 晴． 0930 研究室着．
• 朝から釧路大会雑用にまきこまれる． いやはや．
• スキをみつけて苫小牧直径成長モデリングをススめる …… があまり進展ナシ．
• 北大構内ネット， mail virus にやられてしまってどうしようもない． 当方の場合， fetchmail (antispam option 設定) → procmail → spamassassin → 自作の点検用スクリプト (in Perl) という縦深防御陣地を構築してるわけだが …… いじょーな流量だよね． 呪われてる， というか．
• 雑用は片づいたけど計算はススまぬまま午前終了． 1345-1415 北大構内走． 曇． 昼飯．
• 昼飯後に Matt の壊れ ThinkPad (保証期間内) を日本国内において修繕できるかどうかを調べる． まず北大生協に問い合わせたところ， IBM PC help center に電話で修理可能性を問い合わせろ， との指示． これから先が， 悪夢のごとき展開となるわけで． だからお電話はいやなんだよねえ．
• まずなかなかつながらない． 先方はごていねいにも電話混雑状況 を図示してくれてるぐらいだ． ふーむ， Desktop 機と ThinkPad で混雑時間帯が異なるのはなぜだろう?
• で， つながったと思ったら機械対応で． IBM ごじまんの音声入力とやらはかなりダメ． で， かとーオフィスの電話端末に [PB] なるボタンがあり， これをもってすれば Pulse → Tone の切りかえができるとわかった． そこで， これで分岐を突破しようとしたんだが …… そうとうに速押ししないと回線を切断されてしまうんだよねえ． むこうの機械がしゃべり終るのを待ってては遅くて， しゃべってる最中に入力しないといけない …… やってる本人にとっては悲劇以外のなにものでもない (が， 客観的には喜劇的であるらしく， Matt はげらげらと笑ってる)．
• で， よーやく人間がでてきて対応してくれるようになった． これがまた懇切丁寧なるもので …… そんなことやってるから混雑するんでは? なぜ web 経由で必要な情報を入力させないのか， そのあたりは疑問だ まあ， PC 壊れてたら web は使えないはず， というロジックなんだろうけど， 企業ユーザーとかの多い IBM 計算機では必ずしも常に真とはいえまい …… ともかく修理できることわかった．
• そのあいだにもぢりぢりとススめてる苫小牧直径モデリング， なかなか難しいねえ． 「光合成によさそう」な気象条件を選抜していく気象フィルター (750 から 1200 とーりぐらいの可能性の中で検討している) で計算される「気象値」， そしてこういう植物生態学の分野ではよく考えられぬまま使われてる 累積気温 (はっきりいってわけわからん) なる量との兼ね合いが難しい． AIC みたいな複雑モデルに対する 罰則項を含んだモデル選択基準であっても， 「ちょっとそれヘンでは」 と思える複雑なモデルが選ばれてしまう．
• うまくいかない原因のひとつは気象要因間の交絡かも， ということで， 気象関係の変数を分離した推定も含めてやらせてみる． AIC 的な意味でよかった順にならべると，
  1. 前年気象値 + 当年気象値 + 累積気温 (1-8 月)
  2. 前年気象値 + 当年気象値
  3. 前年気象値 + 累積気温 (1-8 月)
  4. 前年気象値 + 累積気温 (1-5 月) + 累積気温 (6-8 月)
  これらは， それぞれの組み合わせの中で最良の気象フィルターを選抜させて， さらに最良の変数くみあわせのもとで得られた AIC を比較したときの 順番になってる． うーむ， いまいち複雑．
• この統計モデルが表現しようと努力しているのは 「この年はこのあたりのサイズの個体の成長がよく， それに対してこちらの年ではこのあたりが……」 ということだ． 推定結果から得られる情報としては， このパターンがもっとも重要なのかもしれない． そういう年変動する現象を， まあ多少はむりやりにこれまた年変動する気象要因あれこれと 関連づけようとしてるわけで．
• 気象フィルターに関しては， いくつか奇妙な性質があるよーで …… たとえば， 累積気温を含まぬモデルでは光量子密度 (PPFD) を重視するフィルターが選択され， 一方で累積気温が入ってる場合には PPFD の影響が最小化されたフィルターが好まれる． なぜこうなるのか， はなんとなく想像がつくんだけど． あ， でも必ずしも常にそうなるわけではないなぁ．
• 探索してる範囲ではどういう気象条件が重視されるか， だいたいめどがたってきたんで， 使用を中止していた一般化線形混合モデルをもいっぺん適用してみる …… ああ， やっぱりこの問題に関して， glmmML() による混合モデル計算はだめっぽいね． Hessian がどーのこーので計算を挫折する頻度が無視できない． そして start.sigma 依存性は気になるところだ．
• ということで， 混合ではないふつーの一般化線形モデルの推定結果， かつ前年気象値 + 当年気象値モデルに関する作図をまとめて， 簡単な報告かき & アップロード．
• 2040 研究室発． 2100 帰宅． 体重 73.6kg． 晩飯．
• [今日の運動]
  • 北大構内走 1315-1345． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0850): 米麦 0.7 合． キャベツ・タマネギ・ショウガ・エノキダケの炒めもの．
  • 昼 (1450): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・タマネギ・ショウガ・エノキダケの炒めもの．
  • 晩 (2210): うどん． キャベツ・ニンジン・エノキダケ・豆腐・煮干の味噌汁．

2004 年 04 月 06 日 (火)

• 0820 起床． うーむ． 朝飯． コーヒー． 0920 自宅発． 曇． 0930 研究室着．
• また釧路大会雑用． 公式サイトの手直しがメインなんだけど， ぜんぶ一気には片づかん． 春の陽気にさそわれてメイルもあれこれ到来しますよ． 講演要旨集だけの申し込みだの， 家族 (accompany person というらしい) のあつかいだの …… 家族の学会 (の大会) 参加って日本生態学会では まるっきり想像の埒外からの 奇襲になってるんではなかろーか．
• なぜか甲山さんの透析液うけとり …… って， ご本人いるぢゃん．
• お茶部屋の壊れ iMacG4， かとーさんが保証書もっておられましたか． ということで， またまた北大生協にメイルをだす．
• うだうだと昨日の作業内容なんぞ書いてたら， これまた時間がたってしまってですね …… などとやってたら， さっそく苫小牧から問い合わせメイルが．
• ひととーりお返事書いて …… う， いつのまにかまた釧路大会ネット雑用にとらっぷされてた． 効率化してるはずなんだけどな． みょーに時間をくってるな．
• 1430 北大構内走． 曇． 昼飯．
• 苫小牧直径成長モデリングと釧路雑用， それぞれの細かいことがちゃがちゃやってたら仕事がススまん． 何度も述べたとーり， 私にはまるちたすくは不可能．
• 北大生協から明日午後に壊れ iMacG4 (お茶部屋マシン) 引き取りにくるとの連絡メイルあり (当方からメイルで問い合わせたため)． 連絡おそい． しかも 「故障の具合など確認のためお電話してもよろしいでしょうか」 とかわけわからんこと書いてあったんだけど， めんどーになったんで「さしつかえありませんよ」 と返信して数秒後にかとーオフィスの電話がなった． 内容はたんに 「明日午後，A 棟 8F 803B 室でひきとります」 だけ， というほぼ無意味なしろもの． えーい， この …… 他人の仕事を呼出し音をして中断せしむれば 自分の問題は片づくと信じて疑わぬ お電話決戦主義者どもめ．
• いやー， 今日は釧路大会がらみのおもしろい質問というか突き上げがばしばしと． 午前の accompany person といあわせに続いて， 午後は「5500 円のうちわけ」問題だ． あまりに楽しいので一部を無断で転載してみよう:

	学会の参加は原則仕事（公費）になります。
	公費参加には参加費の内訳概算が必要です。
	お忙しい所申し訳ありませんが、5500円の内訳を教えてもらえないでしょうか。
	当然概算で結構です。
	要旨集が3000円と書いてあるのは見ました。
	残り2500円の内訳概算をお願いします。
	会場借り上げ代、当日配布資料代、文房具などなど
	よろしくお願いいたします。

• さぁー， 大会実行委員会の皆さん， どう回答されるか興味しんしん． 事情わからんヒトのために書いとくと， このあたりの金勘定はまるっきりいいかげんで 「毎年これぐらいの金額だから」 というような理由で決まってしまう根拠 free な世界なわけで． うちわけも何もありおりはべりいまそかり．
• 要旨集一冊 3000 円ナリ (これは「大会には参加しないけど要旨集だけ購入」 という場合の特別値段) ってのもむろん前例踏襲． 原価はその一割弱らしい． この値段設定は 「要旨集を高く売りつけて少しでも赤字補填」 ってハナシをどこかで聞いたよーな気がするなぁ． そんなにうまくいくものなのかしらん? 何かの間違いで 100 部も売れちゃったとしても 30 万円弱． これって参加者が予想より 50 人少なかったらそれでもう帳消し， という程度にしかなってないよね． 母分散はもっとでかいはずだ．
• 直径成長モデリングもススめたいし， そろそろアカマツ光合成補償モデリングの作文にも着手したい …… しかし釧路大会雑用を優先せんといかんよーだ． もうなんかわけわからん理由で 「ポスター賞」 が再開されることになってしまって， これが申しこみ方式にも影響するんだよねえ． ということで， このあたりの受けつけプログラムを変更せんといかん．
• えーと， すでに一般講演申し込まれた皆さんには， 後日この件に関してメイルをさしあげます， ということで． ああ， このあたりも私が準備せんといかんわけだなぁ．
• さーて， そのあたりのプログラミングにかかる前に …… 直径成長推定プログラムのたぐいを .tar.gz してアップロード， 苫小牧に 「計算プログラムとかちょっと見てもらえますか」 連絡． R プログラムの迷宮に入りこんでもらえれば， それで時間かせぎになるわけで．
• この撹乱戦法が有効なうちに， 「ポスター賞の審査対象にしますか」 選択肢の追加， それからシンポジウム発表者専用の申しこみペイジの派生， といった釧路大会プログラミングを片づけてしまって， と．
• ポスター賞， ね． どうなっちまうんでしょうかね． 少なくとも (私ではない) 何人かはかなりたいへんなことになるんでしょうね …… というか， そのへんの犠牲者がまだ決まってないあたりも大いなる難点なんだよなぁ． で， これって 「若手」 が対象なんだけど， その明確なる定義はどこにもない (私は学生会員だけにしたら， と提案したんだけど却下された)． たぶん私の倍の年齢のヒトが応募しても OK． ということで， 親切なるネット雑用掛が人々の便宜のため， 以下のような注意がきを追加したのはぢつに利他的な行為といえるだろう．
  • あなたが 「若手研究者」 なるものに該当するかどうかの 客観的な決定基準は何もありません． ご自分の責任において判断してください．
• 修整作業でばててしまったんで， 2050 研究室発． 2100 帰宅． 体重 74.2kg． なんか塩分とりすぎが体重に反映してるんでは， という気もするなぁ． 晩飯．
• 本日はなぜかしら遠来のメイルあれこれ． 修行の旅にでてる料理人氏とか， 東京湾岸の雑居房に収監されて外部との 有線連絡を禁じられてる某 SE とか ……
• [今日の運動]
  • 北大構内走 1430-1505． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0840): 米麦 0.6 合． キャベツ・ニンジン・エノキダケ・豆腐・煮干の味噌汁．
  • 昼 (1520): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・ニンジン・エノキダケ・豆腐・煮干の味噌汁．
  • 晩 (2130): 米麦 0.7 合． キャベツ・ニンジン・エノキダケ・豆腐・煮干の味噌汁．

2004 年 04 月 07 日 (水)

• 0840 起床． どうもよろしくない． 朝飯． コーヒー． 0930 自宅発． 晴． 0940 研究室着．
• 昨日の「大会参加費 5500 円のうちわけ」 問題について 竹中さん (前回の つくば大会 の実行委員長) からいろいろと助言いただく． まあ， てきとーに回答するしかないようで． そもそもいちいちこういううちわけを尋ねる 「お役所」って， じゃあ航空運賃とかの構成も航空会社に問い合わせてるのかねぇ …… 「札幌-羽田間 24000 円ということですが， このうち燃料代・空港使用料・乗務員手当・保険料・飲物や新聞など消耗品代は それぞれいくらになるんでしょうか?」 とか．
• 釧路雑用さっさと片づけてしまおう …… と思ったら， 苫小牧からのお問いあわせあれこれと． まずはこちらの応対が先だな． さて直径成長問題も次なるステップにススむわけだが …… どうやら鍋嶋さんのほうでうまくこちらのプログラムうごいたようなんで， 肩ならしも兼ねての 「こういう推定計算やってみたらどうでしょう」 というような提案してみる． いえいえ， これは「時間かせぎ」とかではなくて， ですね．
• 釧路大会の参加申込者データとむきあう． このあたりのデータの handler class に手をいれる必要がある． まずは， 軽く一般講演分類表を生成するプログラムを作って， データ構造などを思いだそうとしてみる …… ふーむ， 現在の講演者数はまだ 17 か．
• 粕谷さんから， 面白くかつ恐ろしくかつ教訓にみちたハナシを教えていただいた． 皆さん， 研究室にでいりするヒトたちには注意しましょう． 下は粕谷さんからおゆるしを得て転載するものである (太字は久保による勝手な修飾です)． 「Y 電機」おそるべし． いやー， 粕谷さんが大声だすなんてよっぽどのコトですよ．

　研究室の主力白黒プリンタ−がときどき不調なので、 エプソンのPostscriptレ−ザ−プリンタ−の安いのを買いました。 昨日、 大手量販店の１つＹ電機が納品に来ました。 知らされていた時間から30分ほど遅れて２人でやってきました。 設置して電源をつなぎました。 レ−ザ−プリンタ−は軽くなったし設置は簡単ですね。 後で考えると、 電源をどこからとったらいいかをこちらに聞かずに勝手にやるなど あやしい兆候がすでに見えていました。 私は接続すべきイ−サ−ネットケ−ブル（100base）はこれであると教えました。 ロ−カルなテスト印刷を終える（順調に終了）あたりから、 ２人はそのイ−サ−ネットケ−ブルがどこにつながっているか調べだしたので、 それはちゃんとつながっている、 何かトラブルがあるのか、 と聞くと、 特にトラブルはないといって一旦は捜索活動をやめるのですが、 私が目を離すとまた始めます。 さらに、 プリンタ−を設置した同じ部屋内のコンピュ−タ−をじろじろ見ているので、 なにをしようとしているのかと聞くと、 二人の片方がやおらウェブサ−バ−にしているマシン （Linuxで動いているＰＣ）に手をやります。 『それはサ−バ−だからさわるな！』という私の制止（かなり大声）を振り切って、 電源を切りやがりました。 『おい、 何してるんだそれはサ−バ−だって言ってるだろう』と言って詰め寄ると、 今度はあわてて電源を入れようとするので、 やめろ！（さらに大声）と制止したのですが、 かまわず電源を入れやがりました。 『客のコンピュ−タ−（しかも無関係）をやめろといっているのに、 電源オフオンしやがって、 サ−バ−いかれたら責任とってもらうぞ』と詰め寄り、 二人の片割れは一応あやまったのですが、 もう一人（電源オフオンの実行犯） は終始無言でまったくあやまりませんでした。 さらに、 二人のうちしゃべる方が言うには、 WindowsだけでないLinux,Mac,Windows混在環境 （ここの研究室はそうです）の設定はよくわからない、 とのことでした。 サ−バ−に損害があったら弁償してもらうということを、 しゃべる方と合意して、 さっさとおひきとり願いました。 あとには、 小さい液晶パネルにネットワ−ク設定をするかどうかの ダイアログが出たままのプリンタ−が残されました（すぐ動くようになりました）。 電機屋さんの配送・据え付けパ−ソンとしては、 最近にないひどさで、 しかも、 こちらがとがめてもかまわず、 無言で調べる様は、 スト−カ−を連想させました。

• 1600 から平林さんハチデータにとりくむ． 本日は R プログラム群 tutorial のつもりだったンだが …… 平林さん手動で修正したというデータ， これが以前のものとはいくつか異なる性質をもっていて， 処理プログラムの想定外の事例がぽんぽんと飛びでて， なかなかたいへんだった． これまた教訓的ではある．
  • 花 A → 花 A という移動パターンが追加された: これは「移動」の定義が以前とは変わったため． 途中の「ゼロ時間訪花」が削除された (いいのかなぁ……) ためにこういうパターンが生じた． 対策は警告だして「矢印」描かない．
  • 花 A → 花 B の移動パターンの組み合わせが増えた: 以前のデータだと最近接花に移動する確率が高かったので， 13×13 とーりの 移動パターンの一部しか使われていなかった． しかし上述の「中ぬき」をやったため， 「ジャンプ」する事例が多くなった． ここで問題となるのは summary(factor(...)) による「カウントわざ」である． じつは summary() には maxsum なるオプションがあり default 値は 100 である． つまり 100 より大きい数のカテゴリーは (Others) に分類されてしまう． この問題に関しては maxsum = 169 と指定しないといけない．
  • 不正規入力への対応: これは汎用性のある小わざで対応できない場合が多いので， たいへん面倒ながら if (...) とかで 「不具合が生じるとわかってる状況」 に対処療法するしかない (元データを書き換えない場合)． 間違い修復処理をやったときにはその旨の警告文をだす．
• ということばかりやって 3 時間がおわった． いやはやー
• 1900 すぎからカエデ論文． いよいよ投稿できるよーで． 本文は dvipdfm で PDF 化すればよい， として …… cover letter を RTF ファイルにしろってのはどゆことのなの?
• ためしに LaTeX の letter document class で手紙を書いてみて latex2rtf で RTF 化してみると …… ダメだ． レイアウトぼろぼろ．
• しょうがないんで LaTeX の article document class で \begin{minipage} & \hspace{} & \vspace{} 濫用しまくった悪しき手紙を …… これもダメ． 基本的に \hspace{} はうまくいかないよーで．
• しょうがないんで生成された RTF ファイルを Linux 用 OpenOffice 上で手作業で修正 …… くつじょく的ですなぁ． CVS 操作もしくじってしまったし． 2200 研究室発． 2220 帰宅． 体重 73.6kg． 晩飯． すでにおそいんだが眠れん．
• [今日の運動]
  • 北大構内走 1400-1435． ストレッチング．
• [今日の食卓]
  • 朝 (0900): 米麦 0.7 合． キャベツ・ニンジン・エノキダケ・豆腐・煮干の味噌汁． タマネギ・卵の炒めもの．
  • 昼 (1500): 弁当． 研究室お茶部屋． 米麦 0.7 合． タマネギ・卵の炒めもの．
  • 晩 (2320): スパゲッティー． レトルトパウチドのカレー．

2004 年 04 月 08 日 (木)

• 1000 起床． なンか終ってますなぁ． 朝飯． コーヒー． 1040 自宅発． 晴． 1050 研究室着．
• で， 仕事ススまぬまま午前終了． ますますオワってルぅー
• なんか学外とのネットつながりが遅い． そういや， 今日は virus たちの DDoS 攻撃の日だったかしらん?
• 釧路大会， 参加登録者データいじり …… とりあえずポスター賞関係の補正は終了． 試しに一般講演もうしこみ状況にそれを反映させてみる． おっと， 早めに一般講演もうしこんでくれた皆さんに 連絡せんといかんかったのだな． これはちょいあとまわし．
• 次． 釧路大会の会計チェック用のリストの生成． 会計は高田さんが担当なんだけど， この checking 下請けぎょーむはここの谷さんがやるわけで …… ということで， Perl & R つかえるヒトには CSV 形式のリストを渡せばよろしかろーと， そういうプログラムを作ってみる．
• カエデ論文， ひとまずは投稿にこぎつけたようで． めでたいことです．
• プログラム書いてる最中にもあれこれと …… 1555 ごろひとまずは終了． メイルで連絡． 少し休憩．
• 「ポスター賞」 でっちあげることが決まる前に一般講演を申し込んでしまった ヒトたちのリストを生成するプログラムを書く． まあ， どうせ私が先方とやりとりせんといかんので， 「いまさらですんませんけど， ポスター賞に応募しますか」 質問の作文． とりあえずメイリングリストに投稿． 時刻は 1655．
• 入金チェック用リスト生成について谷さんと あれこれやりとりしてるうちに， ふりがなをカタカナでいれてるヒトがいることに気づいた． めんどーなので， 受けつけ段階で強制的にカタカナ → ひらがな 変換を噛ませることにする． Perl における変換関数の例は下のとーり． ただし Perl スクリプトじたいが EUC-JP ぢゃないと 正しく動作しないと思う．

sub kata2hira # convert katakana -> hiragana
{
	my ($string) = @_;
	my $from = 'ア-ン';
	my $to   = 'あ-ん';
	return Jcode->new($string)->tr($from, $to);
}

• げ， そもそもこの CGI プログラムってば 呪われ「半角カタカナ」ちぇっくをいれとくの忘れてるよ …… これは Jcode->new($string)->h2z->euc てなかんぢで変換できる．
• で， 毎日 0200 ごろに「請求書」を生成して アップロードする設定をすませてしまう． いちおー暗号化はしておいて， と．
• 1930 ごろ上記の「ポスター賞応募しますか」メイルを出してしまう． これで， ポスター賞がらみの問題で当方で処理すべき件は， ひとまずは終了． あとは 「ポスター賞に応募します」 というのを受けつけていけばいいだけだ．
• 今日は研究がまったくできなかった． そしてさらに悪しきことに， 釧路大会雑用 4 月前半分はこれで終ったわけではない． 2030 研究室発． 2040 帰宅．
• 私は使ってないんだけど HINES web mail ってのが version かわってますますヘボくなったようで． これってやらなくてもいいのに， browswer check やったあげくに， 現在確認できてる範囲では MSIE on Windows 以外は受けつけない， という阿呆設定になってる． 終ってるねえ．
• 生活周期がなかなか復旧しない．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0920): シリアル．
  • 昼 : 食ってない．
  • 晩 (2200): 米麦 1.2 合． ホウレンソウ・タマネギ・ブナシメジのシチュー． トマト．

2004 年 04 月 09 日 (金)

• 0850 起床． 朝飯． コーヒー． 0940 自宅発． 晴． 0950 研究室着．
• 1000 研究科の入院式に出る (地環研が学部がなくて大学院のみ)． 正門わきの学術交流会館 2F にて． 「これ読んでみたら」 宿題として送られきた年輪気候学論文をよむ． 苫小牧気象データとかでモデリングやってるおかげで， こーいう論文を読んだときに以前よりは理解しやすくなっている． で， この論文の使いどころとしては， 当年の気温と樹木サイズの関係が直径成長に与える影響を見るべし， というあたりのよーで．
• 研究科のそれは 40 分ぐらいで終った． で， いったん研究科のガイダンスのあいだ， われわれは退場して生態環境科学専攻科の ガイダンスなるものが始まるのを待つ． 今回は一時間ちかく締め出されていた． そのあいだ他の人たちと雑談したり．
• 1145 ごろから専攻科 (学科のごときもの) のガイダンス． 自己紹介が終ったらさっさと退散してしまった．
• 入院式の最中に読んだ論文についてのコメント書きとか．
• 1330 より講座のガイダンス． 全員出席．
  • 公費による物品購入
  • めんどーになった薬品管理
  • ネットワークあれこれ
  • 計測機器の使いかた，管理
  などを分担して説明．
• 昼飯．
• 1500 から 1630 ごろまで 新M1 と教官の会議． 甲山さん talk 30 分ほど． ハナしてる最中に連想が拡散して蛇行していくさまは， これはじつのところ一種の芸になってるのではなかろーか …… 聞かされてる M1 たちにすれば半分もわからなかったと思うが．
• 新 M1， 4 名中 3 名は甲山さん弟子， もひとりは露崎さん弟子． 工藤さん弟子は今回はナシ． というのも新 M2 は 3 名とも工藤弟子であるからだ． それはともかく， 甲山弟子 3 人ってことは 1.5 年後に 何かたいへんなことが生起する予感が ……
• 近ごろ Matt とギロンしてると， かんぢんな場面で聞きなれない単語が使われることがある． 私の耳にはそれは 「たーぴゅれと」 と聞こえる． 「……の条件のもとで『たーぴゅれと』すると……」 というような使いかただ． 辞書で調べてもわからんので先方がヒマそうにしてるときに質問してみる (ギロン中は正確な意味がわからなくても質問してるよゆーはない)． ご本人ともども「何だろう」と懊悩したすえに， これは ``interpret'' であることが判明した． なんというか， 耳がくさってるなぁ …… ter にアクセントがつくわけか．
• 平林さんからハチ実験データの結果を格納しているデータ構造に関する質問 …… とりあえず基本的なところを説明してみる． 今日のところはこれでわかってもらえたような気がするんだけど …… documentattion は面倒だから逃げたい， とは思うけれど逃げられない気もする．
• 1810 研究室発． 1830 北 15 の中華料理屋「美食楼」で歓迎パーティー． 今日は Matt が横にいたので， うだうだと敵国語会話． 周りの大学院生は敵国語会話に慣れているよーだ． 私が大学院生だったときと比べると． 2120 帰宅． うーむ， 今日は研究どころか雑用すら進捗しなかったな．
• [今日の運動]
  • うんどう休養日になった
• [今日の食卓]
  • 朝 (0900): パン． ホウレンソウ・タマネギ・ブナシメジのシチュー．
  • 昼 (1430): 弁当． 研究室お茶部屋． パン． ホウレンソウ・タマネギ・ブナシメジのシチュー．
  • 晩 (1830): 「美食楼」 3000 円コース． たくさんの皿がでてすっかり満腹してしまった．

2004 年 04 月 10 日 (土)

• 0900 起床． 朝飯． コーヒー．
• 朝から釧路大会雑用とか． 作業効率よろしくない．
• おお， 苫小牧直径成長モデリングのキビしい一手が …… そういや， もう 4 月も 10 日ばかりすぎてしまった． 計算にまだケリはついてないわけで …… とはいえ， またしても時間かせぎ的な応じ手をさしてしまった．
• 1520 自宅発． 晴． 1530 クラーク会館地下の 理髪店 着． 二ヶ月ぶりぐらいにアタマを刈ってもらう． 1630 同発． 1635 研究室着．
• また釧路大会雑用． 一般講演 & 参加申し込みで不正入力が発生した場合の対処とか． しかし実験してみてわかったんだけど， 確認ペイジでブラウザーの警告無視して reload すると， 二重登録になるな． この段階でふせげんこともないんだろーけど， めんどーなので放置．
• 公募シンポジウム 12 コ & 自由集会 23 コのデータベイス化． こんなに応募がある (しかも全部採択される) とわかってたら最初からちゃんとした入力インターフェイス作っとけばよかったかも …… Perl わざを駆使して 4 時間ばかりかけてカタチをととのえる．
• これでまだ作業の第一段階だ． つぎにこいつをオンラインで編集できるプログラムを作らんといかん …… しかしデータいじってるうちにどういうもの作ればよいか， という見とおしはたってきた． 集会企画者たちのオンライン編集内容は即座に HTML 化できるようにしておき (というか最初から HTML にしてある)， ついには印刷組版用の LaTeX 変換につなげる． この半自動化わざ連鎖こそが我が「逃走経路」だ． 35 集会の動的データを 「間違ったらどうしよ」 「こんな面倒な手作業を」 ストレス free で管制するにはこれしかなかろう …… しかし本日はもう遅いので明日以降に． 生活周期がずれるのもイヤだし． 2140 研究室発． 2200 帰宅． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 3 × 3 回．
• [今日の食卓]
  • 朝 (0930): パン． ホウレンソウ・タマネギ・ブナシメジのシチュー．
  • 昼 (1330): 米麦 0.8 合． 昨晩の美食楼 (北 15 の中華料理屋) の残りもの， 棒々鶏 & 鶏からあげ．
  • 晩 (2230): 鶏肉・ニンジン・キュウリの炒めもの (昨晩の美食楼の残りもの) を利用した焼きそば． トマト．