ぎょーむ日誌 2001-07-02

2001 年 07 月 02 日 (月)

• 0700 起床． 寝たのが 0030 すぎ． なかなか 0600 起床の生活に戻らん．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー． ひげそり．
• 0750 自宅発． 晴れ． 特急に乗って 0803 京急平和島発． 普通に乗り換えて 0830 京急上大岡発． 0845 研究所着． この時刻ですでにすごく暑い． さらに京急の車内は冷房がキツいので (かなり寒い)， 外はよけー暑く感じる．
• 「ニセ」アカマツ改良の続き． 「ニセ」枯れ上りとはいえ， みょーな小細工 (制約の多い親子茎関係とか 「減衰率」とか) は廃しても， 個体として バクハツもしなければ速やかすぎる自滅もしないんで， この路線で押してみよう．
• 現在の茎の親子関係は以下のようになっている． それぞれの「ニセ」枝 (ホントの枝なんて無いといったら無いんだっ) 先端部はこういう構成である．

シミュレイション過程のほうでは，
1. main の子供は親の複製である (したがって， ある rank 内の茎の数に比例する)．
2. sub も当年生茎の複製であるが， 所属する rank をランダムに選ぶ (したがって， ある rank 内の茎数とは無関係 …… なお， 「ニセ」枝形成のため， ランダムに親となる茎を選んでいる)．
という処理をしている． ところで， アカマツ解体データから個体を でっちあげ的に再構成するときに， すでに 1. のような処理はやってるんだけど (といってもかなり怪しい)， 2. に関してはいきなり親茎をランダム選択してたので， まずは rank を選ぶように改める． つまりこのあたりの手続きを できるだけ矛盾無いようにみせかける努力である．
• 「ニセ」枯れ上りを導入することで， いまのパラメーターセットでは 「個体自滅」という現象が発生しなくなった． rank_number が相対的に低い， つまり樹冠下部の「ニセ」枝たちをばさばさと 惜しげなく切り捨てていることで， 資源を浪費するだけの「高齢茎」の発生を 最小限に抑えているためだ …… とりあえず小林さんに 簡単な現況報告メイルなど出してみる．
• しかし， この all Perl な「ニセ」アカマツ， 計算機のメモリをバカ食いするなぁ． 茎の数が一万から二万で破綻する． つまり， これぐらいになると BookPC の 192 MB の RAM を使い切って キャッシングが始まり実行速度が急減速する． Perl のリファレンスカウンターは ちゃんと動いていて， 不要になった領域は開放されてはいるようだけど．
• やはり all Perl な樹木シミュレイターは難しかったか …… いや， データ保存のやり方の問題か， などと考えつつ昼飯食う．
• かな漢字変換 Wnn7 の案内． 7 月 7 日発売． いまさら完全 free の Canna から 有料製品の Wnn には戻れんなぁ (Canna は相変わらずとぼけたやつですが) …… と思いつつ Omronsoft のペイジ をのぞいてみると …… へー， Wnn って近ごろは携帯電話とかにも使われてるんだ． 京急の中でメイルを打ってる 高校生の女の子たちも Wnn 使いなのか!
• 「ニセ」アカマツ，

  (生産量) ∝ (生産器官の重さ)

  とすると茎数の増加ますます著しい． どうでもいいんだが， 三次元図に描いたときに「ニセ」枝の トラ刈りぶりがますます目立つ． やはり

  (生産量) ∝ (生産器官の重さ)^α

  とでもして， α を 0.5-0.9 ぐらいにしてみるか．
• とうぜん， 対応するパラメーター (比例係数) 推定も この生産モデルを使う．
• しかし残念ながら， このαそれ自身はデータから決めようがない値である． というか， 当然ながら実際の樹木では生産力はこの部分は 環境 (光環境とか) 依存だろう． ま， α小さいほうが計算速く終ってラクではある．
• コマンドラインオプションで α指定できるようにして実験 …… なるほど， α小さくすると主幹以外の枝がほとんど落ちるけど， 茎の数は増える …… これは新規生成茎の サイズ分布がどんどん小さくなるからだな …… 観測データのある三年生茎を固定した母集団にするか．
• としてみるt， 重い茎が増えて小さな茎が減ったから 計算は速くなった． 茎のテンプレイトが常に固定されてるので， たぶんよりロバストになったはず． それはいいんだけど， しかしα=0.5 の計算結果はかなりみじめだ．

[まだ生きてます] 30 年生． そうか， 生産器官が小さいときには αが小さいほうが効率いいのか．

• たとえば， 処理区の第 2213 号ほんものアカマツのデータから作った 「ニセ」アカマツ個体だと α=0.8 では茎数が 2000 ちょっとで安定 (収支がバランスしてしまう)， メモリも食いつぶさない計算も速い (30 年生まで 80 秒) …… とういことで， このあたりを中心に調べればよいか．
• ちょっと全コードを Postscript ファイルにしてみると …… 全部で 42 ペイジ? psnup コマンドを使って 一枚に 4 ペイジ印刷する． うーん， いつのまにか合計 1800 行以上に増えている． でも， モデルが複雑化したわけではないからいーや． 時刻は 1515．
• 印刷した all Perl な樹木シミュレイターコードながめる． 解体データから個体を再構築するのと， 将来の状態をシミュレイトする部分が混在していて …… 分量としては， 「ニセ」枝形成の三次元関係も少なくない． 「ニセ」とはいえ， 三次元作図は動作チェックには欠かせないものだ． ふう， なかなか面倒なものを作ったなぁ．
• X 上で動作する Emacs-20.7 のうっとうしい スクロールバーを削除する方法を試行錯誤． おそらく $HOME/.emacs などに

	(set-scroll-bar-mode nil)

と書けばよいのである．
• 北大<鳳翔>用のバックアップ HDD が 壊れたんだけど， かとー大先生が新しいドライヴにつけかえてくれた． バックアップのやり方もいままでよりスマートになる．
• Princeton から戻った甲山さんからメイルいただく． 彼の地の森林動態モデリングまにあっくすな ヒトたちにはウケたようで， よかったよかった．
• プログラムの挙動調査も兼ねて， あれこれと条件をかえて計算． 「ニセ」アカマツたちの成長する様子を見る． 1750 研究所発． 普通に乗って 1805 京急杉田発． 快特に乗り換えて 1810 京急上大岡発． 普通に乗り換えて 1830 京急川崎発． 1905 帰宅．
• 晩飯食ってでれでれする． う． 悪しき生活習慣．
• (現在判明している) 残された問題， パラメーター推定問題をどう処分しようか， 迷う．
• 今日の食卓
  • 朝 (0715): 米 0.6 合． キャベツ・ニンジン・エノキダケの味噌炒め． キュウリ．
  • 昼 (1240): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1930): 米 0.8 合． 朝の残り．