ぎょーむ日誌 2001-11-(11-20)

2001 年 11 月 11 日 (日)

• 1020 起床．
• 昨日書いてた，

	確率  1/2 性比 1/1 男
	確率  1/4 性比 1/2 女男
	確率  1/8 性比 1/3 女女男
	確率 1/16 性比 1/4 女女女男
	……

に関して， よゆーの博士論文執筆を楽しんでおられる 長谷川宣教師殿 から「子供の性比は 1:1 になるんじゃないの」 という御指摘をいただく (月曜日になって 竹中さん からも同様の御指摘)．
• つまり， 一家族あたり男の子供の期待値が常にであるのに対して， 一家族あたりの子供の期待値は 1/2 + 2/4 + 3/8 + ... = 2 となるので， 性比は 1:1 ではないか， ということである． これは 遺伝学概説 (JF Craw著、木村資生ら訳) にも紹介されてる問題だそうで．
• たしかに集団の子供全体の性比はこのように 1:1 になる． では， 私が計算した 7:3 とかいう怪しげな数値は何なのだろうか? これはぢつは 「集団中から一家族をランダムに選んだときに， その子供たちの性比の集団全体の期待値」 がこの数値になる， という …… まぁ， あまり使いみちのない無意味な計算結果なのであった．
• なんか週末ばてなので何もやる気しない． NetHack もやっていない． 前にも書いたけどあれは知力・気力・体力をかなり要求されるんで …… 掃除とかやらんといかんのだけどなぁ． さすがに食糧切れには対処しなければならず， ATM 閉鎖時間前に現金をいくばくか下ろして安い米を探索する． 新潟県産「ひとめぼれ」精米日 11 月 9 日 10Kg 3200 円． このへんで妥協するか．
• 今日の食卓
  • 朝 (1040): パンケイキ． コーヒー．
  • 昼 (1430): ニンニク茎・ニラ・シイタケの焼きそば．
  • 晩 (2000): 米 1.0 合． ハクサイ・レタス・モヤシ・タマネギ・豆腐・コンニャクの いためもの．

2001 年 11 月 12 日 (月)

• 0700 起床． 眠い． ここ二日ほどでれでれしてたんだけどなぁ．
• 朝飯・弁当の準備． 朝飯． コーヒー． シャワー．
• 0750 自宅発． 雨． 特急に乗って 0801 京急平和島発． 普通に乗り換えて 0830 京急上大岡発． 0845 研究所着．
• くだんの「48 時間以内に返せ」校正刷りの真の deadline が接近しつつあるので， 朝からいきなり集中してそれにとりくむ． やってみたら 2 時間半ばかりで終った． 速達便にして投函．
• 昨日はついに一行も書けなかった「ぎょーむ日誌」を書いてみたり． とかやってると午前中は終った．
• 昼飯食って， このへんの文書をでれでれと書く．
• 午後から発芽確率の最尤推定に取り組む． 温度によってこの確率がふらふらと変わる． こういうしんどい推定には mlfitting を使うほかない． さらに処理がふたとーりある． このあたりは mlfitting でもさばききれないので， Perl なんかによる統括が必要だな．
• 1600 ひととーりできた．

縦軸は確率で横軸は温度． 点で示してるのは「参考のために」計算した 単純な比率 (発芽数/全種子数) である． mlfitting はかかる比率なんぞは斟酌しないで， 各処理温度における 一個一個の種子の発芽した/しなかったの集積である 二項分布に合致するような 曲線をでっちあげるようになっている．
• いつものごとく人間に対する説明はむずかしいもので． ごくごく簡単な説明だけなのに 30 分ほどかけて説明メイルを書いて送信．
• 京大農学部で Perl を使いこなしておられる牧さんから Date::Calc という CPAN モジュールが 日数計算とかに便利ですよ情報いただく． さっそく試運転 …… おお， これでよーやく system "date 複雑怪奇なオプション" というケガれた使いかたから免れるかも ……
• さて明日からは「ニセ」アカマツに専念できるかと思いきや …… げ， 事務からわけのわからん書類の堤出予告が． 詳細まだ不明． 1810 研究所発． 普通に乗って 1820 京急杉田発． 快特に乗り換えて 1830 京急上大岡発． 普通に乗り換えて 1850 京急川崎発． 1910 帰宅．
• もと過激派の某宣教師殿からよくわからぬテロ予告があったんだけど， 早めに寝てしまう． しかし夜中にめざめる．
• 今日の食卓
  • 朝 (0710): 米 0.6 合． ハクサイ・レタス・モヤシ・タマネギ・豆腐・コンニャクの いためもの．
  • 昼 (1220): 弁当． 米 0.7 合． 朝と同じ．
  • 晩 (1930): 米 0.7 合． 朝と同じ．

2001 年 11 月 13 日 (火)

• 0700 起床．
• 朝飯・弁当の準備． ラジオのニュース聞いてると …… なんかまた飛行機が落ちてるな． 今回はエアバス． 昨晩のテロ予告はこれだったのか? 朝飯． コーヒー． シャワー．
• 0745 自宅発． 晴れ． 朝から晴れてるのはしばらくぶりのような気がする． 特急に乗って 0754 京急平和島発． 普通に乗り換えて 0820 京急上大岡発． 0835 研究所着．
• これまた昨晩予告のあった事務の面倒書類． 金曜日まで， か． また腐れ M$-Word なんかで準備してるし． VBA virus にやられてしまえというかんぢ． PDF ファイルあるけどこれは加工できん． LaTeX で模造するか．
• 日本お役人が愛してやまない無意味罫線タテヨコに引かれた 阿呆書式を tabular でマネしてみる． 贋作専門家の私にもこれはつまらん作業である． ほんと罫線って無意味なばかりでなく， むしろ美観と機能を損なってると思う． しかしながら， そのようなふつーの判断を連中に要求するのは酷といえるだろう．
• ばかばかしいのだけど定規でワクの大きさなど測定し， 1 時間ちょいもかけてだいたい似たような見てくれの偽物を作成． これを埋めていく．
• まだ片づいてないけど昼飯． あーん? スーパーカミオカンデ も壊れちゃったのか? これも某宣教師殿のテロ活動だろうか …… と思っていたら， また同宣教師から犯行声明． いろいろと多忙な御仁だ．
• 昼からも間抜け書類作成． ふと最新の latex2html で HTML にコンヴァートしてみたらどうなるだろう …… と寄り道してしまって …… VineSeed から latex2html_2000.1-0.1vl2_noarch.rpm なんぞを apt-get install してしまった． ぢつは， こいつは tetex のパス設定を破壊するのである． VineSeed おそるべし． けっきょく tetex 一式は再インストール (この拙劣なる応急対応も間抜けだが)．
• これは latex2html_2000.1-0.1vl1_noarch.rpm では問題ないみたいね．
• 贋作したのは自己評価とかいう書類なんだけど， どっちかというと内容よりレイアウトに苦闘しつつ， 書き上げて PDF 化する． 今日はじめて知った LaTeX ワザ．

	        \begin{minipage}[c][43mm][t]{150mm}
			(……内容……)
		\end{minipage}

これで高さ 43mm 幅 150mm の minipage 環境 ができる． 文書モードの tabular だの 数式モードの array だの に埋め込むには便利である． まぬけ罫線対策．
• あとは悪の組織・地球ふろんちあの中ボスである 生態系領域長の東大の安岡さんに回せばほぼ作業はおわりなんだけど， その前に小ボスであるアーキテクチャグループ長の 北大の甲山さんに検分していただくか， ということでそちらに連絡． どこが「自己評価」なんだか． 1600 作業終了． 最初の〆切は金曜日なんだけどイヤな雑用はさっさと終らせるに限る． ひとやすみ．
• ようやくこまごまとしたことを終えたので本業に復帰． しかし， どうして私は並列処理がうまくできんのだろうか． なにかコツがいるのか? 「ぎょーむ日誌」に途中経過を core dump すれば 複数ぢょぶを切り換えながら処理できるのだろうか．
• さて「ニセ」アカマツだが …… 前回, 北大に行ったときに現モデルの致命的な弱点がわかったので 計算のやりかたを少し変えねばならない． たった 1 個体の「ハズレ」のために ……
• 少しまよったけど， 現在の redpine2 は凍結保存する． さらに redpine2→redpine3 と改名． 環境変数まわりを整備する．
• まだ中身に手をつけてないけど撤退するか． 1810 研究所発． 普通に乗って 1820 京急杉田発． 快特に乗り換えて 1830 京急上大岡発． 普通に乗り換えて 1850 京急川崎発． 1920 帰宅．
• なんか近ごろ自宅における作業効率わるいなぁ …… と思ってたんだけど原因わかった． 室内が寒いのである． ということで小型セラミックファンヒーターの運転開始． おお， 文明的な生活だ．
• しかしここ数日にわたって睡眠周期が乱れているので眠い．
• さっさと寝よう …… しかしどうせまた真夜中とかに目がさめて寝れなくなるんだよね．
• 今日の食卓
  • 朝 (0715): 米 0.6 合． レタス・モヤシ・タマネギ・ニンニク茎の炒めもの．
  • 昼 (1220): 弁当． 米 0.7 合． 朝と同じ．
  • 晩 (1940): 米 0.7 合． 朝の残りで雑炊．

2001 年 11 月 14 日 (水)

• 0600 起床． 夜中に目がさめて眠れぬまま朝になってしまったのである． そろそろこの悪しき周期がリセットされるか? しかし今日一日は眠いままで終りそう． コーヒー．
• 朝飯・弁当の準備． 朝飯． シャワー．
• 0700 自宅発． 今朝も晴れ． 特急に乗って 0710 京急平和島発． 普通に乗り換えて 0730 京急上大岡発． 0745 研究所着．
• 朝から頭ばててるのに， なぜか面倒な三次元アルゴリズムの 改良案が脳裏から去らない． 関連ありそうなアルゴリズム解説ペイジをネット上で調べてみたり．
• vim6 のかしこい autoindent/smartindent 機能は便利なんだが HTML ファイル書くときは邪魔なので …… 野蛮にも /usr/share/vim/vim60/indent/html.vim のファイル名を _html.vim などと変えて読み込まないようにしてしまった． $HOME/.vimrc で :set noautoindent などとやってもうまくいかないので．
• まずいな． びみょーに風邪気味だ． 用心しなくては．
• さっさと「ニセ」アカマツ改良すべきなんだが， ちょっとだけ Abies PipeTree を手直し． こいつはバカみたいに計算時間を食いつぶすことがあるんで， 独房むかいのサーヴァー室 Dell RedHat マシンで ゆっくりと培養実験というか飼育してみようかなぁ， と目論んでいるところ． まぁ， のんびりと．
• どうでもいいが RedHat7.x 系で rpm package をまとめて 管理してるユーティリティの名前は gnorpm． めったに使わんのでよく忘れる．
• Abies はのめりこむとキリがないのでやめて， 「ニセ」アカマツにとりくんでみようとする． さて．
• 昼飯食いながら vim6 情報あさってると …… どうもユーザー定義インデントは :set indentexpr=... とかで設定するみたいね． なんかややこしそう． やはり/usr/share/vim/vim60/indent/ の中をいじるのは蛮行だったな． あるいは $HOME/.vim/ftplugin/html.vim とか．
• とりあえず $HOME/.vimrc (これは /usr/share/vim/vim60/virc_example.vim を写したもの) の「それらしそうなところ」に一行だけ

	autocmd FileType html   set indentexpr&

と追加してみた． これで HTML 用のインデント設定がちゃらになる． あるいは強引に

	autocmd FileType html   set indentkeys=o,O

などとすると， indent が作動する状況を限定できる． しかしこれよりは set indentexpr& としたほうがよさそう．
• 前にも「ぎょーむ日誌」に書いたけど vim のhighlight 設定の変更は やはり $HOME/.vimrc にて

	highlight Search ctermfg=Magenta ctermbg=Black
	highlight Comment term=bold ctermfg=Green

てなかんぢで． あと vim6 だと :set colorscheme ... で blue だの murphy だの evening だの色スキーム指定できる． これは使ってないけど．
• 「ニセ」アカマツの変更点の整理を再検討．
• おっと， お客さんからのメイル． はい， ホテル久保， 一名様のご予約をうけたまわりました． 当ホテルがお客様に誇れる三つの特長は:
  1. 羽田空港まで 25 分， 品川まで 15 分の絶好のアクセス．
  2. 隣接して設置しております 東京銭湯「第二松の湯」(第一は未だ発見されず) に特別ご優待．
  3. うーん? …… 久保ふう朝食?
  ああ， 掃除せんといかんなぁ． いやいや宿泊客が来てくださるからこそ， 当家もたまには掃除される， そう考えるべきだろう．
• 今日の午後は静かだ． しかし風邪疑惑はますます濃厚． ちゃんと対応したほうがよさそうだ． とりあえずうがいないど． 喉の痛みがマシになる．
• 「ニセ」アカマツデータ， 光合成に関するものを新しく送っていただいので， これもまとめて Perl モジュールで扱えるようにしてみようとする． あと作りかけだった Table2 ほぼ清書版を出してみたり． 茎重に関してはどれも分散でかい． たしか対数正規分布に近かったよーな気がする． そうか， これに葉っぱ重もつけないといけないのか．
• 独房内闇ルーター兼自分専用メイルサーヴァー兼数値計算機 BookPC は Yahoo! オークション 6 万円落札品とは思えぬほど よく働いてくれてる． しかしながらその闇側 ethernet カードは 1000 円ぐらいで買った安物であり， これが月に一度ほど動作不良してカーネルごと死んでしまう (正確にはカーネルは生きてるんだけど ユーザーとの連絡が絶たれる)． で， 今日はその不運な日でした． 復旧に数分．
• どうもはかどらないなぁと思いつつ撤退． 1800 研究所発． 普通に乗って 1810 京急杉田発． 快特に乗り換えて 1820 京急上大岡発． 普通に乗り換えて 1840 京急川崎発． 1900 帰宅．
• うがいして晩飯食って暖かくして早々に寝る． しかし夜中に目がさめたので風邪薬など服用してみる． えーい， 宿泊客迎撃準備ができない．
• 今日の食卓
  • 朝 (0610): 米 0.6 合． コマツナおひたし． エノキダケ卵とじ．
  • 昼 (1220): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1920): 米 0.4 合． ネギ・ニラの味噌雑炊．

2001 年 11 月 15 日 (木)

• 0700 起床． 風邪のせいか風邪薬のせいかアタマがぼんやりしてる． おそらく徐々に回復中．
• 朝飯・弁当の準備． 朝飯．
• 0755 自宅発． 晴れ． 特急に乗って 0801 京急平和島発． 普通に乗り換えて 0830 京急上大岡発． 0845 研究所着．
• アタマがぼーっとする． なんか寒いな． セーター着てくるんだった．
• かくのごとく virus に呪われぎみ， というかおそらく私の免疫系の奮闘ならびに 外来の anti-virus の副作用のせいで機能不全ぎみ． 呪われと言えば， 先天性免疫不全ぎみの M$ 系サーヴァーたちから hosho への感染アクセスの間抜けな試みがときどきあり， あの懐かしの Nimda だの CodeRed だのにやられた連中が いまだにネット上をのそのそ徘徊してるのがわかる． たぶん飼い主からも見捨てられているんだろう．
• sim_redpine.pl その他の部品である otc_data.pm Perl モジュールを直そうとしてるんだけど なかなかすすまない． Perl のおぶぢぇくと指向的な用法はたしかにケガれているけど， モジュール (ま，クラスとほぼ同じです) の切り分けが適切であれば温泉旅館的な拡張はすごく容易だ． しかしぼけ気味なのでぽつぽつとしかコードが増えないのである．
• 甲山さんから論文原稿の PDF ファイル． なぜか私も共著者の一人になっているんだけど， ダウンロードして拝見しても何を計算したのか思い出せない． 直線回帰とか F 値とか割算とか私があまり好きではなさそうな解析ばかり． たしかに Serimbu.tar.gz という凍結プロジェクトが 手元に保管されてるから， 何か計算やったんだろうけど …… ひょっとしてスクリーニングだけ， とか．
• 拡張はかどらぬまま昼飯． 意外なことに弁当食べると意識が少しはっきりした． 血液は消化器に取られてるはずなんだけどな．
• DeltaGraph の PostScript 出力→PDF は他機種では読めません， と注進したら， DG 依存症の甲山さんが DG フォント埋め込み出力とやらの 怪しげなるワザを使って原稿ファイルを作り直された． なるほど Linux 用 AcrobatReader 4.05 で見てもバケてない． xpdf ではなぜかダメなままなんだが．
• 「ニセ」アカマツ改良の続き． えーい， 簡単な機能追加に何をのそのそしてるのやら．
• さてさて …… ようやく葉っぱの重さとか最大光合成速度の推定値， とかひっぱれるようになったんだが …… この OTC 実験を使って個体を初期化する tree_initialization.pm という Perl モジュールにわりとトンデモな間違いを見つけてしまった．
• 材の密度を計算するために， 個体の茎の総重量を計算するところが， 何を考えたのかこれに針葉の重量 (かなり重い) を加えて さらに 2 倍 …… なーんてコトをやってたんだよね．
• いまごろになって推定された数値がかなりヘンだと気づいたんだけど， たとえば材密度 2.0-4.0 × 10^-3 (g/mm³) って立方センチに直すと 推定材密度 2.0-4.0 (g/cm³) となるから …… アカマツの材の比重は水の数倍， などという恐ろしい憶測結果をならべていたわけだ．
• ということで， 複雑なデータ構造をもつ観測データベイスからデータを 引き出すときには慎重に． 計算結果が常識から逸脱してないかの点検を怠らない． とはいえ科学者の科学常識なんて怪しいものだ， といま示されたばかりであるが．
• 直してざっと再計算してみる． ここらは間違った場合・モデルを変更した場合に備えて 自動化されてるからラクなのである． 以前と同じように奇麗に整形された新 Table1 印刷出力を見ると …… おお， こんどはちゃんと水に浮く木になってる． しかし， かんぢんのtree.production_rate という憶測値のパターンは変わらんな． 処理区に一本「はずれ」個体がいるとかね．
• とりあえず小林さんに連絡メイル． ここまで派手な間違いはさすがの私にもめずらしいな． 実害はほとんどないけど．
• さーて， 今後はどうしようかね …… と考えがまとまらぬままに撤退準備． 今日こそは宿泊客迎撃準備をしなくては． 1750 研究所発． 普通に乗って 1805 京急杉田発． 快特に乗り換えて 1810 京急上大岡発． 普通に乗り換えて 1830 京急川崎発． 1900 帰宅．
• 晩飯食ったあとに一休み． うう， 体うごかん． 2100 すぎに掃除開始． モノが少ないので思ったよりはラクだ． しかしホントにたまには掃除しなくては． 体に悪そうな環境を構築してしまった． でれでれ作業して 2300 ひととーりなんとかなる． 2400 まで気合で開店してるスーパーマルエツで買いもの． 宿泊客にさっさと寝てもらうためのアルコール飲料を物色． いつものことながら自分で飲まない酒を選ぶのは奇妙な感覚だ． 帰宅して料理の準備など．
• さて， こっちはなんとかなった． アタマまだふらふらするなぁ． 明日， 都立大のセミナーにお邪魔でもしようかと思ってたんだけど， どうしよう． うーむ， このように体力が低下した状況で多摩川を溯上するのは …… ひどく危険なことのように思える． なぜなのかはよくわからないが．
• 今日の食卓
  • 朝 (0710): 米 0.3 合． 昨晩の雑炊の残り．
  • 昼 (1220): 弁当． 米 0.8 合． コマツナおひたし．
  • 晩 (1940): 米 0.6 合． 煮干・コンブ・ダイコンの葉っぱの雑炊．

2001 年 11 月 16 日 (金)

• 0720 起床． 風邪からなかなか回復しない．
• 朝飯・弁当の準備． 朝飯．
• 0800 自宅発． 晴れ． 特急に乗って 0813 京急平和島発． 普通に乗り換えて 0840 京急上大岡発． 0855 研究所着．
• ここまで来るのにけっこうばてた． しかも独房群のネットワーク転送速度， いつもにもましてのろい． 回線腐ってるのか?
• 昨日は一日中寒かったので， 今日は厚いセーターなど着てみる． それでもやや寒い …… などと書いていたら甲山さんからメイル． うーむ， 今日の脱柵・遠出をススめるものだ． 体力温存してでかけるしかないか．
• とは言ってみたものの， 仕事進まん． これはいつもの怠業ではなく， ホントにばててるのである．
• だめだ． 午前中は「ニセ」アカマツコードを見直しただけで終った．
• 何やら京大農学部でインストール作業をやっておられる牧さんと 昨日から dual boot メイルのやりとり続いてる． この学部はおそろしいところで， Ruby v.s. Perl の布教合戦だのがあったり， 超然たる Kylix 使いがいたりする …… らしい． なんかすごい． 私は dual boot なんて面倒なコトはやめて Window$ なヒトたちには (ちかごろずいぶんと進化した) cygwin とか使ってもらってはどうでしょうねぇ， などと水をさしてるんだが． 間抜けデータシート作って喜ぶ ゑくせりあんが減るなら何でもいーや， という心境． 近ごろは．
• まいったな． 多摩川遡行するなら 1415 までに ここを脱柵せんといかんのか． どうしよう．
• いままで「ニセ」アカマツにおいて無情にも無視されてきた 「葉っぱの重さ」 …… これが問題解決のカギとなるか． package RedpineTree として定義されるクラス (らしきもの) に OTC 葉っぱデータを読ませるようにする．
• ああ， もう時間か． たぶん今日は独房内で服役していても 作業はそれほどススまぬだろう． そして宿泊客捕捉にもそれこそアタマ痛くなる 困難生じうる公算大． ならばこちらから出むくまでだ． 多摩川遡行を決断．
• 1410 独房発． なんか歩くとくらくらするような． 一番内側の関門は無人の電子錠式なので難なく突破． 問題は独房群正面出入口カタめる鉄壁の …… 壁際に隠れて状況をうかがうと …… ああ， いた． 受付をよそおった看守のおねーさんが脱柵者に目を光らせている． どうしよう．
• 以前から脱柵ルートのひとつとして検討していた 某搬入口に迂回してみるか． 施錠されてませんように …… よし． ここは無意味に警戒のカタい独房群のセキュリティホールみたいなもんだな． 独房群腐れファイアーウォールが telnet 受けつけるのと同じで．
• 1425 JR 新杉田発． 横浜線 JR 橋本まで 780 円． 1447 JR 東神奈川で乗り換えて 1538 JR 橋本着． 京王線相模原線に乗り換えて 1542 京王南大沢着． ここまで 140 円． 1550 都立大植物生態学教室着．
• 今日のセミナーの特別講演者・北大かとー大先生は 1600 ちょうどに着． 午前中に甲山さんに足止めをくらった遅れを回復すべく， 携帯電話で路線検索サイトを見ながら羽田からここまで到達したとのこと． 「これが最短でした」という経路はかなり複雑怪奇なもので， おそらく大先生は路線検索サイトを 親指一本でクラッキングすることによって その尋常ならざる最適解を算出せしめたにちがいあるまい．
• 1615 セミナー 開始． ハクウンボクの開花要因， 開花への資源投資が個体群動態に与える影響， オオモミジの雌雄開花フェノロジー． (かとー先生にとっては当たり前のことなのでまったく 強調されなかったのであるが) 取り扱いの面倒このうえない個体・個花序の経年データを Perl で手際よくデータベイスした上で R であれこれと解析されている． 「割り算」を過剰に多用する北大式統計学なんで (これを使わないと学位審査会に派遣されてくる 総長勅任の伝統保全監察官がクラーク先生にお詫びしろ と激怒するのである)， そちら方面はちょっとアブないかんぢ (ふつーの統計学の見地では)．
• 北大式統計学の難解さゆえに何度か昏睡しつつも， 途中でマヌケな質問などしてみる． 最後のオオモミジ話は 「メスによるオス資源の獲得→『ニッチ』分割」 というアイデアでうまくモデリングできるな， と確認． これは私の発案ではなく， ふろんちあ小ボスとして本日の多摩川遡行を命じるメイル出した 甲山さんが考えたものである．
• セミナー後， 京王南大沢横のパセオ内 天龍門 で皆さんとお食事． 研究室出る前に私がだらしなくふらふらしてたので， 鈴木準一郎さんが漢方風邪薬 葛根湯 をめぐんでくださった．
• 2120 ごろ解散． 今晩のホテル久保の宿泊客たるかとー大先生を 多摩川下流の当家までご案内する． 記録ないんでわからないんだけど 2130 ごろ京王南大沢発． 京王稲田堤で下車して JR 南武線の駅をめざす． 2200 ごろ JR 稲田堤発． 2230 JR 川崎着． 2250 京急川崎発． 2300 京急大森町着． 2310 帰宅．
• 風呂はけっこうとのことなので， 今日は第二松の湯ご案内シークエンスはカット． 買っておいたワインなど飲んでもらいつつ， 当家の ADSL 闇ネットなどらちもなく自慢してみたり．
• いざ高脚式ベッドで寝てもらおうとすると シーツがなかなか見つからなかったり， といろいろ不手際． 私は畳部屋でエアマット上の寝袋にもぐりこむ． これはこれで快適である． とはいえばてた．
• 今日の食卓
  • 朝 (0720): 米 0.3 合． 昨晩の残りの雑炊．
  • 昼 (1220): 弁当． 米 0.8 合． 鶏ひき肉・タマネギ・エノキダケの卵焼． なにぜいたくなもの食ってやがるんだよう， と思われるかもしれなけど， これは宿泊客迎撃演習の一環なのである．
  • 晩 (1830): 「天龍門」にて台湾料理いろいろ． ピーナツをカラごと蒸籠で蒸したのがおもしろかった．

2001 年 11 月 17 日 (土)

• 1020 起床． 寝る前に飲んだ風邪薬のせいでアタマぼーっとしてるけど， 寒気やふらふら感は減じたかんぢだ．
• 朝食の準備をして， ホテル久保の宿泊客 かとー先生 と一緒にたべる． 内容はたいしたものではないのだが， 今回は朝食準備に関しては成功と評価してよいのかもしれない． これまでの宿泊客の皆様のとうとい犠牲があればこそ， というべきだろう．
• コーヒー飲みつつ雑談などして， 1220 自宅発． かとー先生を京急平和島まで送る．
• 午後は風邪療養と称してでれでれする． うん? これではいつもの週末とかわらんか．
• なんか北大全体とネット接続の具合よろしくない． また DDoS 攻撃でもくらってるのか?
• 夜になって体力が少々回復してきたような気分になったので， 「ぎょーむ日誌」など書いてみる．
• 今日の食卓
  • 朝 (1040): 米 0.6 合． ダイコン・ニンジン・ネギの味噌汁． 煮干だし． 鶏ひき肉・タマネギ・エノキダケの卵焼． コマツナおひたし． キムチ．
  • 昼 (1530): 米 0.4 合． 味噌汁ののこりで雑炊．
  • 晩 (2030): 米 0.4 合． 昼と同じ．

2001 年 11 月 18 日 (日)

• 風邪療養と称して夜まで寝たきり． おかげでだいぶよくなった． 問題は週末時差だな．
• 1930 ごろにのそのそと起き出して 大森町の商店街に買いもの． 野菜など．
• あいかわらず北大ネットは反応おそい．
• 今日の食卓
  • 朝 (1130): クラッカー．
  • 昼 (1430): 雑炊． コマツナ． 醤油味．
  • 晩 (2100): タマネギ・シイタケ・ニンジン． 焼きそば．

2001 年 11 月 19 日 (月)

• 0700 起床． 風邪状態からまだ完全には回復してないけど， だいぶましになった．
• 朝飯・弁当の準備． 朝飯． コーヒー．
• 0750 自宅発． 晴れ． 特急に乗って 0801 京急平和島発． 普通に乗り換えて 0830 京急上大岡発． 0845 研究所着．
• 風邪りはびりと称して， 計算機設定など． 近ごろよく考えずに apt-get install してしまったために生じた問題点をいくつか改善してみる．
• うっかり qpopper まで RPM (つまりバイナリファイル) で更新してしまう． これは自分専用メイルサーヴァー BookPC で動かしてる qmail に対応するように ソースコードから再構築しなければならない． といっても， qpopper.spec で configure 呼び出すところで --enable-home-dir-mail=Mailbox を追加するだけ．
• X 上での BackSpace キーの挙動がおかしい． stty -a すると， ちゃんと erase = ^H とでるし …… よくわからんので ネット上の文献など調べてみて， あらたに $HOME/.Xmodmap ファイルを作り， その中で keycode 22 = BackSpace と明示的に定義してから， X を再起動すると とうまく動作するようになった．
• 北大ネットワークの応答速度， 多少マシになったか?
• どうでもよい変換フィルター書いてみる． いままで私は ( ) や { } の前後をこういうふうに空けて いたんだけど ……

	sub SetRankNumber {
        	my ( $self, $number ) = @_;
	        if ( defined $number ) {
	                $self->{ "rank_number" } = $number;
		}
	}

…… 残念ながらこれはあまりコード書きなヒトたちから 好まれない作法らしいんで改めることにした． Perl で nospace というのを作って， これに通すとそれほど嫌われない様式になる．

	sub SetRankNumber {
        	my ($self, $number) = @_;
	        if (defined $number) {
	                $self->{"rank_number"} = $number;
	        }
	}

たしかに空けなくてもそんなに見苦しくはない …… かな．
• などと愚にもつかぬことやって午前中終了．
• 弁当食ってると浦口さんから mlfitting (汎用性それなりのつもりの最尤推定プログラム) に関する質問メイル． これまでと少し違うパターンの問題なのでやや苦闘されてる模様． 弁当食いつつ， 当方でもその浦口さん製 plugin_function.so (動的に読み込まれるモジュール) 動かしてみたり．
• 結局， データみながら plugin_function.cc で定義される 確率分布とその期待値の関数型， parameter.txt で設定される パラメーター初期値をことごとく変えるとうまくいった． うーむ， 「汎用」への道は遠い． なんかノンパラメトリック回帰だの local likelihood/regression だの やってたほうがマシじゃないかと思えてきましたよ． そっちはそっちで鬼哭啾啾という気もするが．
• なにげなくそちら方面の恐ろしげな手法について調べていると， 今度は本ものデータつきメイル． うーん， カエデの直径成長か． これはこれで難しそうだ． まじめにやろうとすると．
• よくわからないので， 時間稼ぎメイル送信． その間にちょっと検討してみるか． えーい， このへんの問題ほっておくんじゃなかった (小川モデリング以前からいい加減にしてた部分である) …… と思ったらすぐにまた来信． 時間稼ぎになっていない． うう．
• 手元の確率分布本をひっくり返してながめてもよい思案うかばず． 腹も減ったので 「上限下限とるなら …… いっそベータ分布はどうでしょうかね?」 とこれまた時間稼ぎメイルを書いてから， すばやく逃走準備． 1820 研究所発． 普通に乗って 1830 京急杉田発． 特急に乗り換えて 1840 京急上大岡発． 1910 帰宅．
• 今日は御役御免というメイルをいただきばたっと寝てしまう．
• 真夜中に目をさまして， ふとデータで遊び始めたのが運のつき， 明け方まで並べかえたり均等分割したり作図してみたり．
• 得られた良い結果としては， やはり世の中には正規分布などそうそうありはしない， ということであろうか． 下の図は最も小さい 百十数個体と同数の最も大きい個体たち． しかしこれって何分布だろうか． 均等分割した標本の傾向としては， 平均と標準誤差は比例関係に近いようだ． 個々の分布はベータ分布とかにあてはまるだろうけど， いまいち釈然としない．

• いかん 0430 すぎてる．
• 今日の食卓
  • 朝 (0710): 米 0.6 合． ホウレンソウ・タマネギ・エリンギの炒めもの．
  • 昼 (1220): 弁当． 米 0.8 合． 朝と同じ．
  • 晩 (1930): 米 0.6 合． 朝の残りで雑炊．

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 合． ダイコン・ニンジン・シイタケの味噌汁．