ぎょーむ日誌 2001-06-15

2001 年 06 月 15 日 (金)

• 0610 起床． やや眠い．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー． ひげそり．
• 0700 自宅発． 今日も雨． 特急に乗って 0709 京急平和島発． 普通に乗り換えて 0733 京急上大岡発． 0750 研究所着．
• さっそく， 下請けぎょーむの続き． Makefile がへっぽこで， 自動処理行程の推定の前で止まっている． 直してそこから処理続行． これは 30 分ほどで片づく． しかし， パラメーターの相図を描く， というようなあまり意味の無いことをやっても …… やはりたいして違いはない．
• 2 パラメーター/3 パラメーターモデル の尤度を比較してモデル選択できるような Perl スクリプトを書く． 一時間ほどかかる．
• やってみると，

	元データ
	# alt_selected: 186
	# null_selected: 285

	xy 無作為化
	# alt_selected: 92
	# null_selected: 379

てなかんぢで， 3 パラメーターモデル (サイズと開空度が影響) が「モデル選択」されるのは 471 樹種のうち 186 樹種にすぎない (ただしこの数値は三次元林冠モデルの形状に 強く依存している)． とはいえ， 樹木の位置を無作為化 (xy 位置のランダム交換) した場合は「選択」される樹種数はさらに半減するから， Pasoh にいる樹木のうち何割かは かなり「明るさ」に依存していそう， と言えそうだ．
• あ， 「モデル選択」 の基準として使ってるのは AIC (Akaike's Information Criteria) です． 尤度比検定とどっちがマシか ……
• 午前中はあとはひたすら作図で終る． Pasoh ペイジにアップロード．
• (更新情報) Pasoh 解析ペイジ をぼちぼちと更新してます． (2001.06.15)
• 昼飯の弁当食ってから郵便物を見に行くと …… うーむ， NASDA 広報め， また無意味なマジナイを ……

[打ち上げ祈りシール] 「この夏 H-IIA 初号機打ち上げに向けて」 と書かれたシールが 30 枚． 一緒についてきた「ぎょーむ連絡」 によると， これを NASDA 名刺にぺたぺたはって お友達に配れ， ということらしい．

まぁ， 気合い入るのもわかるケドさ． しかし， ロケット打ち上げというバクチのオッズ どんどんあおってだいじょーぶなのか? そりゃあ， 成功すれば盛り上がるだろうが ……
• また停電か． どうしていつも全館停電なんだ．
• さて午後は 2030 までかかってさまざまな作図作業をしては， 上記 Pasoh ペイジにアップロードし， 作図にはコダわりのある甲山さんの修整指導をうける， のくり返しであった．
• いろいろと重要な指摘だの， 作図ワザの開発などあったりしたわけだが， なにしろ「ぎょーむ日誌」書いてるヒマもないほどの 作業量だったので， その多くは闇の中に永久に消えてしまったのである． ふだん私がでれでれとばかりしてるのは， まさにこういうことにならないため， と期せずして正当化されてしまったではないか．
• ひとつだけかすかに覚えている gnuplot ワザについて． gnuplot の出力が Postscript (PS) をもって 基本とすべきである． これが一番マトモであり， 論文だの口頭発表だのに使いたい品質が実現するからである． しかし set term postscript の制動はだいたいにおいて面倒といってよく， このあたりでキレてしまうと 長谷川宣教師や牧さんのように PS に直接さわるような ばてれん黒魔術 の世界を邁進していくことになる．
• つね日頃から常識というものを重んじている私は もっとふつーな解決策しか思いつかないのである． たとえば， 甲山さんから PS から PNG に変更するのはケッコウだが， もう少しその解像度をあげてくれ， という請求があったとしよう． やってみればわかるが， これはそれほど簡単ではない．
• gnuplot で PS ファイルを生成， ImageMagick の convert で PS から PNG に変換． これは変わらないところだ． しかし convert コマンドでは単なる 「拡大」はできても それは解像度あげていることにはならない．
• 正解のひとつは， gnuplot と convert の あいだに psutils の psresize を噛ませる， というものである．
  1. gnuplot の set term post landscape color "Helvetica" 18 というような指定で fig.ps という PS ファイルを生成．
  2. psresize -Pa4 -pa3 fig.ps fig_large.ps というふうに A4 から A3 にサイズ変更． これによって解像度が約 1.4 倍増える．
  3. 最後は以下のように convert する．
     convert -page A3 -filter Point -crop 0x0+10+10 fig_large.ps fig.png
• 上のワザをもうちょっと工夫すればいろいろな解像度にも 対応できそう． あるいは convert -page ... を工夫すれば psresize ナシでも問題解決できるのかもしれない．
• どうみてもまとまりそうになかった Pasoh 解析について， 甲山さんの工夫で ``Ashton Forest'' v.s. ``Species Hubbell Forest'' v.s. ``Spatial Hubbell Forest'' という新しい文脈でまとめられることになる． 森林内樹種たちの機能的分化を主張する Peter Ashton と そのランダム性というか「中立」性を重視する Steve Hubbell， それぞれの立場を借用したウマい整理である． よくそんなことを思いつくモンだ， と文献をほとんど読まない私などは 感銘うけたのであった．
• 進捗したのかしてないのかよくわからぬまま撤退． 2030 研究所発． 雨． 普通に乗って 2043 京急杉田発． 快特に乗り換えて 2050 京急上大岡発． 普通に乗り換えて 2110 京急川崎発． 2140 帰宅．
• さっさと寝てしまえばよいのに， なぜか nethack を 0100 までやってしまう． うーむ?
• 今日の食卓
  • 朝 (0625): 米 0.7 合． ワカメ飯． キャベツ・ニンジン・タマネギ・シイタケの 野菜炒め．
  • 昼 (1240): 弁当． 米 0.7 合． 朝と同じ．
  • 晩 (2210): スパゲッティー． 朝と同じ．