ぎょーむ日誌 2003-07-(11-20)

2003 年 07 月 11 日 (金)

• 0750 起床． 朝飯． コーヒー． 0840 自宅発． 曇． 0850 研究室着．
• 起動状態よろしくない． 自分の．
• ここしばらくトラップされてる， 観測値 = 乱数 + 乱数の推定問題 …… 右辺がどちらも正規分布だと， ひじょーにわけわからんのではというやつね． In All Likelihood - Statistical Modelling and Inference Using Likelihood にそのあたりの解法があるわけだが (お， 著者 のペイジだ)， 部分的にモーメント法を使ったりしていてなかなか面倒だ．
• というふうに午前が終ったのだった． 北大構内走．
• もどってきたらいきなり会議． 来年 8 月の生態学会釧路大会関係． とりあえずのメモ．
• 今日も仕事がススんでない． 1900 研究室発． 北大図書館に立ち寄る． 2020 帰宅． 体重 71.6kg． 晩飯．
• [今日の運動]
  • 北大構内走 1240-1315． ストレッチング．
• [今日の食卓]
  • 朝 (0810): 米麦 0.7 合． タマネギ・卵の炒めもの．
  • 昼 (1630): 弁当． 研究室お茶部屋． 米麦 0.7 合． タマネギ・卵の炒めもの．
  • 晩 (2200): 米麦 0.7 合． レトルトパウチドものの麻婆春雨．

2003 年 07 月 12 日 (土)

• 1010 起床． 朝飯． コーヒー．
• 昨晩から図書館で借りた 旧石器遺跡捏造 を読んでる． 驚くほどちゃちな手口のワンパターンな 「ねつ造」が繰りかえされていたと報告されてるわけだが …… 考古学会の研究者たちも東日本各県も こういうのに 30 年間ちかくもだまされてたわけね．
• 以前に， われわれ職業的うそつきであるモデル屋の職業倫理は ただひたすらにばれにくくてだましやすい うそを創造することにつきる， と述べた． で， この ゴッドハンド こと 藤村氏 はそれをいちおー成し遂げたわけではある …… もし生態学周辺の分野で 30 年間も通用するナニごとかを 見つけることができたなら， それはかなりスゴいことだよねえ．
• ねつ造は科学の他分野でもひんぱんに行われている (最近ではたとえば J.H. Schoen 事件 とか MOX 事件とか)． とはいえ， 多くの科学者は (少なくともあからさまな) ねつ造発見のたぐいには手をそめようとはしない …… いったい， ねつ造という作法のどのあたりがよろしくないのか? このへんのまぢなハナシなんぞは 「自分は誠実だ」 と確信して疑わぬまぬけな連中にでも まる投げでもしてやることにして， 一介のプロを自称してるうそつきならばどう評価すればいいのか?
• ねつ造によるナニかはたしかにわかりやすくて (何しろ非常に好都合なものが「見つかって」しまうんだから)， 分野と時代によってはなかなかばれない状態を維持できるのかもしれない． しかし， うそは有限時間内に確率 1 でばれるものであり (永久にばれないうそは誰からも無視されたうそだ)， いったんばれてしまったら ねつ造の上に組み立てたハナシは実にもろい． 決壊した堤防のごとく， あるいはドミノ倒しのごとく虚構体系がくずれていく ……
• ということで， 旧石器遺跡ねつ造事件から学ぶべき教訓は， われわれは「奥行きのある」うその構築をめざさなくてはならぬ， ということにちがいない． いくつかのうそ部品が破綻したとしてもシステム全体は そう簡単にはぐらつかない冗長性と頑健さを備えている …… これこそが「総体としてばれにくい」うそなのだろう． こーいうのを設計する方法論， それは私にはまだわからぬ部分がある． しかし「一点突破で壊れるうそ体系」は「すぐばれるうそ」 と同じく素人じみた手抜きとして忌避されるべきものであり， プロの仕事ではない．
• 1420 自宅発北大構内走． つごうよく曇天． 1510 帰宅 体重 71.2kg． 昼飯．
• 1700 自宅発． 曇． 1710 研究室着．
• しかし研究すすまず． まぬけだ． 2150 研究室発． 2200 帰宅． 晩飯．
• [今日の運動]
  • 北大構内走 1420-1510． ストレッチング．
• [今日の食卓]
  • 朝 (1030): バナナ． トマト．
  • 昼 (1520): 蕎麦． 納豆．
  • 晩 (2300): スパゲッティ． タマネギ・マイタケのカレー．

2003 年 07 月 13 日 (日)

• 0910 起床． 朝飯． コーヒー．
• しかしこの休日だらだーらモードはどうにかならんのかねえ．
• 1210 自宅発北大構内走． 今日もよい具合に曇天． ちかごろどうも直射日光にヨワくなってるような気がする． なんというかバテ感がつのるというか． 1300 帰宅． ひざとか冷却． シャワー． 体重 72.0kg． 昨日は水分とりすぎたか? 昼飯．
• 1430 自宅発． すこし晴． 北大図書館 に立ち寄る． 日曜日も 0900-1700 開館している (そして北大関係者でなくても登録証があれば利用できる)． なかなか素晴しい． 大学におけるさまざまな勉強支援ぎょーむの中にあって， 図書館が他のあれこれ (授業とか計算機網とか) と比べてなかなかよろしーのは， 運用とかに関して 他大学との情報交換がそれなりに活発だからではなかろうか． 1530 研究室着．
• 休日ぼけの作業療法ということで， 軽めの労働として下請け原稿ファイル整備を少々． Table まわりがだいたい片づいた． 時刻は 1750． 一休みしますかね．
• 植物学会大会 のシンポジウムの要旨がなかなか書けん． どういうヒトが聞きにくるのかよくわからん， というのと， 今回の眼目のひとつになりそうな計算 (地下部での水のとりあい) をまだまったく計算してない， というのと． 後者はまったく私自身の責任だが．
• しかしいつまでもうだうだしていてもしょうがないんで， 書いてしまう． で， さらに英語版の提出も必要とされてんだよねえ．
• 夜になるとどんどんやる気が失せる．
• 近ごろでは敵国語作文のときに検索エンジン google が手放せなくなった …… 「葉っぱなんかがかせいだ資源をその場で使うこと」 を何と言えばいいか， と考えていてアタマの中でふと ``in situ resource utilization'' なる怪しげな連語を思いついた． こんなヘンてこな造語を使ってる莫伽はどれぐらいいるんだろうか， と検索してみると …… 1160 件 hit． つまりすでにあちこちで使われているコトバだったようで． いろいろな使用例が見つかるわけだけど， これなどはたいへん感動的な in situ resource utilization の一例．
• 2230 ごろによーやくできた． 締切は明日なんだけどネット経由で提出できるんで， 明日に提出しましょうかね． 日曜日の夜も苦闘する大学院生たちがいる …… 植物学会の要旨がまとまらなかったり， 呪われ C-N coder から解放されなかったり， データ隠匿の方法論をひたすら工夫していたり． そういうヒトたちに 「いやー，本日の仕事は片づいてしまいましたよ」 と自慢してくやしい気持ちにさせておいて 2320 研究室発． 2340 帰宅．
• [今日の運動]
  • 北大構内走 1210-1300． ストレッチング．
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (1000): 米麦 0.7 合． タマネギ・マイタケのカレー．
  • 昼 (1340): 蕎麦．
  • 晩 (2410): パン． タマネギ・マイタケのカレー．

2003 年 07 月 14 日 (月)

• 0730 起床． 朝飯． コーヒー． 0830 自宅発． 晴． 0840 研究室着．
• 植物学会大会 のシンポジウムの要旨を投稿しようとしたところ， 大会委員会側から警告のあった 「投稿したはずの内容が途中で切れる」 問題に遭遇する． しばらく挙動を観察して原因がわかった …… これは英文だろうが和文だろうが， 要旨文中に &があらわれると， そこから後ろは全部すてられてしまう， という要旨投稿用 CGI プログラムのバグだ． & を and と書き換えて再投稿． これで問題ナシ．
• 7F に行って雪野さんと科研費の使いかたについていろいろ教えていただく．
• 午前はこういったことどもとメイルのやりとりと 粕谷さんから送っていただいた計算メモの検討とりかかりで 終わっちまった …… うーむ．
• で， いつの間にか講座セミナー． 本日は D2 小菅せんせー． 本講座の草本院生たちはつくづく MineSweeper 苦手型のよーで …… 林床の多年生草本の 光合成を測って繁殖への資源配分を調べます，と地雷を踏む → それやるなら生産生態学の方法にのっとってきちんとやれよな (地下部の重量とかまで調べるつもり?) といつもの甲山さん bomb (といっても立腹してるわけではないけど) → ゲイム終了 (小菅君の場合は投影画面を見て苦悩する)． そろそろ このへんの定石ぐらいは覚えてもいいんではないのかのう．
• これまたよくあるパターンで工藤さんから 「それってそういうハナシではなかったんでは」 とモノ言いというか救命浮輪の投入というか …… で， その浮輪にすがるべきか否かと， またハムレットのごとく激しく懊悩してみせる小菅せんせーはともかく， 聴衆たちはなるほどそういうわくぐみだったのか， と方向性が少しわかる． セミナー室のうしろのほうで 「隠匿」 「今回は意図のほうの」 といったコトバがささやかれる．
• とかやってるうちにうやむやのまま終劇． この新観測計画 (なンと D3 の春から着手，ですよ) のとーとつ感の強さを不思議に思った観客の中には， これまで継続してきた他の観測ぎょーむで 何か素晴しい成果が得られたのではないか， この新計画はそもそもそちらから注意をそらすための …… といった偽装計画説もウワサされた．
• 甲山さんが屋久島の直径成長データ． なるほど variance が あがったり下がったりするのは面白いかもしれない． 例によって迂遠というか間接的なアプローチなんで， こういうパターンが時間相関によるものかどうか， 実際のところよくわからない．
• 北大生協クラーク館に行ってアルバイト求人広告を眺める …… ソフトウェア開発だと時給 1000 円の自由時間労働． いや私が職さがしやってるわけでなくて， 科研費で学部学生さんを傭兵化しようかな， ということで札幌の雇用相場を調べている． 地域生態の院生たちは …… 私の下請けやってるヒマあるんなら自分の仕事ススめろというかんぢ． というか私が彼ら彼女らの下請けをやってるわけで． あ， かとー先生に働いていただく， という一手はあるな．
• なぜ全学を視野にいれてヒトさがしをやってみたいと考えるか …… 何もない北辺の大学ではあるけれど， 学部学生 10000 人・大学院生 6000 人の中には 私などより (たとえば) プログラミングとかに優れた学生はいくらでもいる． それが理由であり， この点はまったく疑いない． 問題はそういう才能を探しだして雇えるかどーか， ということで． 探すのも雇うのも難しそうだ． どうやったもんかなぁ．
• 雪野さんに教えてもらった学生雇用条件あれこれ (地環研事務室の方針)．
  • 可能なかぎり地環研所属であること (なんじゃこれは)． 他の学部の学生・院生の場合は 学生証のコピーを提出 (うーむ)．
  • 時給 950 円が上限 (いやはや)．
  • データをあつかう労働は地環研建物内でのみ許可される …… 雇用側が「監督」するためだそうで (いんたーねっととか知らないんだろうな)．
  このあたりのばかげた規制は， 科研費を裏金化する 効率を悪くしよう， というつもりなんだろうなあ． もちろんぜんぜんその目的とするところを達していないわけだが (なんとしても裏金つくりたいヒトにとってはほとんど障害になってない)． こういう下らぬ「校則」どもが生成されてきた背景は …… どこかの大学で不正流用がばれて報道される → 文部科学省から「対策たてろ」とお達しがくる → 上のような制約とか下のような意味不明な料金表を作って 「本学はこのような対策たてました」と上表する， という阿呆らしいお役所仕事を繰り返した帰結だろう．

• 人買いの策も思いつかないので， とりあえずの予算使用ということで液晶ディスプレイ発注． いま使ってるのは数年前に Yahoo! オークションで買った私物なんだけど， 一日に何回かちらつき解消の調節をやらないといけないんで．
• まぬけなことに仕事はあまりススまぬまま 1920 研究室発． 1940 帰宅． 体重 72.2kg． うーむ． 晩飯．
• 電気代 1415 円． 電話代 1900 円．
• [今日の運動]
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0800): 米麦 0.7 合． タマネギ・マイタケのカレー．
  • 昼 (1320): 弁当． 小菅セミナー拝聴しつつ． 米麦 0.7 合． コマツナ．
  • 晩 (2030): 米麦 0.7 合． キャベツ・タマネギ・ショウガ・鶏レバの炒めもの． ワカメ． トマト．

2003 年 07 月 15 日 (火)

• 0710 起床． 朝飯． コーヒー． 0810 自宅発． 曇． 気温低い． 0820 研究室着．
• うーむなかなか仕事にとりかかれん． 問題の最尤推定計算をつらつらと見直してみたり．
• で， この日の午前は作表下請けぎょーむに逃げてしまった．
• 曇天下の北大構内走． 昼飯．
• 1400 から 2050 まで M2 大澤君とシウリザクラデータハンドラーの Perl モジュール作り． なにしろ元ファイルがゑくせるなので， かなり呪われきってるんだが (同じ内容を記述するのに何とーりも書きかたがあったり) …… ともかく全データをきちんと格納するところまではできた． 5ヵ月ぶんの葉っぱ・シュート増減データ (in 半規則的フォーマット) をいろいろ間違いチェックをやりながら 読み込むのに Celeron 500MHz 鈍足 ThinkPad240Z でも 1 秒． Perl はやはりこういうぎょーむは速い．
• たぶん 2002 年度データを自動解析する上での諸問題において データ処理技術側面に限定すれば難題はこれですべて解決した …… とはいえ，統計学的・生態学的側面に関してはまだ「これから」 なんだけど (そして 2003 年度データはまた異なる 様式に格納されてしまっている!)． やはりある程度以上に複雑な時間・空間・階層の 構造をもつデータが得られる場合には， 腐れゑくせるなんぞは経由しないでいきなり データベイスに格納する経路をつけるべきだな (入力端末としては web browser とか?) …… この問題に関してはすでに手遅れなんだけど．
• なんとなく大学院生たちと晩飯を食いにいくことに． 2105 研究室発． 北 12 西 3の担々麺専門店 喜和苑 で坦々麺． 2155 帰宅．
• [今日の運動]
  • 北大構内走 1240-1310． ストレッチング．
• [今日の食卓]
  • 朝 (0730): 米麦 0.7 合． キャベツ・タマネギ・ショウガ・鶏レバの炒めもの． ワカメ． トマト．
  • 昼 (1320): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・タマネギ・ショウガ・鶏レバの炒めもの．
  • 晩 (2120): 北 12 西 3 の 喜和苑 の坦々麺 (辛さ 1; 720 円)． 評判の店だけになかなかうまく， しかも普通もりですごく満腹する． もうちょい塩味がウスければなあ (しかし東日本人には最適な濃度なのだろう)．

2003 年 07 月 16 日 (水)

• 0710 起床． 朝飯． コーヒー． 0820 自宅発． 曇． 0830 研究室着．
• ぎょーむ日誌書きつつ Perl データマンジング本 などつらつらと見返す．
• さっさと仕事すべきなのに， 科研費で買った 実用 Perl プログラミング をぱらぱらと見てしまう …… かなり良い． この本を読んで， SWIG (simplified wrapper and interface generator …… swig には「がぶ飲み」という意味もある) というちょー便利そうなモノがあることを知った (いまさらながら)． 他にもデータ永続化とか有用な情報いろいろ． 1997 年原著出版ということで， やや古めかと思ってたけどそんなことはない． Perl をぎょーむで使わねばならぬヒトにとって， この実用 Perl プログラミング (Advanced Perl Programming) と Effective Perl は必携だろう． どちらも応用ワザのための基本を重視してるんで．
• さーて， いつまでもふらふらとしてるわけには …… お， /.J にも R の紹介 が! ほらほら， ゑくせる統計関数 の壊れ具合がどんどん報告されてますよ (cf. 青木さん (群馬大) の 「ひどい話です!」)．
• ということで， ゑくせるやめますか， それともデータ解析やめますか．
• それはともかく， まだ解決しない粕谷さん宿題にとりくむ． ある小集団の平均値が確率変数として与えられ， かつその小集団内の標本もまた確率変数 (平均ゼロ) である， という混合モデル的状況 (後者の確率分布の平均がゼロなので random effect 的状況 というべきなのかもしれないが)． で， とりあえず確率分布がどちらも正規分布としている． 分散はありとあらゆるところで異なる．
• こういう計算やったらいいんじゃないでしょうか， ということでたたみこみ尤度方程式を書いたメモを先週つくったんだけど …… これではうまく最尤推定できない， ということがわかった (集団内分散と集団間分散がうまく分離されん)． Searle, Casella & McCulloch (1992) Variance Components なる書物に上と似たような問題における 尤度方程式が掲載されてると 粕谷さんが調べてくれたんだけど， 私のそれとは少しばかり異なっているんだよねえ．
• で， このあたりの相違が何によるものなのかよくわからん， という状況が続いている． よくわからぬまま午前終了． 北大構内曇天走． 昼飯．
• 午後も出典 Variance Components 数式をひねってみるけどよくわからん． そこでいったん自分の定式化をゼロから再検討してみる …… と， わかってしまいましたよ． 粕谷問題の想定をうまく定式化できていなかった． この nested 状況はよく考える必要があり， 小集団の平均は確率変数だけど， その小集団にとっては定数のごときものであり， そこを中心にして小集団ごとに異なる標本セットが得られる， というわけだ．
• というふうに数式を書き直してみると …… どひゃー， 多量の正規分布密度関数が積のカタチになって登場か! 正規分布ぎらいな私にとっては， わけわからんムシどもがひしめきあってうごめいている 小部屋にほうりこまれたような気色わるさを覚える …… かとーオフィスから逃げだし北大構内逃避の旅に出る． ただし紙束はかかえて．
• 北大生協の自販機紙コップコーヒ飲みつつ， 「えーと， 正規分布の積 → 和の正規分布 に帰着できるハズだ……」 とココロ落ち着けさせつつ計算問題に取り組む． そういう式変形を行い， かつ再生性を悪用するカタチにもっていこうとすると， おのずと小集団の標本平均と標本分散が導出されてくる．
• えらく難しい問題やってるなと思われるかもしれないけれど， 実際の難度は数式頻出分野でいうと 学部学生レヴェルの演習問題程度なんだが (じつは微分積分ができる高校生でも問題ナシ …… ただし初等的な統計学の理解が必要)．
• どうやら， (途中をはしょると) 下のようにまとめられるんでは， ということがわかってきた． f(., μ, σ2) が Gauss な確率分布密度で， Mi,j と Vi,j はそれぞれ小集団ごとの標本平均・標本分散．

• これは Variance Components 本のそれに近づいている． Mi,j と Vi,j は十分統計量 (sufficient statistic) というやつだろうか， と思ったんだけどそうではなさそうだ． これらが導かれてきたのは f(., μ, σ2) が正規分布であることに強く依存している． 十分統計量なら確率分布にはよらないはず， と理解してるんだけど …… うーん， 実は f が指数関数族なら何でもいい， ということなら十分統計量と言えるのかもしれんなぁ． どうなんだろう．
• それはおいとくことにして， かとーオフィスにもどる． とりあえず粕谷さんに連絡メイル． それから， ここまでの計算を LaTeX でまとめてみる．
• どうやら粕谷さんのほうも小集団ごとの 標本平均・標本分散を使った定式化に気づかれたようで …… 私のやりかたとは少し異なるようだけど．
• カエデ下請けぎょーむで CVS リポジトリの挙動不審な点を 浦口さんと解明してみる …… ディレクトリ追加した場合は cvs checkout せんといかんの? 追跡調査が必要．
• 最尤推定問題は簡単化された尤度方程式が出る書きかけ状態で make dvi && make pdf してアップロードしておく． また連絡メイル． 今日はここまで．
• と思ったらお茶部屋に白板にはってあったパズルなどにはまってしまう． さらに M2 森君とマイクロサテライトによる父親探しの 確率論的モデルの検討． 2230 研究室発． 現金がないので買いもの行けない． 2240 帰宅． 体重 71.6kg． 晩飯． 計算あれこればて．
• [今日の運動]
  • 北大構内走 1235-1310． ストレッチング．
• [今日の食卓]
  • 朝 (0730): 米麦 0.7 合． キャベツ・タマネギ・ショウガ・鶏レバの炒めもの． トマト．
  • 昼 (1320): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・タマネギ・ショウガ・鶏レバの炒めもの．
  • 晩 (2300): 米麦 0.7 合． レトルトパウチドのカレー． トマト．

2003 年 07 月 17 日 (木)

• 0710 起床． 朝飯． コーヒー． 0810 自宅発． 今日も曇． 0820 研究室着．
• 昨日の計算のみなおし．
• かとー先生とか露崎さんとかふだんならもっと遅くに来るヒトたちが 9 時台に登校してくる …… 聞けば本日は地環研事務による何ちゃら会計監査がある， ということで雪野さんが招集かけたとのこと． 監査うんぬんはわけわからんけど， 研究室の中では雪野さんの権力が最もスゴイとはわかる．
• 昼すぎに現金調達に出る． 晴． 過去数十時間にわたって文字どーりの無一文状態が継続してるので． 帰りがけに北大生協書籍部に立ち寄る． 昼飯調達． 研究室にもどって昼飯．
• 粕谷さんから文献おしえていただく． 統計的検定は有効か　−有意性検定論争− (このリンクは福井県立大・佐野さんの要旨ペイジ) …… 古い本ではあるけれど (原著は 1970) なかなか面白そうだ． 調べてみるとここ地環研図書室にも蔵書されている． あとで借りてこよう．
• 統計メモ， 計算みなおしやって notation を整理して， スコア関数 (対数尤度方程式をパラメータで偏微分したもの) をつらつらと書いてみるだけなのに， 時間がかかる． 途中で宣教師殿とお茶部屋で 樹木モデリング雑談なんぞをやってたせいか?
• 夕方になるにつれ， かとーさんたちが朝から待機させられてる 「監査」とやらの実態がわかってくる． 実施してるのは文部科学省， らしい． 北大内のいくつかの研究室が今回の対象． 目的は不正使用の発見． この研究室が対象になってるのは 甲山さんのインドネシア調査科研費の額が大であるから， らしい． かとーさんの給料がそのへんから出てたので 事情聴取を受けねばならなかったそうで． 私が関与してないのはこのお金と無関係だから …… なんだけど， それ以外にも理由アリとのウワサもあり．
• なんかいろいろややこしいことをやってる．
• 早めに撤退． 1850 研究室発． 1900 帰宅． 体重 71.6kg．
• Yahoo! オークション (またかよ) で落札したでぢかめが到着． 以前のやつが行方不明になってから半年が経過したんで (とはいえ，まだそこらにあるような気がする)． Linux (kernel 2.4.20) でも何の問題なくデータやりとりできる (via usb-storage kernel module)． 脱力するほどの簡単さだ …… というのも， 私のカーネルはそのあたりがなぜかしら すでにキチンと準備していたもんで， カーネルまわりの設定変更ゼロ． mount -t vfat /dev/sda1 /mnt/...， まぁ， このへんの詳細はまた後日にでも．
• [今日の運動]
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0730): 米麦 0.7 合． 納豆．
  • 昼 (1310): 研究室お茶部屋． 北大生協ポテトサンドイッチ．
  • 晩 (2230): パン． キャベツ・トマト・マイタケの炒めもの．

2003 年 07 月 18 日 (金)

• 0840 起床． ねぼうした． 0850 自宅発． 晴． 0900 研究室着． 朝飯． コーヒー．
• 今日も朝からかとーさんが待機している． で， 事務からの電話に呼びつけられていった …… 数十分後に戻ってきて曰く 「けーさつ官のような狂った眼つきのヒトたちでしたよ ……いろいろねちねちと聞かれて」． もちろん， はっかーかとー先生が 中央官庁の派遣監査官ごときにヘマをやるはずもなし．
• しつこく計算みなおし …… 私の尤度方程式と， 粕谷さんが別経路で導出した尤度方程式は別経路で導出したのは， 合致してる， と考えてよさそうなんだが．
• 一週間ほど前に生物統計メイリングリスト biometry に 「正規分布とある確率分布の『近さ』を定量的に表現するには， どうしたらよいか」 という質問投稿があったので， そんなの Kullback-Leibler 情報量 (KL-divergence) つかえばいいじゃんという投稿してみた． すると KL-divergence って何をあらわしてるんでしょう， と質問がまた返ってきたので， 教科書的に

  直感的な意味合いとしては， ある確率分布 X を (この問題の場合) 正規分布の代替物として使った場合， どれぐれらいの情報量が損なわれるか， をあらわしたのが KL-divergence (distance) です．

  と回答してみた (この場合は，あえて正規分布などと書いているけれど， KL-divergence はどんな確率分布にも適用可能)， するとまた質問者の原嘉隆さん (九州沖縄農研・土壌) から「情報量ってなんでしょう」 という私信いただいた． こう， 質問されてる内容は 私がちゃんとよく理解してない領域に徐々に近づきつつあるわけだが …… このような説明を書いてみた．
• 情報量は情報理論の定義される量です． うまくいくかどうかわかりませんが， 我流の説明をしてみます．
• ある計算プログラム S が， 完全な正規分布の乱数を生成できる， とします． このときにこのプログラム S を「情報源」 と呼ぶことにします． S から出てくるのは乱数なので予測不能なでたらめさをもった情報源です． ただしその出力乱数は正規分布に従ってる， というのがわかっているので， まったくのでたらめではなく規則性もあるわけです．
• つぎに S を改造したプログラム S' があって， それは正規分布に似ているけれど 少し異なる未知の分布 X にしたがう乱数を生成している， とします． この場合， あたかも正規分布情報源である S から出てくる値に 何かわけのわからない「雑音」 を足したり引いたりして歪めてしまった値を S' は生成している， と解釈します． つまり S に比べて S' は予測不能なでたらめさが増大している， ということです…… S から出てくる乱数は正規分布にしたがっている， とわかっている， しかし S' はそれをどう歪めているかよくわからない， という状況ですね．
• KL-divergence はここで「もとの情報源である S に比べて S' はどのぐらいでたらめさが増大してしまったか」 を表しています． つまり私が前のメイルで「情報量の損失」 と言ってる部分は「でたらめさの増大」 と言い換えられます．
• つまり「でたらめさの増大」⇔「情報量の損失」 だから， 情報量ってのは「でたらめではなさ」 をあらわしてると言いたいわけ． あるいはでたらめさは予測不可能性で， 情報量は予測可能性 …… というような言い換えばかりではいかんのかな． ともあれ， Claude Shannon ってのは， よくもまぁこんな概念を考えついたもんだ (この著明な数学者が 亡くなった のはわりと最近である 2001 年 2 月末 …… なぜか「ぎょーむ日誌」でも言及している)．
• というふうに午前が終わって昼飯． 直射日光つよそうなんで北大構内走はあとまわし， の予定．
• 地環研図書室で 統計的検定は有効か　−有意性検定論争− を借りてくる． なるほど， これは原著者も訳者も少しばかり異端なかんぢなのかもしれない． とはいえ， 批判そのものは健全なものという気がする (もし仮にいくらかまとはずれであったとしても)． そもそも科学の道具としての統計学には， 根本原理から運用実態にいたるまで さまざまなアヤしさがあるのだから (数学としては首尾一貫してるにもかかわらず， あるいはそれゆえに)． まぁ 統計学に限らず， 科学というのは方法論からその帰結にいたるまで， 永劫の懐疑をつきつけるというのが本質なんだからしょうがない．
• 昨日の午後の仕事があまり進捗してなかった理由を思い出した …… 以前にやった Pasoh データ怪しげ trade-off 解析 (共犯久保) をなんとか救援しよう， と甲山さんが言ってきて …… そのやりかたというのが， ほとんどクロに近い灰色被告側弁護士の法廷戦術のごときもの， だったんで 「そりゃダメでしょう」 とか議論してて時間を費したんだよなぁ．
• しばらくは放置できそうだけど， Pasoh データ解析からは逃げられないんだろうな． となると， かかる口先のごまかしではなく， もう少し巧妙な策が必要という気がする． いわば， Hubbell の熱帯林中立共存説のごとく 誰でもいちゃもんつけられるけど， より人為的な対抗モデルしか提示できない (「本家」たる木村資生の遺伝子の中立説もそうなのか?)， というようなヨワいふり→誘出撃滅型のワナを内包したわくぐみ が必要なんじゃなかろーか? ふりも何も実際に当方はすさまじくヨワっちいわけだが．
• ナニやら某研究科事務室のおぢさんたちがぞろぞろと かとーオフィスまでやってきて，
  「かとー先生はこちらでしょうか」
  「まったくあの会計検査院の役人のヒドさときたら……」
  「昨日の◯学部の監査のときみたいになってしまったらどうしようかと」
  「よくぞあれほど手際よく…… 助かりました」
  ということで， かとーさんが日本国政府から派遣されてきた勘定目付を 短時間で退治した御礼ということでこういうの持ってきた． さすがわ，かとー先生．

[かとー先生祝捷献上品]

岩見沢の菓子屋 赤いリボン の suger maple 菓子つめあわせ (といっても岩見沢にサトウカエデがあるわけではない)． JR タワー (JR 札幌駅ビル) 地下で調達されたとおぼしきもの． で， その費用はどこから ……

• 他人に雑用を押しつけるのが巧妙なる D2 赤坂君の陰謀にかかって， 北大図書館まで文献複写に． 記名してから書庫に入場し， 不可思議な構造をもつ内部をさまよう．

[北大図書館書庫 -- 迷宮]

「これこそは北日本最大の……」 というわけでもないが， この階段のヘンてこさは迷宮的な雰囲気がある (しかも階段周辺の天井はすごく低い)． 前の階段を下ってむこうの階段を昇ると …… さて， どこに行くのでしょう．

• 北大図書館の基本ワザ． 本館書庫内の書籍・雑誌を「借りずに」 複写する方法．
  1. 地環研図書室で 「複写機使用伝票」 なる三連 (A/B/C) の紙切れをもらう
  2. 北大図書本館に行き， 書庫入室の記名をする
  3. 係員に頼んで「地環研」用の コピーカードを借りる
  4. 書庫内で必要な文献をさがす
  5. それを「書庫内の複写機」 (書庫 4F 西側に設置) でコピーする
  6. このときにカウンターの変化 「複写機使用伝票」 を記録しておく
  7. 文献をもとの位置にもどす
  8. 書庫退出時にコピーカードを返却し 「複写機使用伝票」 を提出する
• 研究室にもどって複写してきたのをぱらぱらとながめる．
• 応用統計学会 ML で「連合大会プログラムできた」 という案内があったんでダウンロードして読んでみる (私は発表申し込みとかやってない)． この「連合大会」というのは 応用統計学会・日本計量生物学会・日本統計学会 (ついでに日本行動計量学会) が「連合」して大会をやります， というもの．
• これだけあれこれくっつけても， 発表者総数は生態学会のそれよりはずいぶんと少ないような気がする …… 一人 20 分間の口頭発表， ということでハナシをじっくり聞きましょうということなのかな． で， 発表者も厳選されてる， と．
• とはいえ， やはり本邦は統計学者の数があまり多くはないんではなかろーか． ちなみに国内のどの大学にも統計学部はもちろん統計学科もない (しいていえば 統計数理研究所 の総合大学院統計科学専攻)． また， 統計学の研究室がある数学科もそれほど多くない (たとえば北大理学部 数学科 にはない)．
• 会期は今年の 9/2-9/5． どんな様子なのかいっぺんのぞいてみたいんだけど (発表しない参加者は会場までいって会費を支払えば参加できる)， 現在の仕事進捗では無理のような気がする …… 9 月末の植物学会 札幌大会 ってのがあってですねえ．
• 2120 研究室発． 2140 帰宅． 体重 71.6kg．
• [今日の運動]
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0910): 研究室お茶部屋． パン． キャベツ・トマト・マイタケの炒めもの．
  • 昼 (1300): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・トマト・マイタケの炒めもの．
  • 晩 (2230): 米麦 0.7 合． キャベツ・トマト・エノキダケのカレー．

2003 年 07 月 19 日 (土)

• 0810 起床． 休日なのになぜかこの時刻に起きてしまった． 朝飯． コーヒー． といってもだらだらしてるわけだが．
• む． この hosho.ees.hokudai.ac.jp につながらんな …… 原因は A 棟内部にある， というところまでつきとめたんだけど， しばらく放置だ． どうせ土曜日だし， ちょっといま工作モードで手がはなせない．
• さきごろ安値で落札したちょー小型でぢかめ Casio EX-S1 用のポケット内保護ジャケットを作っていて， ですね …… これ， とても軽くてポケットにつっこんでいても ほとんど気づかないほどなんだけど， 前面のレンズ (またしても固定焦点) と後面の 1.6inch 液晶パネルはちょっと何かでカヴァーしてやりたいなあ， ということで．
• もちろん既製品はいろいろ (これら以外にもサードパーティー製多数) あるわけだが …… もはや調達資金ナシ． ということで自作している． 材料は NASDA バイブル (NASDA 所属時代に配給されたロケット&衛星スペック一覧本， いわゆるバイブルサイズ) の合成皮革製のカヴァーを再利用してみる． 工作といっても採寸して裁断して折り曲げて二箇所を縫うだけなんだけど．

[でぢかめ自作上着]

材料は NASDA バイブルの合成皮革カヴァー． これを装着したままで撮影できる． ところで， この写真はここに写ってる EX-S1 じしんで撮影したもの． どうやったんでしょう? …… 答はこちら．

• 1110 自宅発． 晴． 1120 研究室着． ネットワーク不調の原因は階 hub にあった． こいつの電源ひっこぬいて入れ直してネットワーク復旧． しかし， どういう hub なんだろうね．
• うーむ， 仕事に着手できん． 学校に来ても休日だらだらモードのままだ．
• と思ってたらちょうど具合よく， 足寄 での調査を終えた九大の 田代さん がみえたので， お茶部屋でカエデと樹形雑談．
• そのまま晩飯を食いにいくことに． 北大正門前 100m のジンギスカン屋 ヤマダモンゴル へ (特に匿名希望する料理人氏も同行)． 年輪画像解析のハナシなども伺う …… Java で plugin 作れる ImageJ (NIH-Image の後継) がなかなかよろしいとのこと．
• 飯のあと解散． 2000 帰宅．
• [今日の運動]
  • 今日は運動してないなあ．
• [今日の食卓]
  • 朝 (0830): パンケイキ．
  • 昼 (1300): 研究室お茶部屋． 「北欧」バゲット． トマト．
  • 晩 (1900): ヤマダモンゴル の熟成生ラム． 今日は極上を頼んでみた (一人前 720 円，ふつーのは 520 円)． ここのふつー生ラムでもあっさりしてるのに， これはさらにあっさり度が高い． 野菜． 飯．

2003 年 07 月 20 日 (日)

• 1020 起床． うーむ …… 朝飯． コーヒー．
• でれでれと時間を浪費してしまってから， 腹筋運動ストレッチング． 1345 自宅発北大構内走． いいぐあいに曇天． 1430 帰宅． 体重 71.0kg． 昼飯．
• 1630 自宅発． 曇． ふらふらと札幌駅周辺の書店をさまよう． 何も買う予定もないし， 何も買わなかったわけだが． いやはや． 1800 研究室着．
• ふと「魔がさして」しまって， ghostscript-7 のコンパイル (srpm → rpm) に着手してしまった …… だって， Vine-Linux 3.0 のリリースがどんどん 遅れて いるんだよね．
• 予想どーり， rpm 依存関係解決に苦闘する．
• 特に匿名を希望する料理人氏と大学院生たちが 晩飯を作ったんで， 私もそのおしょうばんにあずかる．

[スイカ名古屋食い]

名古屋ではこのように種子を並べるのが正しい作法らしい．

• 晩飯後もコンパイル・コンパイル． で， できたんだけどいまいちありがたみがない．
• さらに何をとちくるったのか， VineSeed による 3.0α 化に手を出してしまう …… これは初心者のヒトは絶対にやってはいけません．
• 帰宅できたのは翌朝 0600 すぎ． どうにかこうにか仕事できる環境にはなったんで． へろへろ． たしかにいろいろ変わって良いかんぢではあるんだが …… ともあれ Vine Linux モルモット状態だな．
• [今日の運動]
  • 北大構内走 1345-1430． ストレッチング．
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (1050): ヨーグルト．
  • 昼 (1450): うどん．
  • 晩 (2030): 特に匿名希望する料理人氏のつくったひつまぶし (名古屋ふうウナギ飯)． キヌサヤ・豆腐のすまし汁． 野菜あれこれあえもの． ワカメあえもの． スイカ．