ぎょーむ日誌 2004-10-29

2004 年 10 月 29 日 (金)

• 0700 起床． 朝飯． コーヒー． また計算法問題にとりつかれる． 0820 自宅発． 曇． 0830 研究室着．
• 予想どーりイタヤカエデ計算には 23 時間を費したか． 残り 6 樹種はまだ計算途中だけど， とりあえずイタヤカエデ結果だけでも図にしてみる作業 …… うーむ， 作図プログラムに 細かい修正があれこれと必要． 以前の手ぬきのむくいなんだけど．
• とりあえず図はできた． 苫小牧ボス構想どーりに， 標本数が増えたことでサイズ依存性がさらに それっぽくはなった．

• しかしながら， この図をみるといくつか気になるところがある． たとえば， 「逆行」してるのは何なのだろうか． たしか直径測定用の金属バンドのつけかえが， というようなハナシを聞いたような気が …… とりあえずお尋ねメイルをだす．
• さらに， これが最良気象値セットとして選抜されたわけだが ……

         ps.prev  ps.curr       mb6      mb7      mb8
   1998 29.32520 28.47498 0.9027496 3.597631 5.870321
   1999 28.47498 27.47836 1.5465118 2.963915 5.808647
   2000 27.47836 27.18689 1.6203430 5.217963 7.027452
   2001 27.18689 30.19222 0.5989664 6.524714 4.136431
   2002 30.19222 27.40877 0.2516792 2.378106 5.018846
   2003 27.40877 29.60587 1.4684730 1.525402 3.319788
   2004 29.60587  0.00000 0.7690110 0.000000 0.000000
  

  当年の天気に依存する代謝気象値は おそらく成長量の集団平均によく対応してるように見えるんだが， 昨年の天気で決まる光合成気象値のほうは必ずしも単純でない ……
• これは何を意味するか， というと …… ここしばらくいろいろ検討している方式の 推定計算高速化がかなり難しい， ということだ． 代謝気象値と光合成気象値を分離して「序列化」できん． たとえ近似的， と制約をゆるめたとしても．
• つまりこれら二つの気象値はセットのまま評価せんといかんわけで …… どうしたもんかな． アイデアなし．
• つまったので， 雑用始末． 来週の出張申請． ひどく腐れたわーど書類なので， A801 の共用ゐんどーづ機で処理する． どうも印刷うまくできぬと思ったら …… 用紙指定だけでなく， 「倍率指定」とやらも A4 にしないといけないのか． 研究科事務で作った原版が A4 を逸脱しているためだろう ケイ線遊びに夢中になってるうちに 紙面をはみだしてしまったに違いない．
• 苫小牧ボスより連絡いただき， 「逆行」問題の一部は R プログラミングの型変換が原因かも， というような． R の「型」 (とゆーか class) って実際のところ， 便利かつ不便なんだよねえ …… それはともかく， 2004 年 4 月に一部の樹木個体の直径バンドを つけなおしたそーで． だとすると「逆行」あっても不思議でないな．
• 1305-1340 北大構内走． 晴． 気温， やや高め． 昼飯．
• 事務どたばたにまきこまれる． 出張書類苦闘． 年末調整あれこれ， 会計係らしからぬ手ぬきでまったく意味不明になってる． 問い合わせ． 先日の毎木調査の謝金問題． 金の出どころをどうすべきかわからず， 自分のところから出すことに． ひとくぎりついたら時刻はすでに 1550． やれやれ．
• jeconet に流れた情報によると， 大阪大会 でのシンポジウム・自由集会の申しこみは「まだまだ」 のよーで．
• 会計事務処理問題解決． くだんの給料明細システムの「横」にあった． 年末調整はほぼ完全にオンライン化か． 現時点では印刷・押印が要求されてるけど， 将来的にはこれも不要になるかも． けっこうカネのかかっているシステムで (ケイ線再現にプログラマーたちの苦闘のあとが読み取れ， また「統合オンライン事務処理」としてはまだまだであるけれど)， 北大内の他のみぢめな事務系ペイジとは一線を画している．
• オンラインな北大サイトで， 一番便利で使いやすいのは図書館． 大学内組織ではないけど北大生協もそこそこに努力している． 上記会計システムをのぞく事務系はかなり劣悪． 一番ひどいのは言うまでもなく情報基盤センター (つまり北大ネットの管理者ども) である．
• とりあえず， 会計関係の書類を始末． 時刻は 1615． 昼飯後はこんなことばっかやってるな．
• 苫小牧ボスから新しい成長データ送っていただく． 「逆行」というか DBH マイナス成長はほとんどないはず， とのこと． 念のため全樹種の自動作図プログラムに放りこんでみる …… たしかに「逆行」はほとんど無くなった (目でみてわかる範囲では)． 1720 推定計算再開． 旧来の計算方式のままなので 80 時間ぐらいかかりそう．
• お茶部屋植物雑談． ゴボウ (Arctium lappa) はアザミ科のは二年草で こんな地上部 つけると初めて知った． 日本には自生しないと書いてあるけど， 北大内にはたくさん生えてるらしい． 昔， 工藤さんはこれをそのへんから掘りだし食卓に供していたそーだ．
• かとー先生が， MacOS X 版の R-2.0.0 は良いと絶賛されている． たしかにユーザーインターフェイスは使いやすそうに見える． MacOS X なヒトたちはインストールしましょう． 私も Linux 版 2.0.0 を数日前から使ってるけど (中間さんのサイト からダウンロードした Vine-3.x 用)， lazy load 機能なんかのおかげで起動が速くていいかんぢだ．
• データ解析言語Rによる統計的プログラミング: 第2回　機能プログラミングとデータ操作 統計学用語にふなれな翻訳者が苦闘していておもしろい．
• 1930 研究室発． 2000 帰宅． 体重 73.6kg．
• 律義なる苫小牧ボスより， さらに完ペキを期したデータ送っていただく． A 棟 8F にもぐりこんで， また計算再開． この方式では 80 時間とかを要するので …… といったメイルやりとりをやってるうちに， また少しマシな策を考えつく． もちろん， くだんの高速化 (手ぬきの別名) に関する．
• つまり， いいかげんな測度をもちいた高速化用選抜なんぞは役にたちそうにない， というのであれば …… 単にすべての気象値セットにたいして事前に glm.nb(growth6 ~ ps.prev + mb6) すればいいだけじゃないか． いや， 厳密には glm.nb() で評価できる負の二項分布モデルとはちょっとズレてるわけだが (しかし glm.nb() で十分かもしれん， というのは検討の余地アリ)． この準備段階にあたる序列化計算のためだけに， おそらく 1458 気象値セットならば 10 分以上を費すだろーけど， 80 時間に比べればぜんぜんたいしたことないわけで．
• [今日の運動]
  • 北大構内走 1305-1440． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0830): 米麦 0.7 合． タマネギ・豆腐のカレー．
  • 昼 (1400): 研究室お茶部屋． 米麦 0.8 合． タマネギ・豆腐のカレー．
  • 晩 (2200): 米麦 1.0 合． タマネギ・豆腐のカレー．