ぎょーむ日誌 2001-11-20

2001 年 11 月 20 日 (火)

• 0745 起床． う． 寝坊した．
• 朝飯・弁当の準備． 朝飯．
• 0810 自宅発． 晴れ． 特急に乗って 0822 京急平和島発． 普通に乗り換えて 0850 京急上大岡発． 0905 研究所着．
• 今朝がたまでやってた並び換え・均等分割・頻度分布計算は あまり当方が先走ってもしょうがないので， split_dist.pl と split_plot.pl というふたつの Perl スクリプトにコメントいれたり， Makefile 書いたりして浦口さんに引き渡してしまう準備をする．
• このへんの連絡あれこれ． いつも 24 樹種 (小川) とか 30 樹種以上 (インドネシア森林) とか 500 樹種近く (Pasoh) とか おおざっぱな直径成長解析しかやってないので， 一樹種を丁寧に見るのはなかなかおもしろい． 全樹種おなじ数式 を使ってパラメーター推定などしたら 検出できないような， びみょーなサイズ依存性とか明るさ依存性とかね． というか， 本来こういうふうにじっくり見るべきなんだろう．
• Perl スクリプトを hosho に置いたり メイルで「うーん，よくわかりません」と書いたりしたせいか， 浦口さんお悩みモードに突入したもよう． 今回は時間稼ぎ成功か? はたまた呆れられただけか．
• 先方が長考してる間に， こちらは (じつはあまり使ったことない) ベータ分布で遊ぶ． こういうややこしい確率分布は R を使って計算するのが 手軽である． といいつつ， ふだん R を使ってないので いろいろ調べていて時間かかった．

[ベータ分布] こういう分布です． x ∈ [0,1]

• 上のような図を生成する R スクリプトは こういうかんぢ． 実行するには R --no-save < スクリプト名 とか． Rplots.ps というファイルに結果出力．

	x <- seq(0, 1, length = 101);
	par(bg = "#d0d0d0", fig = c(0.25, 0.75, 0.15, 0.85))
	for (p in seq(1, 19, length = 10)) {
		plot(
			x, dbeta(x, p, 5),
			xlim = range(0:1), ylim = range(0:5), type = "l"
		);
		par( new = TRUE )
	}
	title( main = "Beta Distribution", font.main=3 )

• 漬物しかない昼飯の弁当など． いつのまにか 2ch のプログラミング関係で void 御大が猖獗をきわめるようになっている (うん? 単体 (のハズ) なんだけど猖獗という表現が似合ってる)． 実は告白すると， 先日 ご著書 を買ってしまっていたり． これは TSUTAYA のわりには計算機関係に気合いの入ってる 京急大森町横のいつも行く店で購入． この本って C 初心者というよりむしろ void watcher に買われてるじゃなかろうか (私は watcher とゆーより圭子サマふぁん)． とはいえ， ちゃんとした内容なので C を勉強したいヒトにおススめできます．
• 水循環の栗田さんから画像処理うんぬんの質問あったので， あれこれまとめて ImageMagick で解決してみては， と簡単な使いかたを実演しつつこれをおススめしておく． 画像処理を専門とするヒトたちでも， このすごく便利なコマンドラインなツール群を 意外と使っておられなかったりする． かくもエラそうに述べてる私は「ぎょーむ日誌」に貼りつける へっぽこ画像ファイルの変換ぐらいにしか使ってないんだが．
• ちょっと統計学関係のことを調べてから， 「ニセ」アカマツ再建作業に取り組もうとしてみる．
• ここしばらく独房群内でも (体感的に) かなり寒かったんだけど， 今日は体調も回復しつつあるのか， セーター着てればまったく寒くない状態．
• 「ニセ」アカマツに葉っぱ情報を記憶させる作業の続き． 動作確認がなかなかはかどらない． あちこち改造しかけのまま
• 考えモードの浦口さんから難しい統計学質問あれこれ． たとえば， Kolmogorov-Smirnov の二標本検定において， 私が調べた Sokal-Rolph 訳本は標本数が同じでなくてかまわないと述べ， 浦口さんが調べた石居本においては小標本においては 数をそろえろと書いてるそうだ． えーい． 著明な教科書を信じるのがイヤになってきたので， 最尤推定法+モデル選択で安直にケリをつけたら， とか ranodmization で信頼区間計算できるでしょう， とか悪の道にゆーわくしてみる．
• カエデ太り問題はよくわからぬまま． 1750 研究所発． 普通に乗って 1800 京急杉田発． 快特に乗り換えて 1810 京急上大岡発． 普通に乗り換えて 1830 京急川崎発． 1910 帰宅．
• 晩飯食ってから久しぶりに風呂に入ってさっぱりする．
• しかし仕事すすまん．
• 今日の食卓
  • 朝 (0750): 米 0.6 合． キムチ． さすがに今朝は料理してるヒマなし．
  • 昼 (1220): 弁当． 米 0.8 合． 朝と同じ． うう．
  • 晩 (1940): 米 0.6 合． ダイコン・ニンジン・シイタケの味噌汁．