ぎょーむ日誌 2006-12-08

2006 年 12 月 08 日 (金)

• 0700 起床． 朝飯． 準備のつづきにとりくむけど， なかなかに進捗しないのですよ． 0850 宿発． 曇． えーと， どこかに有栖川宮公園方面に登っていく坂があるはずだ， と左をみながら歩いていくと …… 昨日の地下鉄広尾駅までもどってしまった． ということで， いつものごとくドイツ大使館と有栖川宮公園の間の坂 (南部坂?) を登攀する経路をとる． 0910 統計数理研究所着．
• 間瀬さん・中野さん・中間さんにごあいさつ． 中間さんに 「R で算時間を要する code を」 (おそらく何かの性能評価用だと思うんだけど) と …… はい， 忘れてたわけではありませんが …… どういうコードがいいのかな．
• 0930 「Rの整備と利用」研究会． 最初は ``Prof'' Brian Ripley 先生， ``Software for Statistical Developments''． 「私はここではソフトウェアそのものの話はしない， それは本日の他の演者がやるから」 というでだしで始まったハナシだけど …… 私にはやや難解だった． 内容的にはぜんぜん難しくないはずのハナシだと思うんだが．
• 強調されていたのはこういったあたりだ． これらは理解できたんだが ……
  • end-users が使う統計学的手法を新しくしていこう (統計学の前線と end-user 世界を近くする)． そのためにはソフトウェアが重要．
  • あるデータ解析のやりかたを考えるのと同じぐらい， そのデータ視覚化の方法の検討は重要であろう (Finding ways to vizualize datasets can be as important as ways to analyze the datasets)． 例として lattice trellis package 使った作図による脳の MRI データ解析． The results can be explained to non-statisticians.
  • 混合モデル (mixed effecs model) の重要性． 上の脳 MRI データの個体差を考慮して library(nlme) で解析．
  • GAM (一般化加法モデル) など ``semiparmetric'' 解析のすすめ． GCV など cross validation 技法とくみあわせて．
  • Classification Tree を open にする過程． 1990 年代後半においてこの問題をあつかうときには CART という (いろいろな意味で) 閉鎖的なソフトウェアが支配的だったので， Therneau たちが rpart という独立した S-plus 用の実装を作った (そして Ripley 先生が R に移植)． 独立に作ってみるともとの CART のバグなどわかった．
  • SPSS も SAS も S-plus も多数の統計学者の verification は受けていない． しかし R は違う． R は統計学専門家たち多数の目にさらされている． ユーザーが bug を見つける → コードが修正される過程もすごく速い．
  ややわかりにくかったのはこれらをつなげているハナシ， とくにハナシの最初につながる導入部のあたり， かな． 「もっとも使われている統計ソフトウェアは何か? ゑくせるである」 といった指摘もあったんだが．

[Ripley 先生]

上の写真は BLAS の後藤さんと話していて， 下は Stefano Iacus さんと． Ripley 先生の印象， がっしりとした体格， 血色よく眉毛ふとく意志の強さがでているかんじ， でかい声で明晰に話す …… にもかかわらず， 講演は残念ながら 8 割ぐらいしか意味が取れなかった． うーむ， 原因かもしれないと思えるものは， r-help 投稿にも見られる「うまい Ripley 表現」 みたいなものに当方がついていけないとか? これと関連してなんだけど， 投影資料 (「パラグラフごとはりつけ型」とでもいうべきスタイル) ながめてると 「慣用句」 の使用頻度がふつーよりやや高くて難しい印象になってるのかも， という気もした． まあ， 私の訓練不足がなによりの原因ですが．

• 二番手はテキサス工科大学から来られた 後藤さん． 統計学のハナシではなく， 最適化 (optimization) のハナシである． 「私は英語はそんなに得意ではないけど， CPU のコトバ (アセンブラ) ならどの CPU のものでも話せる」 というヒトによる CPU の「実験大脳生理学」とでも言うような内容．
• 午前の部は終わり． 谷村さん と昼飯にいき， 空間統計学のハナシなど．
• 1330 から Stefano Iacus さん， イタリアからきた R core team の一人． また optimization ハナシなんだが， こちらは R の中での話なんで さすがに私でも応用手段が考えられる内容でした． いわば ``Effective R 序章'' というような． たとえば， x <- x + i という操作において何回のオブジェクトコピーが発生するか考えてみよ とか， R の関数は pass-by-value なんだからさあとか， また「型変換」でメモリの使用量 (これは gc() その他で調べられる) がぜんぜん変わったり， などなど． ふーむ， external pointer ワザだの， garbage collector よけの PROTECT ワザだの …… ふーむ． これはちょっとふつーのユーザーの範囲をこえぎみですが．
• で， このあたりで「明日の準備がまだまだ不足」 というのが気になって撤退． 1445 宿にもどる．

[有栖川宮記念公園]

今回は統数研の正面 (北側) にある公園の中を突破して戻る． 今朝みつけられなかった「公園に登っていく坂」は 北条坂 と名づけられていた．

• しかしどうもばてぎみで進捗せず． 一時間ほど寝たらかなり復旧した． 作業つづける． 自己紹介の一枚を作るの忘れてた． 何しろ「生態学って何?」なヒトたちがほとんどなので． これを「エンドユーザー」なる keyword の導入とするか．

• そうか OpenBugs は Component Pascal なる言語で実装されてるのか (``BlackBox'' Component Builder)．
• この一枚も新しくする．

• 2300 晩飯に出る． この時間に開店してるところを探すのはめんどうなので， ここと地下鉄広尾駅のあいだにある ジョナサン へ．
• 晩飯くいつつ R News と R2WinBUGS 論文 (Sturtz et al. 2005) 読んでよーやく理解できたんだけど， WinBUGS は 2004 年に開発が停止し， その正式な後継 project として OpenBUGS が始まった， と．

  OpenBUGS is also expected to run under Linux.

  だそーだが， 現状では達成されておらず，

  OpenBUGS is still in development and suffers frequent crashes.

  これは現状でも維持されている， と．
• 2410 宿にもどる． 2730 準備終了． いやはやー
• [今日の運動]
  • うんどうしてない ……
• [今日の食卓]
  • 朝 (0720): バターロール．
  • 昼 (1220): (昨年 も行った) 統数研前坂くだったところにある 元 でゴマだれ鉄火丼． 1050 円． コーヒーつき．
  • 晩 (2320): タラコスパゲッティー． サラダ， ハーブティーセット． 1304 円， 税 + 深夜料金 (118 円) こみ．

[深夜の晩飯]

金曜日夜のぢょなさんは深夜でも客がたくさんいてやかましい． 若い女性も多いけど， (酔っ払ってもいないのに) みょーな説教してるおやぢなんかも． その説教論法は私には「東京ふう」に思える …… my ルールを他人に押しつける， という点において． これはどの地方でも共通なのか?