ぎょーむ日誌 2004-06-30

2004 年 06 月 30 日 (水)

• 0900 起床． 朝飯． コーヒー． 0950 自宅発． 曇． 1000 研究室着．
• 昨晩というか今朝というか Perl プログラムで生成した釧路なデータを先方のサーヴァーにアップロード．
• ついでに， ということで以前から 放置しておいた釧路雑用にちょっと手をだしてしまったんだが …… これで午前中が終ってしまった． 最悪だ． オモテからは見えない部分の変更あれこれあったにせよ …… やっぱ時間かかるわ．
• 雨， か． 北大構内走でもやって気分転換したいところなんだけど．
• ちょっと次なるぎょーむを列挙してアタマの中を整理してみるかい． 研究方面では
  • 苫小牧直径成長モデリング --- 計算処理・作図の自動化をススめて次の一手に備える
  • カエデ論文の修正 --- とくに時間変動けしさりの証明
  • アカマツ論文の修正 --- まずは折れ線回帰の説明を， さっさと ……
  • ぱいぷ樹木論文 …… まだ何もススんでない!
  • 釧路大会の 自由集会 の準備， など
  なんつーか能力の限界とゆーコトバが脳内で明滅しますなぁ． 最善のばあいでもこれらが全部が片づくかどーか， というような． なんだかいきなり苫小牧な leaf miner 問題とか入ってきそうで恐い． で， ザツなやつらとしては，
  • 8 月あたまの おーぷんゆにばーしちー のポスター (対応する現在の Web 版は これ ) …… まあ， てきとーに
  • 釧路大会雑用 (1) J-STAGE な要旨集の体裁を整える (7 月上旬): これは傭兵まかせ
  • 釧路大会雑用 (2) こちらの Web な要旨集 (7 月下旬): これは印刷版をつくるついでに， 傭兵必要
  • 釧路大会雑用 (3) 印刷版要旨集 (7 月下旬?): これで釧路大会雑用という牢獄からは釈放される， 校正傭兵が必要 --- 昨日の雑用処理で要旨データはだいたい 「ほどけた」 状態にあるのでそれほど面倒にはならんだろう， と期待したい
  • 釧路大会雑用 (忘れてた) けーたい端末 web: まあこれは論文書きとかにいきづまったら， ということで
  あいかわらずスケジューリングが破綻してる． いかん．
• 苫小牧直径成長のほうを先にどうこうしたいところではあるが …… 行程進捗のよめんカエデ論文の 計算の説明変更などをまとめておいたほうがよいのかな．
• そういや， 最近は釧路大会妨害のための 多種多様な自爆メイルによるテロが少なくなったな． やっぱりイラク情勢とも何か関係があったのだろうか? くだんのちょー大型弾はまだ不発状態だったな． これはあと 2-3 週間ぐらいは放置できる．
• などと思ってたら突発的な釧路大会雑用メイル， 小さいのが二連発． J-STAGE がらみでぐだぐだ言ってきたのは， たぶん先方のサーヴァーがへぼいせいだな． もう一方については， うーむ， とアタマをひねってるうちにプログラム編成の大原さんからお電話とかですね ……
• うーむ …… J-STAGE ぐだぐだに関しては， こちらがアップロードするファイル内の各行の途中で \t\t として「なにもない列」 いれろなる不思議な作法を強要されてたんだけど (下に示してる例では講演番号のうしろ; これは変更ずみ)，

  	000028	r			0
  	000029	r			0
  	000030	r			0
  	000031	a	O1-V30		0	680
  	000032	a	P3-069		0	511
  	000033	a	O2-V03		0	684
  	000034	a	P3-142		0	584
  	000035	a	P3-143		0	585
  	000036	a	P3-034c		0	476
  

  ついつい \t としていた， というオチだった． そのうしろに重要な列がきてるんだから， エラーぐらいだせよ． とりあえずプログラムなおして make して再アップロード．
• 朝からくだらないことばかりで時間が失われてしまった． イヤな作業だからといって， こういう細部のテをぬいたらますますイヤなことになってしまう， というとても教訓的なハナシでした． イヤな雑用のときこそ いよいよ慎重に (二度とやりなおさなくてもいいように)． 時刻はすでに 1430． 昼飯．
• なぜかしらやらなくてもよい講座内雑用を． あした， Trendy セミナー あるんだけど， その reminder を半自動で発信するようにしてみる …… うーむ， ezmlm をあざむくことに失敗してしまった． だめだな．
• カエデ論文については， アンサブル平均で時間変動消去の式による証明にまちがいナシ， の確認しかできんな． 繁殖量を条件つき生涯繁殖量 ( cR0 とか Rc0 とか 書けばよいのか?) にするかどーかも， よくわからんし．
• というようなことも整理できたんで， 苫小牧直径成長問題にとりくむ． まずは作図の自動化だな． ややこしい混合モデルにしてから， 作図プログラムをそちらに対応させずに放置してきたんだよねえ ……
• で， 放置してきただけあってなかなか難しい． 以前のやつのも， こうアタマがいたくなる方法で対処してたからなぁ． 何が難しいといって結果出力を解読させないことには センを描画する式がきまらないのである． しかも 1998-2003 年の 6 本を別々に描かんといかんし．
• と思ったら敵国語メイル． ふーむ， PopWave のソースコードが見たい． そりゃ， おみせしますけどこれって Perl で書いてあって， しかも悪いことに (このプログラムだけは) なぜか日本語でコメントとか書いてるんですけどねえ ……
• とうだうだ返信してたら， 苫小牧ボスからなんだか計算のほうは進捗してないわね (久保脳内訳) ということで， ミズナラ & アオダモ & サワシバ の観測データが送られてきた． どうやら， これまでの樹種にくらべて， いろいろ難しい問題があるらしい …… 当年の光合成量が影響してるかも， とか． ひー
• で， 鍋嶋さん謹製の R による自動作図の結果がいつのまにか wiki 上にアップロードされていたので， それを拝見する …… うーむ， じつはイタヤカエデとオオモミジが別種のくせに年変動が同期してる， ってことのほうがヘンだったわけかしらん． まぁ， 1999 年ってのはノストラダムスが予言していただけあって， どの樹木にとっても「良い年」だったみたいだけど． いや， アオダモは逸脱してそうだなぁ． ミズナラとアオダモは 環孔材 か．
• しかしこのように図をつくっていただくと， 統計モデリングにはすごく役にたつわけで …… かかる年変動 & 月変動を生成しうる確率論的モデルとは， という視点でアタマの中を整理できる． すると新しいアイデアもでる．
• どうも光合成気象値と代謝気象値の影響のしかたが月によって異る， あるいは樹種によってこの両者の相対的な重要性が異る， というようなモデルで説明できる範囲にあるような …… で， 現在のモデルでは前年の光合成のほうを重視しすぎてるのかも．
• とりあえず， 光合成気象値の計算を前年の 6-9 月に限定する． 小さなバグがあってちょっとてこずってしまった …… で， これだけでアサダの成長を説明させようとすると …… やはり「暑いほどダメ」なるヘンな結果がでるけれど， 今まで得られたものとは異る． ここまでは予想どおりだ．
• どうすればよいか? 前年の影響単独で何もかも説明させようとするからこうなるわけで， 前年・当年の影響を同時に推定させればよい． 当年の影響が相対的に重要である， と示せたならばそれはたとえば「可塑性の高い樹種」 などと呼ぶことができる …… いままでやってきた計算でもこういう推定結果はだせたんだけど， 光合成気象値を計算させるところで 「自由度」 が高すぎた． これを制約すればこうった描像がより鮮明になるかもしれん．
• で， 6-7 時間かかりそうな推定計算開始． 明朝には終ってるでしょう．
• さて， 浦口さんとカエデ論文相談． いつものごとく計算は神速にて片づけられてしまっている …… やはりプロセス間相互作用はここでは「無意味」とわかった． めでたい． 問題は開花モデリングだが …… 推定結果をながめてるうちに関数型を変更すればよさそう， とわかった．
• で， 今日はへっぽこ雑用と直径成長仕事ばかりやっていた私には 宿題ばかりが積みあがっていくわけでして ……
  • 開花モデリング， glm() では解けない 最尤推定問題を解決する関数を構築する
  • 条件つき生涯繁殖量に関する appendix を書き直す
• 2350 研究室発． 2400 帰宅． 体重 73.0kg． 晩飯．
• [今日の運動]
  • うんどうしてない ……
• [今日の食卓]
  • 朝 (0910): 米麦 0.7 合． ハクサイ・ニラ・エノキダケ・豆腐のカレー．
  • 昼 (1330): 研究室お茶部屋． 食パン．
  • 晩 (2420): 米麦 1.0 合． ハクサイ・ニラ・エノキダケ・豆腐のカレー．