ぎょーむ日誌 2003-12-(01-10)

2003 年 12 月 01 日 (月)

• 1020 起床． なんでこんなにまぬけな時刻におきてるか，と言いますと …… 明けがたまで額の奇病 (これは本日めでたく疱疹つまり herpes と認定) が痛んで眠れなかったんで， そういや自宅医療箱に抗炎症剤が残ってたなぁ， と医者の処方もうけずに服用してしまったのであった (じつにアブないことを)． これはキくんだけど， 同時にかなり眠気をもよおすクスリなのでありました． ともあれ朝飯． コーヒー． 別の人からも「ヘルペスでは」ご指摘をいただいたので， あきらめて皮膚科に行くことにする． ちっとも回復してないしね． 1110 自宅発． 晴． 今日も気温たかい．
• 1120 北大地環研前の中野医院着． 20 分ほど待つ． 診療室に入ってすぐに「ヘルペスですね」 で診断終了． 帯状疱疹かどうかわからんとのこと． 抗ウィルス軟膏ぬってもらう． 目や耳がやられる可能性もあるから， 気をつけろ (すぐに医者にいけ) とのこと．
• さて， 本日の請求額はじつに思わぬ出費となってしまった …… 10840 円! えーと， 診察料 2800 円・投薬量 31590 円・処置料 1660 円 の合計 36120 円の 患者負担割合 30% で 10840 円， と． クスリの量が多いんだよね ……
  • 外用抗ウィルス剤 アラセナ A 軟膏 (持田製薬) × チューブ 2 本
  • 内服抗ウィルス剤 Valtrex (GlaxoSimthKline) × 一回 2 錠 × 一日 3 回 × 一週間分
• 1200 中野医院発． 1205 研究室着． ネットで調べると， 単純疱瘡は完治しない一生もんのびょーきだそうで． 気長につきあうしかないのはあきらめるとして， 発症するたびに一万円超とられるのはイヤだなぁ． 単に「ヘルペス」といった場合， 帯状疱疹のほうをさすんのか? こっちだと治るとかいう記述もあるな． 原因は子供のころにわずらった (はずなんだが覚えていない) 水疱瘡だそーで． Virus にもいろいろ都合あるんだろうけど， ヒトの神経細胞内で生活するのはヤメてもらいたいんだが．
• まあ， いづれにせよ， あとはヤク物と私の免疫系の活躍に期待しつつ静観， だ．
• 明日こられるんだろう， と思ってた 紺野さん が突然あらわれる (北大で用事があったそうで)． いきなり R の使いかたに関していろいろと質問 …… このあたりの熱心さって， 混合モデル的には個体差で説明できるのであって， 年齢という変数はあまり重要ではないのではなかろーか．
• 紺野さん去られたあとに昼飯． 食後に抗ウィルス剤なるものをのんでみる …… なンかおそろしげな名称だな， と思いつつ anti-virus の生体インストール．
• 授業の準備 …… じゃなくて， 明日の一般化線形混合モデルの informal (informative ではないかも) 解説の準備． これまた， やはりぢりぢりとしかススまん． ということで， 私は説明の準備がてきぱきできないヒト， という事実が判明しつつあり．
• 先週で生物多様性論 I が終り今週から II が始まる． これはもちまわり授業で今回次回は甲山さん． C-104 講義室の設定はわかりにくいんで初回だけ手伝う． PowerBook G4 がなかなか液晶プロジェクター認識しなくてアセる．
• informal 準備の続き．
• 終らん． ばてた． 2030 研究室発． 2100 帰宅． 体重 71.4kg． 運動不足ふとりか? 晩飯． またクスリ． 服用すんのは気分的にイヤなんだけど， たしかにこれ効きますな． 眠くなるとかの副作用もないみたいだし (胃や腎臓は痛めつけられるらしい)．
• informal 準備 …… はススまず． また別のかたからヘルペス情報いただく． ありがたいことである． 2430 消灯．
• [今日の素読]
  • ヘルペス呪われってことで
• [今日の運動]
  • こっちも
• [今日の食卓]
  • 朝 (1030): 米麦 0.7 合． シュンギク．
  • 昼 (1530): 弁当． 研究室お茶部屋． 米麦 0.7 合． チンゲンサイ・タマネギ・エノキダケの炒めもの．
  • 晩 (2140): マカロニ． トマト・タマネギ・ブナシメジ白・豆腐・煮干のシチュー (?)

2003 年 12 月 02 日 (火)

• 0850 起床． 朝飯． コーヒー． 0930 自宅発． 晴． 0940 研究室着．
• informal 解説準備のつづき． スケジュールの関係上， 長びいてはいかんのでともかく内容を限定する．
• 内容限定してるおかげで 1300 になんとかひととーりの準備できた． 昼飯． informal 解説まで混合モデル文献みなおし．
• 1500 より informal な 講座セミナー で一般化線形混合モデル (GLMM) の簡単な解説． なぜかあちこちからぞろぞろと聞きにきたヒトたちが． 解説に使ったファイルは統計学 授業 ペイジに番外篇としてつけたしておきました．
• 1550 予定どーり早めに終了． ばてた． もうしばらくは統計学授業のマネごとなんぞやらんぞ．
• 1630 本日のメイン event たる Trendy セミナー (私のハナシはこれの前座であり，来週の講座セミナーの前座)． 今日は帯広畜産大 (紺野研) の M2 佐藤さんが樹木の繁殖とそのサイズ・明るさの関連を調べたもの． いつものことながら， 順序変数 (ここでは「明るさ」) のあつかいは難しい．
• セミナー後に苫小牧から来られた鍋嶋さんと直径成長モデリング雑談少々． 私が考えるようなアヤしげな方法に かなり不安感をもっておられるようなので (とモデル業者は顧客の顔色から内部状態を憶測する癖がある)， とりあえずはふつーっぽい推定計算からスタートしたほうがよいんでは， という方向に． なるほど， そちらのほうが無難でしょう． じつはそのふつーっぽい方法の延長として アヤしげなやりかたが控えていたり ……
• 2040 研究室発． 2050 帰宅． 体重 71.4kg． 晩飯． 抗ウィルス剤なんぞ服用すると， アタマの働きがいつもより悪くなるような気がする， と称してナマける．
• 2330 消灯．
• [今日の素読]
  • Anti-virus ばて …… などとほざいていると， なかなか再開できませんなぁ
• [今日の運動]
  • こちらも
• [今日の食卓]
  • 朝 (0910): 米麦 0.7 合． 目玉焼き．
  • 昼 (1320): 弁当． 研究室お茶部屋． 米麦 0.7 合． トマト・タマネギ・ブナシメジ白・豆腐・煮干のシチュー (?)
  • 晩 (2130): 米麦 0.7 合． トマト・タマネギ・ブナシメジ白・豆腐・煮干のシチュー (?)

2003 年 12 月 03 日 (水)

• 0820 起床． どうも 抗ウィルス剤 をのんでるとやたらと寝てしまうような気が …… そういう副作用については書かれていなかったんだが． 朝飯． コーヒー． 0910 自宅発． 晴． ちっとも冬らしくない． 0920 研究室着．
• 午前中は Matt の熱帯林樹木死亡率モデリングの こんさるたんとで終った． かなりばて． 例によって例のごとく近傍個体との相互作用とかを， どう定式化するか， というようなハナシ．
• 昼飯． 食後に散歩がてら北大生協に． 講座お茶部屋のコーヒーなどがキレてたので， そのおつかいということで．
• 午後はなぜかマルハナバチ論文の統計学解析の解読に取り組む． これは M1 平林さんが「何やってるのかわかりません」 質問されたんで， なーにこんな G 検定なんぞ簡単簡単 …… と読みはじめたんだけど， じつに何をやってるのかさっぱりわからない． 数理的なデータ解析がこんなに理解できない論文はひさしぶりだ …… 統計学の適用とかが不適切なだけの論文だったら， 読めばすぐに 「あ，またまた○○ってパターンの間違いやってるな」 とわかることが多いんだけどなあ ……
• 夕方すぎまで解読に取り組んで， よーやくこの著者たち何をかんちがいしてんのか， ってのがわかってしまいましたよ …… もはや感情的な悪口を述べてしまいそうなので， 出典は書かないけど American Naturalist に今年掲載された論文である．
• American Naturalist って誌名みただけでは 何やらわからんでしょうけど， 生態学・進化生物学の論文のせてる学術誌で， 「高邁さ」をウリにしてることになってる． 読むほうとしては， まさかあの American Naturalist で， うるさそうな査読者が 2-3 人アラさがししたあとで， かくも完全意味不明なる変てこデータ解析作法が残存してる， とは考えないものだ． ということで， American Naturalist 掲載だからって， それだけでは何のアテにならぬものである．
• 解読不可能だったのは， 説明がまったく書かれていなかったこと． さらに説明がないのに heterogeneity なる語を集団とか環境の「不均質さ」とはまったく無縁な， 生物学的にも統計学的にもおそらく前代未聞なる 意味不明な使いかたしやがった --- データを「プールしてない部分」の意味なんですよ． 我ながら， たった数時間ばかり苦しんだだけで， よくもまあこんな正気などというものは激烈に飛び越えた語釈 に成功したもんだと， 感心してしまいますよ． いやはや． (後記: 少なくともこの記述は 完全に間違い --- Sokal & Rohlf (1995) section 17.3 は heterogeneity G-test あるいは replicated G-test に関して丁寧に説明しています)
• かかる難解論文はなぜ生じてしまったのか?
  1. 著者は統計学って何をやるのかよくわかっていない
  2. しかし創意工夫の才能はある (これが無いヒトは， とてもわかりやすい誤用をするので， 読むほうとしては助かる)
  3. すると一見もっともらしい外見の解析結果だけど， よく見るとまるで意味不明なデータの流れが生じる
  4. 著者もうすうす「これヘンだな」 と気づいてるので Method 節などから必要な説明すらもすっとばされる ……
  査読者どもは理解できてないくせに， わかったフリをして OK 出した公算大． 素直に「わけわからん」とか言えよな …… 書いててますます気分わるくなってきたので終了． (後記: 例によって， このあたりの指摘はかなりまと外れなんでは， とご教示いただく …… 明日以降のぎょーむ日誌参照)
• 実験とかについてはまっとうな論文だったのかもしれないけど， せっかくそういうハチ訪花の室内実験をきちんとやって 莫大な観測データを得ておきながら， こういうわけわからん解析やるもんかねぇ …… なんというか， 苦労してかせいだお金を怪しげなる新興宗教にじゃんじゃん貢いでる， というような構図といえばわかりやすいだろう．
• というようなハナシを質問者たる平林さんにやって一件落着． ふつーにやればどうということのない推定・検定にすぎないんだよね …… 例によって「検定やる状況ではない」に該当すると思うんだけど． あまりにもばてきってしまったので， お茶部屋でごろごろと転がっていると， 来週に博士論文予備審をひかえ ますます反社会性を亢進させつつある某氏 (他講座) がやってきて， 予備審関係者の搭乗する航空機その他の遭難を神仏悪魔に祈念したり， 自主的なるテロ活動にて混乱惹起させておいて予備審をうやむやにせん と画策する凶々しい現場に立ち会ってしまって， ますますばて感が悪化していく．
• どーにかこーにか起き上がって 2030 研究室発． 2050 帰宅． 体重 71.4kg． 晩飯．
• で， 今度はしばらくおいておいた， 某こんさるがらみのサカナ統計モデリング (無料御奉仕期間中) に関するメイル書き．
• 今月末に四国にお邪魔しようかなと小林さんとご相談．
• 今日は， 熱帯林・マルハナバチ・サカナの統計学， で一日が終ってしまったなぁ． 2430 消灯．
• [今日の素読]
  • ぅ-
• [今日の運動]
  • ちょい様子見中 ……
• [今日の食卓]
  • 朝 (0840): 米麦 0.7 合． ワカメの味噌汁．
  • 昼 (1300): 弁当． 研究室お茶部屋． 米麦 0.7 合． トマト・タマネギ・ブナシメジ白・豆腐・煮干のシチュー (?)
  • 晩 (2130): 米麦 0.7 合． ハクサイ・タマネギ・エノキダケ・鶏レバの炒めもの．

2003 年 12 月 04 日 (木)

• 0700 起床． 朝飯． コーヒー． 0830 自宅発． 晴． しかし今シーズン初の真冬日になりそう， ということで寒い． 0840 研究室着．
• 皇太子妃殿下も 帯状疱疹 ということで， 「へるぺすにやられました」 と言ってもみょーな顔されなくてもすむかもしれない． いままで周囲の人々からは 魚類のびょーき あつかいだったからなぁ ……
• 予算消化に必要なる作業あれこれ …… 買っておきたい本の検索とか． しかし amazon.com って 2005 年出版予定の本の予約とか受けつけてるわけ? 著者はまだ一行も書いてないんでは．
• きりがないのでてきとーにして， しばらく放置しておいた釧路大会雑用にとりくむか， と思ったら苫小牧の村上さんからメイル． この時期に村上さんから， というと一年前の年末年始におよんだ悪夢のごとき 借金とり追求生活が思い出されるわけだが …… R の MacOS X Aqua 版の 1.8.1 がなかなかよろしい， との近況メイルだったのでひと安心． 私も試みにお茶部屋 iMac にインストールしてみる (ダウンロードサイトはここ) …… たしかに package 導入とか， かなり簡単かも． しかし同じ MacOS X 用でも X11 版のほうがいいなぁ， などと考えてみたり． これら二つは整合性とかとれてないのか?
• にしても， MacOS X の描画はおそい --- RAqua も X11 で描画してるからだろうか?
• CGI ぎょーむのつづき．
• 年末香川旅行の航空券の手配， とか．
• かとー先生にいろいろとご教示いただき， MacOS X 上の X11 ターミナルからでも /Applications/StartR.app/RAqua.app/Contents/R で RAqua 起動できるとわかった． .bashrcとかに alias で設定すればいいんだが ……
• MacOS X 上では chsh できん． これまたかとー先生に教えていただいて， Netinfo マネージャ なるユーティリティ使わんといかん， とわかった． うーむ， 不便な Unix 機め ……
• また釧路大会ぎょーむに取り組む ……
• 粕谷さんから親切なるご教示メイルいただく …… いやー， またまた大誤爆というか大誤読していましたなぁ． われながら， やれやれ， だ． 昨日のぎょーむ日誌にとりあえず 後記 をいろいろと入れてみる．
• どこから説明とゆーか釈明とゆーかおわびを始めたものか …… えーとですね， ``heterogeneity'' なる語ですが， これは replicated G-test においては当たりまえのごとく使われる語です． ちゃんと， Sokal & Rohlf (1995) の Box17.4 に懇切丁寧に説明されてます (しかも粕谷さんによると Sokal & Rohlf 以前の太古の昔から 使われつづけている， とのことで……)．
• もうちょいだけ補足しておくと …… 反復ごとに G を計算してそれを足し合わせたものが GTotal， つぎに反復を「まとめ」てから計算したのが GPooled， でこれらを引き算して得られたのが GHeterogeneity = GTotal - GPooled となる．
• いやー， もと論文では完全ランダムという帰無仮説のもとで replicated G-test やってるのはわかるとして， 特定の花は避けるというのも帰無仮説にしちゃってるんだけど， こっちでは GPooled だけ しか計算してないんだよね (しかも事象生起確率ゼロはイヤなので 0.01 という値使います， とかいうヘンな尤度計算)． このあたりにだまされてしまった …… 論文中には書かれていないので， GPooled だけが計算されている理由を憶測すると， GPooled ≤ GTotal という関係が常に成立いるから， GPooled が十分に大きければ他の G は計算しなくていいでしょう， というロジックであるよーだ．
• つまり， 誤読した私も不注意きわまりなくてよろしくないけれど， ちょー不親切な著者たちも少し責任を分担してくれたらなぁ …… と． いやいや， すべては herpes virus が悪いのです， ということにしておくか．
• 実際のところ， Sokal & Rohlf (1995) とかに載ってるよーな， 実験科学者にとって有用な方法， というのはほとんど勉強してなかったことに， 今回あらためて気づかされてしまった． いやー， 植物生態学の野外調査データなんて ホントの意味で replication と呼べるものなんて無いんですよ …… などと放言してみたら， こんどは植物のヒトたちに怒られてしまうんだろうか?
• ともかく夕方から夜にかけてそのあたりの勉強で終る． 昨日に， ウソ講釈をしてしまった M1 平林さんあいてに， replicated G-test の再履修的な講釈 (講釈師がわの再履修)．
• 帰ろうとしたら， なにから良からぬことを企んでる大澤君のあやしげ シウリザクラわーるどに捕獲されてしまって ……
• 今日もばてた 2245 研究室発． 2300 帰宅． 体重 71.2kg． 晩飯． あいかわらず抗ウィルス剤 Valtrex をのむ日々 …… これって「特効薬」(とうたうわりには効き目がゆっくりしてるよーな) ってことになってるから， 皇居の中でも服用されてるのかもしれませんなぁ．
• 2430 消灯．
• [今日の素読]
  • 誤爆誤読は時間のロス
• [今日の運動]
  • なんか， Valtrex のんでると主観的には 病人モードになってしまって 「運動したらいかんのでは」 と思ってしまうんですよね …… さっさとクスリやめたい
• [今日の食卓]
  • 朝 (0730): 米麦 0.7 合． ハクサイ・タマネギ・エノキダケ・鶏レバの炒めもの．
  • 昼 (1300): 弁当． 研究室お茶部屋． 米麦 0.7 合． ハクサイ・タマネギ・エノキダケ・鶏レバの炒めもの．
  • 晩 (2330): 米麦 0.7 合． ハクサイ・マイタケ・豆腐・煮干の味噌汁．

2003 年 12 月 05 日 (金)

• 0750 起床． 朝飯． コーヒー． 0900 自宅発． 曇． 0910 研究室着．
• 1000 論文セミナー， 本日の担当は M1 立花さんで， 量的遺伝モデルで説明されそうな植物の資源配分可塑性， なかでも環境によって配分の重点を変更できるかどうか， ということを調べた論文の紹介． Allium Vineale なるユリ科の 植物は， 両性花・地上部に形成されるむかご (bulbil)・地下部の小鱗茎 と三つも繁殖方法があり， これから genet ごとに栽培実験してみると …… という内容． 環境応答的な資源の傾斜配分ってのは なかなかデータからは特定できん， という結果． つまり， まぁわれわれモデル屋が妄想するよーな 安直なスイッチングを現実の中から探しだすのは， 手間かかるだけでなくなかなか困難そう， ということがわかった． うん， 昨日なやまされたハチ訪花論文ぢゃないけど， 動物の行動とかだともう少し人間っぽいのかもしれない．
• 昼飯後， 釧路大会 CGI プログラム改造ぎょーむ． 竹中さん たちが第 50 回 (つくば) 大会用に作成されたのを直してるだけなのに， えらく手間どるな．
• 1600 釧路大会の旅行関係手配・エクスカーション・懇親会は 近畿日本ツーリストという旅行代理店にお願いしてるんだけど， その担当者のかたがこちらまで来てくださる． とはいっても， 札幌で開催された会議のついで， だそーで． ネット雑用掛下請けたる私あいてなので， 先方の web site との関連についての連絡・確認あれこれ． エクスカーション …… たしかにコースが 5 とーりもあるけど， 区別つきませんなぁ． そもそも生態学研究者ってのは， 大会終了後に釧路湿原の半日見物行とかに参加するもんなんだろうか?
• しかし， 大学ってのはつくづくお客さん迎撃むけではないな …… あ， 院生その他の飼育小屋的な 8F ではなく (私が居候してるかとーオフィスについても社会逸脱度は低いと言えぬ)， 7F の甲山さん部屋横の図書室兼応接室とか使わせてもらえばよかった (雪野さんが片づけておられるので)．
• また釧路大会参加申し込み CGI プログラムいじり． 私の場合， プログラミングは嫌いでも苦手でもないんだけど， こういうユーザーインターフェイス (しかも不特定多数が使うもの) が重視されるやつは， じつに慣れていなくてしかもやっていて楽しくない， ということがわかってきた． うう．
• 2200 研究室発． 2220 帰宅． 体重 71.6kg． 運動不足ぶとりか? 晩飯．
• 2400 消灯．
• [今日の素読]
  • スケジュール破綻な日々
• [今日の運動]
  • やはり Valtrex のんでると， カラダがびみょーにだるい
• [今日の食卓]
  • 朝 (0820): 米麦 0.7 合． ハクサイ・マイタケ・豆腐・煮干の味噌汁．
  • 昼 (1240): 弁当． 研究室お茶部屋． 米麦 0.7 合． ハクサイ・マイタケ・豆腐・煮干の味噌汁．
  • 晩 (2300): 米麦 0.7 合． レトルトパウチドのカレー． さつまあげ．

2003 年 12 月 06 日 (土)

• 0940 起床． 抗ウィルス剤 Valtrex ばて (?) で眠りこけてた． 朝飯． コーヒー． だいぶ具合もよくなってきたので， Valtrex やめてみる． ヘルペス用軟膏はつけるけど， ばん創膏やめ． 傷口も完全にふさがってるんで．
• いつものごとくでれでれと．
• 1430 自宅発． 曇． 1440 研究室着．
• また釧路大会参加うけつけ CGI プログラムの改造ぎょーむ．
• 小菅せんせーとあやしげ音楽ファイルやりとりのため (と勝手に断定)， 井田君が来ていたので統計学バナシなど． 2210 研究室発． 雪． これは積もりそうだ． 2230 帰宅．
• 2430 消灯．
• [今日の素読]
  • とほほ
• [今日の運動]
  • まだ復活できづ ……
• [今日の食卓]
  • 朝 (1010): 米麦 0.8 合． 干しシイタケ・ワカメ・豆腐の味噌汁． スクランブルドエッグ．
  • 昼 (1330): 米麦 0.8 合． 干しシイタケ・ワカメ・豆腐の味噌汁．
  • 晩 (2310): 米麦 1.0 合． 干しシイタケ・ワカメ・豆腐の味噌汁．

2003 年 12 月 07 日 (日)

• 0920 起床． 朝飯． コーヒー． 積雪 20cm というところか．
• でれでれ …… してるわけにもいかんので， CGI プログラミングぎょーむの続きとか． ユーザーインターフェイスの手なおし ってのはキリがないからイヤだ．
• 1610 自宅発． 曇． 歩道上も部分的に凍結． 車どもがよたよた走っていて， じつに危なっかしい． 1625 研究室着．
• また北大ネット (HINES) の POP サーヴァーが落ちてやがる …… なんで， この連中はメイルサーヴァーぐらいまともに管理できんわけ?
• 日曜日だとぜったいになおさんし …… さらに問題なのは 月曜日になったからといってなおる保証はどこにもないことだ．
• とぶちぶちと文句いいながら釧路大会参加申し込みペイジの 敵国語版などつくってみたり． 講座トップペイジで使ってる jesplit.pl なる Japanese-Enemy lang 分離 HTML ファイル生成プログラム． 自作ながら， これ便利なのか不便なのかいまいちよくわからん．
• ペイジなど書いたりしてると， あほあほな間違いにいろいろと気づく．
  • ``51th Annual Meeting'' って何だ? ``51st'' だな …… ちくせう敵国語め
  • ナゾの一日: 4/30 までは 6000 円， 4/30 以降は 7000 円 …… じゃあ 4/30 はどっちなんだ?
  しかしメイルが使えんので連絡とれん．
• しかしまあ， これほどメイルサーヴァーがあてにならん， という状況では参加申し込みの複写をメイルで取る， という手法はあまり有効でないよーな． どうしたもんかのう．
• えーと， 明日に D 論予備審のある谷さんは， 予備審対策弁当を作っていただいたそーで． いいのう ……

[予備審対策弁当]

谷さん公認撮影のもの． アサリやホウレンソウといったビタミン・ミネラルの摂取に 細やかな配慮のうかがえるもの …… なお当方では左の写真の高解像度版もきーぷしております． このお弁当を作ったヒトが将来なにか 有効活用 していただける機会もあるかもしれませんので．

• どうでもいいが， 釧路大会ぎょーむに時間をとられすぎている． 2220 研究室発． 歩道はばりばりに凍結． 2240 帰宅． 晩飯． 除雪車の音が外から聞こえてくる．
• 2500 消灯．
• [今日の素読]
  • Salsburg, D. 2001. `` The Lady Tasting Tea -- How statistics revolutionized science in the twentieth century''. Owl Book.
  • Chapter 21. A Genius In The Family
    • I.J. Good
  • One of I.J. Good's major achievements was the contribution he made to the development of empirical Bayes and hierrachal Bayes methods which he derived from the work he did at Bletchley. He emerged from his wartime experiences with a deep iterest in the underlying theories of mathematical statistics. He taught breifly at the Unieversity of Manchester, but the British government lured him back into intelligence, where he became an important figure in the adoption of computers to deal with problems of cryptanalysis. The power of the computer to examine vast numbers of possible combinations led him into investigations of classification theory, where observed units are organized in terms of different definitions of ``closeness.''
• [今日の運動]
  • 雪つもっちまいました ……
• [今日の食卓]
  • 朝 (0940): シリアル．
  • 昼 (1430): 米麦 0.7 合． 干しシイタケ・ワカメ・豆腐の味噌汁．
  • 晩 (2330): 米麦 0.7 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．

2003 年 12 月 08 日 (月)

• 0910 起床． ヘルペス対策の抗ウィルス剤 Valtrex は夜だけはのんでるんで …… とゆーのは遅おきの言いわけとして利用できないんだろうか? 朝飯． コーヒー． 0945 自宅発． 曇． 0955 研究室着．
• 1015 より谷さん予備審． 1140 ごろ終了． まあ， 意外と無難におわったといいますか．
• しかしそのあとお茶部屋でひたすらデータ解析ぎろんが 1330 ごろまで． これまた気象データ解析問題だ．
• とりあえずメイルサーヴァー復活はしてるようだが ……
• またまた釧路雑用ぎょーむ． いいかげんに終らせたい．
• 気分転換に Vine Linux の apt-get． 今回はめずらしく 10 コぐらい package の入れ換えが発生した． くだんの rsync がらみ， とかですな．
• 参加登録もうしこみの複写をどう取るか再考する． 竹中さん方式では Cc メイルだ． しかし北大のメイルサーヴァーはかくのごとく信用できん． さらに容量制限がいよーにキビしい． scp とかどうだろう? ssh で パスワードなしでバックアップとったりする場合には， パスフレイズなしで設定すればいいわけだが (参考 ) …… CGI プログラムを駆動してる架空ユーザー apache は ssh 使えんな．
• ちょー大容量メモリをちょー高速検索できるかとー先生に ご教示いただいたところによると， apache には SuEXEC なるオプションがあって (しかし， まあ， いつものことながら， かとーさんはなぜ自分が使わない機能あれこれも知悉してるのだろうか?)， これは個別ユーザが CGI プログラムを実行できます …… というもの． つまりは apache のコンパイルやりなおせばいいだけなんだが …… たかが参加登録複写のためにそこまでやるか?
• 結局はやはりメイル方式に． ただし北大には送らん (釧路大会用メイルアカウント @ees.hokudai.ac.jp はすでに取得してるんだけど)． 私が加入してる ISP である asahi-net ではメイルボックスが複数つかえて， しかも一個目の追加は無料だ． asahi-net だと北大よりあきらかにサーヴァー事故は少ないし (とゆーか，北大ネットみたいにコケまくってたら廃業ですな)， 格納できるメイルの量は無制限だ (とはいえ，予定としては合計 0.1MB 程度にしかならんはず)． これを申請してみる …… アカウント発行は明日になるのか．
• 複写用のメイルボックスを学外に出したところで， 北大ルーターや言うところの「メイルゲイトウェイサーヴァー」 なるものがコケたら連絡不能になるところは依然として同じなわけだが ……
• ところでメイルアカウントだけ を Perl の正規表現で match させるにはどうしたらよいだろうか? いただいたプログラムでは /^[^\s@]+@[^.\s]+\.[^\s]+$/ というチェックが入るようになってる． これでかなりのでたらめ入力がふせげるんだけど， k*u*b*o@ees......hokudai#./ac\.jp!!!... というようなわけわからん文字列も通過してしまう． ちょっと改良を試みてみる …… 現実にメイルアカウントに使われるのは Perl で言うところの alphanumeric にくわえてせいぜい _ とか - だろう， と強引に仮定して， /^[\w-\.]+@[\w-]+\.[\w-\.]+\w$/ てなところでどうだろう．
• もちろんこれでもアナだらけなんだけど (ees.....hokudai などはとおる)．
• で， いまふと思いついたんだけど， /^[\w-\.]+@[\w-]+(\.\w+)+$/ だと少しマシになるのかな． ees.....hokudai.ac.jp はさけられる． ku..............bo@ees.hokudai.ac.jp は通してしまうけど …… よく知らないんだけど， RFC 的にはこういうヘンなアドレスも OK なのか? (後記: これではダメで /^[\w\-\.]+@[\w\-]+(\.[\w\-]+)+$/ とする …… なお - にもエスケイプが必要みたいだ)
• また北大ネットのメイルサーヴァー不調 …… 「温厚なる」かとー先生もお怒りで「だめ」出し．
• はらへったので撤退． 1930 研究室発． 1940 帰宅．
• 季節変動する気象データのことを考えて眠れなくなってしまった …… 2630 消灯．
• [今日の素読]
  • Salsburg, D. 2001. `` The Lady Tasting Tea -- How statistics revolutionized science in the twentieth century''. Owl Book.
  • Chapter 21. A Genius In The Family
    • Persi Diaconis
  • Suppose we prodce 500 measurements at each clinic vist for each patient, and taht the patients are sent 10 different visits over the cource of the study. Each patients provides us with 5,000 measurements. If there are 20,000 patients in the study, this can be represented as 20,000 points in a 5,000-dimensional space. Forget the idea of travel through a mere four dimensions that preoccupies much of science fiction. In the real world of statistical analysis, it is not unusual to have to deal with a space of thousands of dimensions. In the 1950s, Richard Bellman produced a set of theorems, which he called ``the curse of dimensionality.'' What these theorems say is less possible to get good estimates of parameters. Once the analyst is out in ten- to twenty- dimensional space, it is not possible to detect anything with less than hundreds of thousands of observations.
• [今日の運動]
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0920): 米麦 0.7 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．
  • 昼 : なんとなく食ってない．
  • 晩 (2020): 米麦 1.0 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．

2003 年 12 月 09 日 (火)

• 0830 起床． ねむい． 朝飯． コーヒー． 0910 自宅発． 曇． 0920 研究室着．
• メイルサーヴァーは稼働してるもよう． 小林さん からメイルいただく …… それによると， 私は 12/24 に香川大学農学部で クリスマスツリー の作りかた， というようなハナシをすることになってるよーだ (pdf, 42KB)．
• たいへん不思議なことに， いつまでたっても終らぬ釧路大会ぎょーむのつづき．
• 院生などを実験台にして， 入力する様子を観察していると「山田(空白)太郎」 というふうに， こちらがとくに指定してないのに， 勝手に (空白) いれたがるやつがいる， とわかった …… このよーな勝手はなんとなく腹だたしいので， それならばいっそのこの姓・名分断してやるか， というような改造．
• 1300 講座セミナー， 本日は M2 大澤君． 実に 260 分にもおよぶ大長篇セミナーでした． しかしまあ， あの複雑怪奇なシウリザクラわーるどを一望できた， ということでよかったんじゃないでせうか． シウリザクラ苦闘もひとくぎり， か． 食害によって発生する opportunity cost の計算とか， 少し訂正せんといかんけど．
• ともかくクローナル樹木に対する過大な期待， 何かある種の実験に適した生物なんでは， という幻想を木っ端微塵に打ち砕いてくれた， 痛快なる研究でありました． 昨年度の経験とあわせて考えるに， 木本で資源の分配とか転流とかの研究なんぞやっちゃいけません， てことか?
• セミナー後， 講座見学にきた人のハナシを拝聴する， など．
• で， 昼飯も食ってないのに夕方というか日没になってしまって， ひどくばてて， ですね． 二日連続ひるめしぬきだ．
• ばててるのに， 釧路ぎょーむの続き． Asahi-net における第二のメイルボックスが使えるようになったので， そちらに複写を送るようにする． 公募式シンポジウム・自由集会の web うけつけ， ひどくいいかげんなものを構築する．
• お， 2/3 に開催される D 論本審の案内が． 午後の最初に谷さんで， その次はなんと特に匿名を希望する料理人氏である． この人物は言うまでもなく 8F 専属の料理人なんだが (ちかごろ怠業ぎみとの指摘もあり)， 独法化前の大学院では講座専属料理人 (無給) の存在は公的には認められてないから 地環研建物への出入りが難しかろう， ということでとある他講座の大学院生をよそおっている． で， その偽装をいよいよ完ペキなものとすべく， 学位まで取るわけだな． なぜか (なぜ「なぜか」なのかは説明しない)， 昨日ではなく本日に予備審あったらしいンだが …… 昨日の段階では「何も準備してない」 とつぶやきつつ建物内を徘徊しており， 昨晩はお茶部屋で谷さん慰労会で酒を飲みつつ かとー先生に「だめ」出しされてた， というウワサもあり． ともあれ本審にまでこぎつけてめでたいことである．
• あと染色体研の長谷川理さんも本審のよーで． 育児におわれながらのたいへんな D 論書きだった， とか．
• 2045 研究室発． 2100 帰宅． 体重 70.8kg． 昼飯ぬきやせ． 晩飯．
• ばてて何もかもすすまん． 2300 消灯．
• [今日の食卓]
  • 朝 (0840): 米麦 0.6 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．
  • 昼 : 食ってない．
  • 晩 (2130): 米麦 1.0 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．

2003 年 12 月 10 日 (水)

• 0830 起床． 朝飯． コーヒー． 0920 自宅発． 雪． 0930 研究室着．
• PD 石井さんが二次元正規分布問題をもちこんでくる． 標本平均のペアがまた別の二次元正規分布から得られた， という設定のちょいめんどうそうな問題 …… でもじつは簡単だった． こういう状況であっても， ある変量の標本平均の推定値は他の変量に依存しない， ってことで複雑化が避けられていると思う (推定問題の種類によっては， このような割り切りができない)．
• 釧路大会雑用ぎょーむの続き． 1400 ごろまでかかって， うけつけ状況を Web で見れるよーにするプログラム作る． こういうかんぢ． むろん公開するつもりだ． 申込者が各自で勝手に確認してくれ， ということで．

• crontab コマンドで一時間に一回， 上のリストをサーヴァー上で更新するようにする．
• 一区切りついて昼飯． 二日前に食べられなかった弁当． 研究室の冷凍庫で凍結させておいたもの． マイクロウェイブは氷結した水分子を振動させることができない …… 考えてみればあたりまえか．
• 釧路大会実行委員会に「参加申し込み」の実験動物になってほしいんですが， というメイルを送ってみる …… 誰もやる気ナシのよーで． やれやれ． ということで， モルモット志願者ぼしうちう．
• そっちは放置しといて， 大会案内のペイジを作っていく． Wiki の生成する HTML は良いのか悪いのか …… けっきょくかなりのタグを取り換えてしまったわけだが． しかし， まぁこれは MS-Word ごときがげろげろと吐きだす絶望的に意味不明な HTML ファイルにくらべれば 10¹⁰¹⁰ 倍ほどマシだが．
• 夜になってちょいカエデ論文ぎょーむ． gnuplot 日本語版が eps ファイルに混入させる複合フォントは まぬけ ghostscript-5.10 には負担が重すぎる …… set locale '' したり eps ファイルを書き換えるスクリプト作ったり．
• 本文はスペルチェックやりなおしだけで終ってしまった． そして， いまだ釧路ぎょーむ終らず …… 2040 研究室発． 2110 帰宅． 体重 71.4kg． 晩飯．
• どうもさえんなぁ． 2500 消灯．
• [今日の素読]
  • うー
• [今日の運動]
  • 腹筋運動 30 × 3 回．
• [今日の食卓]
  • 朝 (0840): 米麦 0.7 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．
  • 昼 (1420): 弁当． 研究室お茶部屋． 米麦 0.7 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー．
  • 晩 (2140): 米麦 0.7 合． キャベツ・タマネギ・ジャガイモ・ブナシメジのカレー． カレー終了．