ぎょーむ日誌 2003-12-04

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 合． ハクサイ・マイタケ・豆腐・煮干の味噌汁．