ぎょーむ日誌 2001-10-22

2001 年 10 月 22 日 (月)

• 0730 起床． ねむい．
• 朝飯・弁当の準備． 朝飯． シャワー．
• 0800 自宅発． 曇り． 特急に乗って 0810 京急平和島発． 普通に乗り換えて 0830 京急上大岡発． 0845 研究所着．
• 昨日の「ノンパラメトリック・複数母集団」問題を 脳内からなかなか消去できない． もし問題がパラメトリックであるなら， 母集団がいくつあろうが近似的な解法ならば 何か構成できるような気がしてきた …… たとえばベイズ推定みたいなものを持ち出す， とかね．
• しかしノンパラメトリックで複数母集団となると …… 何ができるんだろう． それらの確率分布の性質として利用可能そうなのは， ハシからハシまで積分すれば 1 になります， ってことだけで他は何も使えそうにない．
• ノンパラメトリックな分布の推定 …… うーむ， ノンパラメトリック回帰， になってしまうんだろうか．
  1. ふたつの標本集団が異なる母集団から 得られたことを示す (すでにこの段階でいいかげんさが混入 …… 避けがたく)． ノンパラメトリックな回帰+クロスヴァリデイション?
  2. それぞれの推定分布を計算． ノンパラメトリックな回帰+クロスヴァリデイション? 次に検定するなら， この段階でその帰無仮説を 勘案しつつ分布を推定してみせないといけない．
  3. それらの分布の統計量 S に関する なんらかの検定を構成する
  …… たぶん， 本職の統計学者からみると穴だらけだろうが， これでしろーとをダマせる程度の枠組にはなってるだろうか?
• 上の 1. と 2. の推定において クロスヴァリデイションやっちまったのであれば， 3. の検定は切ってもいいところだと， 私は考える …… しかし粕谷さんメイルでは， 何と言いますか分布そのものだけ見るのは不十分で (あるいは必ずしもいつも最重要というわけではなくて)， むしろ分布から得られる記述統計量こそを重視する立場 …… のようなニュアンスだったんだよね．
• たしかに， 分布は斟酌しないのに記述統計量となるとありがたがる 生態学者はたくさんいる (粕谷さん自身は分布もすごく重視するヒトであるが) …… じゃあ， もし仮に上の 1-3 の手続きを経て 「異なるふたつのノンパラメトリックな確率分布から それぞれ得られた標本集団間において， (たとえば) ある順位統計量 RS は等しい， という帰無仮説は棄却できない (あるいは，できる)」 という検定結果が得られたときに， 彼らはどういうふうにその結果を利用するわけ?
• …… とアタマの中の思いつきを「ぎょーむ日誌」に core dump しておいて， 「ニセ」アカマツに取り組むのであった．
• …… と思ったら北大の甲山さんから計算ぎょーむの発注． インドネシアから送られてきた Sebangau の森林調査区の成長解析ね …… こういう国内・国外の計算下請けぎょーむは 修士のころから多種多様にこなしてるんだよなぁ． いまのところ金銭報酬をふんだくるカラクリを考えついてないので， まだ国際びぢねすに発展してないだけである．
• たいしたことやるわけでなく， 前回 Perl スクリプト+最尤推定プログラムの連携ワザで 作っておいた自動解析システムを手直しして新しいデータを ほうりこむだけのことなんだが …… Perl なフィルターがデータ異常の警告． たぶん今回から新しく死亡個体もカウントしてるんだろうな， これは．
• そのあたりの質問メイルだして， 返信待ちつつ昼飯． ふーん， 弁当の飯を炊くときにうまいぐあいに 酢をまぜると 味を損なうことなく菌群集の増殖を阻害できるんだ．
• 甲山さんからさっそく返信いただき， あれらはやはり死亡個体とのこと． フィルター直して，

	make filtering
	make estimation
	make plot

でおしまい (さらに打鍵の手間を惜しむなら make all とか)． 1300 関連ファイルを FTP サイトにおいて 熱帯林ぎょーむ終了．
• 1500 からここの事務主催の説明会というのを聞くべく大会議室にいく． 給料はどう決まるとかどうやったらクビになるか， というハナシである．
• いやー， うかつなことに説明きいてはじめて知ったんだけど ポスドク傭兵の雇用期間は 3 年なんだね． 毎年毎年ぎょーせきチェックされるから (つまり来年からいきなりクビ， というのがありうる)， とくに期限はないのかと思ってた． 再応募するのはかまわんとのこと (この点が日本のふつーのポスドク傭いとは違うかも)． あるいはふろんちあ抜けられたら困るヒトは 指名再雇用とか …… しかし 生態系領域 の下っぱである ここあーきてくちゃー組の場合， グループごとお取り潰しになってもまったく不思議はないね． すごく異端なもんで．
• …… さてさて， いよいよのときは生態学計算ぎょーむの国際びぢねす化でも 考えねばならんのかと思いつつ独房にもどる． どうすれば商売としてなりたつだろうか． 需要がないわけではない， はずなんだが．
• 甲山さんから Sebangau の新プロットの成長ヘンじゃないの， という御指摘メイル …… しかし元データみても新しく追加されたプロットがヘンなんだよね． ということで赤道の向うに問い合わせメイル出していただくことに．
• 「ニセ」アカマツの simulation procedure 書こうとしたら， なぜか書けない …… よく考えてみると二ヶ月放置してたのでホントに 忘れてしまっているのだ． Perl スクリプトと原稿をよくみなおす．
• 1750 研究所発． 雨． 普通に乗って 1803 京急杉田発． 快特に乗り換えて 1810 京急上大岡発． 普通に乗り換えて 1830 京急川崎発． 1900 帰宅．
• またしてもうっかり AC アダプターを独房内に残置してしまって， ThinkPad の自宅における稼働時間は 1 時間ちょっとに． 悪いことは重なるもので， 昨晩 apt-get でうっかり rp-pppoe を更新してしまったら， ADSL 接続できなくなってしまった． やれやれ． 明日ゆっくり直すしかないなぁ． ということで今晩はメイル読み書きもできなければ 「ぎょーむ日誌」アップロードもできないわけである． もはやここまで．
• 今日の食卓
  • 朝 (0740): 米 0.6 合． 豆腐・タマネギ・シイタケの炒めもの．
  • 昼 (1230): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1910): 米 0.7 合． 朝と同じ．