ぎょーむ日誌 2007-02-(11-20)

2007 年 02 月 11 日 (日)

• 0730 起床． だらだらとアミノ酸問題にとりくむ． 朝飯． コーヒー．
• しつこく豊平川アミノ酸濃度の階層ベイズモデリング試行錯誤 …… これは今までやってきたのとずいぶんと異なるハナシ， のようで． この測定前処理依存なえらーってのは ……
• 1700 自宅発． JR 札幌駅まわりふらふらしてから 1830 研究室着．
• アミノ酸問題， ちょっと気分をいれかえてかなり大規模なリストラを試みる． この実験系で生じている現象から時系列的な何か を推定するのは不可能ではないか， と ……

• そっちの方向への改造によって見とおしがだいぶ良くなってきた． 2110 研究室発． 寒い． 2140 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • またまたうんどう休養日 ……
• [今日の食卓]
  • 朝 (1030): 自宅から最寄りのパン屋 (しかしなぜか今まで行ったことなかった) 「イソップ」のパン． ここの クロワッサン はたしかにすごい．
  • 昼 (1430): 「イソップ」パンつづき． 海鮮焼きそばパン． バタール．
  • 晩 (2210): 米麦 0.7 合． 麻婆豆腐． ネギ・マイタケの吸い物．

2007 年 02 月 12 日 (月)

• 0720 起床． 朝飯． コーヒー． 洗濯しつつ計算試行錯誤をススめる．
• 木金土日月と雪山ラッセル連続アタック & 連続撤退的な 「豊平川アミノ酸濃度の階層ベイズモデリング」 試行錯誤， よーやくオワリが見えてきた ……
• 昼飯の準備しつつしあげの計算 12000 MCMC step とりぷるチェイン 620 秒 …… まあ， 収束とかまだまだがたがたしてるところあるけれど， 各種ややこしいばらつき錯綜している この観測データセットをシメあげて 人間が引っぱり出せる情報はこんなもんでしょう． 計算終了． 昼飯．

• 今回の解析あれこれ．
  • 水槽効果とでもいうべき random effects があちこちに入っている (図を描くとよくわかる)
  • さらに， 測定機器にかける前の処理によっては 測定値がすごく下ることがあり， その変化はかなり大きい (これも統計モデルに組みこむ試行錯誤が必要であった)
  • 以上のような状況なので 水槽ごとのアミノ酸組成の時間変化， みたいな細かい変化は抽出できないはず …… なンだけど， この組成変化をあらわすような パラメーターが必要である (なぜそうなのか理由はよくわからないんでけど MCMC 計算における混交ぐあいがかなり改善される …… 「動かないパラメーター」をゆさぶるバネとして 機能している?)
  • 今回も多変量確率分布の事前分布は使わなかった …… アミノ酸濃度間の分散共分散 (17 + 17 * (17 - 1) / 2 = 153 個) ぜんぶを推定できるデータではなさそうだったんで …… というか， 多変量事前分布の使いどころをまだまだ理解できていない
• …… てなコトを考えつつ 1340 自宅発． 吹雪だけど晴． 気温は低い． 1355 研究室着．
• 脱力したせいか仕事ススまん． 窓の外は猛烈な吹雪． と思ったらいつのまにか晴れていたり．
• アミノ酸問題に没頭というかかまけていてこの数日間を費してしまったので …… やるべきことは多々あるんだけど， とりあえず 2/15 (木) 〆切の査読原稿にまずはとりくんでみる． こんな状況でどうして査読なんか引き受けるんだ， と思われるかもしれませんが， 世の中なかなか断わりにくい依頼というものがあって， ですね ……
• 先日 30 分ばかりざっと見て (このやや奇妙な) 計算の核心部分は理解できたんで， 本日はなぜかくのごとき状況になってしまったのか， あたりを解明する読解．
• …… するつもりが， 1.5 時間ほど熱帯林解析ハナシ． このハナシ， 次に何をやるにしても計算が格段にたいへんになりそう．
• 査読原稿よみ …… ああ， (私のところにまわってくるだけあって) コレもけっこうハズレだ． つっこみどころ満載． アタマくらくらする． さてさて， どうしたもんだか．
• 1920 研究室発． 1940 帰宅． 体重 68.6kg． 晩飯の準備． 晩飯．
• ちょっとデータ解析こんさるメイルかき， とか．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0800): 米麦 0.6 合． 卵焼き．
  • 昼 (1310): 明太子スパゲッティー．
  • 晩 (2130): 米麦 0.7 合． 納豆． ネギ・シイタケ・豆腐・ワカメの煮物．

2007 年 02 月 13 日 (火)

• 0740 起床． 朝飯． コーヒー． 査読原稿よみのつづき． 0955 自宅発． 雪． 1010 研究室着．
• 1030 より 研究室セミナー． 本日は矢澤さん & 宮田さんで， 私の下うけ記録だとそれぞれ 今月始め (サンプリングに関しては 昨年 10 月ごろ) とか 昨年末 あたりに苦闘してたハナシ． 昨日も書いたけど， どちらの解析も「次の一歩」がたいへんかも ……
• 1200 ごろ終了． 3/6 (火) - 3/7 (水) の ECOAS 仙台会議の案内くる …… おお， 私は何もやらなくていいんだ． こりゃらくちん． しかし来年度はアリ科研とこの ECOAS で何かやらんといかんおだろうね． いやはや．
• 昼飯． 昼飯後は雑用ばたばたと． 来年度から助教， という身分に変更されてしまうので， とりあえず 11 月 (??) あたりに 7 回ほど統計学授業， ということになるかも ……
• メイルなどによるデータ解析こんさるが断続的にわりこんでくる中， 例の査読原稿にとりくむ． はあ．
• …… 当方に読みまちがいによる誤解とかあってはマズいんで， かなり一生懸命よみなおしてみたんだけど， やっぱりこれってヘンだよね， というかかなりむちゃくちゃだよ …… と確信できたのはすでに日没時刻すぎ． うう．
• 却下の方向で査読報告文案かんがえつつ (これはかなり苦痛な作業) 撤退． 1810 研究室発． 1830 帰宅． 晩飯．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0800): 米麦 0.5 合． ネギ・シイタケ・豆腐・ワカメの煮物．
  • 昼 (1320): 研究室お茶部屋． イソップのバタール．
  • 晩 (2300): 米麦 0.7 合． ジャガイモ・ニンジン・タマネギ・牛肉のシチュー．

2007 年 02 月 14 日 (水)

• 0820 起床． 朝飯． コーヒー．
• いきなり ThinkPad X31 が起動不能状態におちいった． 私見ではこの製品は問題があり， しかも (IBM 製) AC アダプターと motherboard 電力まわりの両方が おかしいと思う． AC アダプターは 二年前 の発火事件以後 third party 製のものを使っているんだけど …… 今回は motherboard 側がおかしくなってしまったみたいだ．
• …… どうやっても回復しない， とあきらめがついたので 1130 自宅発． 曇． 1145 研究室着． 事後的に 4 時間の有給休暇を申請．
• 予備機 ThinkPad 240Z を およそ 400 日ぶり に起動． 「戦力化」作業 …… つまり モスボール 状態からの復旧．
• まあ， データについてはこまめに rsync で必要なぶんだけ高速バックアップ (つまり差分バックアップ) してるので身のまわりの計算機 2-3 台が壊れようが何も困らないんだけど．
• 1400 から 1800 ごろまで， 上田研究室 の石沢さん (卒研生!) あいてにアミノ酸濃度の階層ベイズモデルの詰めこみ教育． いろいろな R 作図プログラムをその場で作成しつつ．

• また忠良かつ堅実ただしさすがに動作の遅くなった予備機 ThinkPad 240Z の整備． そしてまたぐたぐだとメイル書き (データ解析こんさるとか) やって撤退． 2015 研究室発． 春一番の暴風雪． 2040 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • うんどう休養日， つづく
• [今日の食卓]
  • 朝 (0830): 米麦 0.6 合． 海藻スープ．
  • 昼 (1340): 研究室お茶部屋． イソップのバタールの残り． リンゴ．
  • 晩 (2200): 米麦 0.7 合． ミツバ・エノキダケ・豚肉の炒めもの． たちぽん．

2007 年 02 月 15 日 (木)

• 0810 起床． 朝飯． コーヒー． う， また朝から質問が …… 1020 自宅発． 晴． 1035 研究室着．
• 来年度の授業分担予定．

  生態系機能学総論(モデル中心）後期前半　月・木４時限
  
  この講義では、生態系の機能を解析する手法として現在用いられている数理モ 
  デル、統計モデルに関して基礎的な知識を概説することを目的とする。まず、 
  前半では動態を記述する数理モデルについて学び、生物の行動や形質の進化に 
  ついて議論する数理モデルについて学ぶ．また、後半では、統計モデルの基礎 
  となる考え方や最尤法について学び、その上で一般化線形モデルやより高度な 
  統計モデルについて学ぶ。
  
  (1-7 担当: 高田)
  1. 数理モデルの概説
  2. 動態編１（基礎モデル）
  3. 動態編２（推移行列モデル）
  4. 動態編３（群集モデル）
  5. 進化編１（最適戦略モデル）
  6. 進化編２（ESSモデル）
  7. 進化編３（共進化モデル）
  
  (8-14 担当: 久保)
  8.  生態学データ解析の統計モデリングとは?
  9.  さまざまな確率分布と最尤推定
  10. 一般化線形モデル (GLM) 1
  11. 一般化線形モデル (GLM) 2
  12. 検定とモデル選択
  13. 一般化線形混合モデル (GLMM)
  14. 階層ベイズモデル
  

• またメイル書きあれこれとか． うう．
• そして本日〆切の例の査読． もはや先方のやりくちは完全に理解できているので， あとはそれがいかにダメであるか， という説明だな．
• 査読報告書はけっきょく 550 words ぐらいで終わった． 現段階では細かい指摘とかしてもしょーがない， という判断から． 1600， web 上で査読報告提出終了 …… した瞬間に新しいデータ解析こんさるメイルがくる， とゆーのはどうなってるわけ?
• A 棟 7-8F 闇ネットの TCP/IP packet loss がヒドい， と長谷川さんから報告をうけて調査 …… といってもよくわからないので， 8F の階 hub 上にのってる闇ルーターの電源をいれなおしてみたところ， なぜかしら通信速度が平常に回復してしまった． ナゾ．
• Peter Dalgaard さんのすばらしい R 教科書 Introductory Statistics with R の翻訳版 Rによる医療統計学 が出版されて， なぜか私にも出版社経由で送っていただいた． ありがとうございます．

[R による医療統計学]

とても良い教科書． 「医療統計学」と書かれているけれど， どの分野の人がみても参考になる．

• さてさて， 本日のデータ解析こんさる業のつづき …… 自分でもよく知らないことについて， データ解析質問うけたんで， R 上で調べてみたんだが …… どうも plot(glm(...)) ってあまり意味ないんじゃないかな．
• たとえばポアソン回帰できるこういうデータがあったとしよう． データ数がこれだけあると glm() はほぼ正確な推定値を計算してくれる．

• で， fit <- glm(y ~ x, family = poisson) して plot(fit) するとどうなるか， というと下のように何だか意味不明ぎみというか奇怪な Q-Q plot 出すんだよね．

• Q-Q plot はもともと正規分布モデルのために作られた， 「あてはまりぐあい」を図示する方法だ． したがって， ポアソン回帰だのロジスティック回帰だのの推定結果を Q-Q plot するのがまちがい， ということになる．
• じゃあ， なんで glm() の出力オブジェクト (glm オブジェクト) でこんな図が書けるかといえば，

  > class(fit)
  [1] "glm" "lm" 
  

  とやればわかるよーに， glm オブジェクト は lm オブジェクト (正規分布の線形モデルの推定結果) を継承して定義されたものであり， glm(fit) とすると lm オブジェクト の plot() メソッドである plot.lm() が呼ばれている， というだけのことなんだよね …… いやはや， ヒトさわがせな．
• 1950 研究室発． 2005 帰宅． 晩飯．
• [今日の運動]
  • うんどう休養日， つづく ……
• [今日の食卓]
  • 朝 (0940): 米麦 0.6 合． ネギ・ブナシメジの味噌汁． 卵焼き．
  • 昼 (1440): 研究室お茶部屋． 米麦 0.5 合． レトルトパウチドカレー．
  • 晩 (2200): 米麦 0.7 合． ネギ・ブナシメジ・ニンニク・ホタテの酒蒸し． ニンニク茎・ピーマンの炒めもの．

2007 年 02 月 16 日 (金)

• 0840 起床． 生活周期ずれぎみ． 朝飯． コーヒー． 0930 自宅発． 曇． 0945 研究室着．
• 岩倉さん修論を読んで， ひたすら石狩浜アリ勉強．
• しかし「熱帯林研究にベイズ統計つかえそう」とサトりつつある 甲山さんと熱帯解析バナシを一時間ほど． なんだかまたカリマンタンだの Pasoh census データの 計算やりなおしをもくろんでおられるよーで ……
• アリ勉強のつづき．
• 1400 から二時間ほど， 東研 のアリ研究 M2 お二人， 岩倉さんと藤原君に石狩浜 エゾアカヤマアリ についていろいろと教えていただく …… うーむ， やはりというか， 詳細を知ってしまうと アリ worker 拉致 & 「敵地」投入実験のハナシは， なかなか素朴な解釈を拒絶する部分があるような気も． しかしアリの感情表現がかいま見える行動データはこれしかないんで， 上手に解析せんといかんね．
• 遺伝子データは例によって F_ST とかそういうハナシ …… 何か使いみちありそうだけど， これはあとまわしでも良い．
• そして体表面炭化水素， これまたやはりというべきか， じつは 25 次元ものデータがあれば， 簡単な線形判別モデルでもかなり細かく「判別」できてしまう (なにしろコロニー数は 7 にすぎないんで)， というのが真相のようで． neural network 使えばさらに改善されてしまうだろう． しかしながら， 行動データみるかぎり， おそらく現場のアリたちはそこまで厳密に区別しとらん可能性がある． とすると， このガスクロマトグラフィーデータを教育データとして使って 「だいたいのコロニー出自はわかるけど， それほど厳密ではない」 アリの神経系 (を模倣している統計モデル) を構築できるかどうか， といった問題になるのかもね ……
• アリばて …… その後， またメイル書きとか統計学勉強とか．
• 2120 研究室発． 2135 帰宅． 体重 68.2kg． 晩飯．
• [今日の運動]
  • ああ， 空腹状態で帰宅するとうんどう休養日になる確率が …… といいわけ
• [今日の食卓]
  • 朝 (0920): 米麦 0.6 合． ニンニク茎・ピーマンの炒めもの．
  • 昼 (1230): アリ勉強つづけつつ食パン．
  • 晩 (2300): 米麦 0.7 合． ネギ・キムチ・納豆いり卵焼き．

2007 年 02 月 17 日 (土)

• 0820 起床． 朝飯． コーヒー． 怠業しつつ洗濯． ひたすらひたすら洗濯． あと， デイパックその他の修理．
• 昼飯． 洗濯ばて怠業． 寝てしまう．
• で， 気づいたらすでに夜． なんたる一日．
• 明日は空港まで Heike Lischke さんを迎えにいかねば ……
• [今日の運動]
  • 嗚呼うんどう休養日
• [今日の食卓]
  • 朝 (0900): 東区のパン屋イソップの 焼きたてクロワッサン (大), 158 円．
  • 昼 (1350): バタール． ジャガイモ・ネギ・アンチョビーのクリームソース．
  • 晩 (2250): 米麦 0.7 合． チンゲンサイ・ショウガ・ネギ・カキの炒めもの．

2007 年 02 月 18 日 (日)

• 0820 起床． 朝飯． コーヒー． 1020 自宅発． 雪． 札幌駅すぐ南の Loft をちょっとうろうろしてから， 1140 快速エアポートに乗って札幌駅発． 1216 新千歳空港着．
• 1235 千歳着 (東京発) NH61 便， 1225 到着に変更? まあ， こういうのはえてして …… ほーら， 遅れた． 現実には 1231 着陸．
• で， 荷物うけとり所の前で待機していて迎撃すべき客人を発見． 先方はきょろきょろしておられるので， ガラスごしに手をふったらこちらを認識された． しかしなかなか荷物がでてこない．
• よーやく 1250 ごろスイスから (地環研主催の何やらアヤしげな) シンポジウム に招かれた Heike Lischke さん と合流できた． いやはや私は COE にも GLP にも何の関係もないんだけど， (甲山さんまわりのヒトで) 日曜日にこういう雑用を命じられるのはだいたいにおいて私， ということになっているので ……
• 昼飯はまだ不要， ということで 1304 また快速エアポートに乗って新千歳空港発． 札幌まで 36 分間， 1040 円． 車中では樹木・森林モデリングのハナシだの， それと統計モデルとの関連といったまにあっくな雑談が ……
• 1340 JR 札幌駅着． ただちに駅から 1 block 北のアスペンホテルにご案内． いっしょに昼飯でも， というハナシになったので Heike さんが部屋であれこれ 20 分間ほどしてるあいだに， 私はホテルロビーでなぜか統計学本とか読んでいたり．
• 昼飯， あまりたくさん食べたくない， ということで札幌駅地下のうどん屋へ． その後， 駅のすぐ南西にある紀伊國屋書店にいったり， 地環研 A 棟 8F までご案内してメイルを使ってもらったり (ホテルにはそういう端末がないようで)， またホテルまでお送りして研究室にもどる …… 1600 すぎ任務終了． ややぐったり．
• データ解析
• 1845 研究室発． ところどころ走りつつ帰宅． 1855 帰宅． 自宅前で待ってた郵パック配達のヒトから X31 マザーボードうけとる．
• うー， これはなかなかたいへんかも．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0920): 東区のパン屋イソップの 焼きたてクロワッサン (大), 158 円．
  • 昼 (1450): 札幌駅地下で野菜かきあげうどん． なんと! Lischke さんにおごっていただいた．
  • 晩 (2140): 米麦 0.7 合． ハクサイ・ネギ・シュンギク・シイタケ・豆腐・鶏の水炊き．

2007 年 02 月 19 日 (月)

• 0900 起床． うう． 朝飯． コーヒー． 昨日ちょっと着手しはじめた母子里原稿修正作業のつづき …… にはまりこんでしまった． 1200 自宅発． 晴． 1215 研究室着．
• 母子里林冠モデル の論文原稿ひたすら修正作業．
• 昼飯くって， ちょっとだけ COE シンポジウム に顔をだす． 吹雪の中の北大の学術交流会館 (conference hall) で開催されてるんだけど， なんというか宣伝してないので， 小講堂内には 20 人ぐらいがぱらぱらといるだけ． 劇ヘン なるホラのなれの果てがこれ?
• 苫小牧ボス・鍋嶋さん (甲山さんに命じられて東京から出張されたとか) の苫小牧樹木直径成長モデルのよく整理された講演を拝聴して撤退． 私が発表してたら， (モデリング方面とかに) すごく偏ったわけわからんモノになっていたであろう …… と考えつつまた吹雪の中を撤退．
• 1520 ごろから二時間ほど， また 東研 の岩倉さんから エゾアカヤマアリ について教えていただく． 本日はアリ行動と体表面炭化水素データをもってきていただいたので， それをあれこれと R などでながめながら石狩浜スーパーコロニーアリモデリングの検討．
  • アリ vs アリの行動連鎖数は 3 まで? 接触してもたいていはすぐに逃げるから (拉致されて敵地にほうりこまれたアリが)
  • away → home かみつきは一回 …… このヘンなデータどうしたもんかね
  • たしかに「一度かまれると次のアリにもかまれる」 傾向あるのかも …… しかしデータ上は「噛みつかれやすいアリと そうでないアリがいる」というモデルと区別できない
  • ガスクロマトグラフィーデータは 慎重に検討する必要あり …… この数値を相対値化してあつかってよいのか?
  • そもそもアリは「組成」を見ているのか? (たしかに neural network を使えば組成を「見る」 ことも可能ではあるけど)
• …… しかしアリモデリングはいったんわきにおき， 共著者のコメントなど参考にしつつ， 夜中まで母子里原稿なおしつづく．
• しかしながらまだ修正おわらず． 隅田さんから Gibbs 林冠説明に関する改善をいろいろ指摘されて， ですね． 2240 研究室発． 2300 帰宅． 母子里ばて．

[ThinkPad X31 基板]

ということで， 壊れ ThinkPad の修理になかなかとりくめない． Y! auction で私費で落札， 29980 円．

• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0950): 米麦 0.5 合． ネギ・豆腐の味噌汁． コマイ．
  • 昼 (1405): 研究室お茶部屋． 米麦 0.6 合． もずくスープ．
  • 晩 (1900): 研究室お茶部屋． イソップ角食． もずくスープ．

2007 年 02 月 20 日 (火)

• 0800 起床． 朝飯． コーヒー． 0900 自宅発． 晴 & 雪． 0920 研究室着．
• ネット雑用． web server まわりの設定を少々． PukiWiki の file attachment 機能を悪用した大量ファイル投入ロボット & ネットワーク形成という「荒し」があって， ですね． PukiWiki 使ってる研究機関， けっこうヤラれていますよ．
• とりあえず file attachment 一切禁止設定．
• で， その PukiWiki 荒しネットワークの中核的 server のひとつがわかったので， /etc/apache2/conf/httpd.conf で Deny from 設定． まあ， これは当方の不快感軽減ぐらいの意味しかないけど． こういうかんぢで error_log に延々と記録されていく:

  ...
  [Tue Feb 20 08:50:15 2007] [error] [client 85.17.103.104] client denied by
  server configuration: /home/ronbun/public_html/index.php, referer: http://
  www.wako.ac.jp/~michael/wiki/index.php?plugin=attach&pcmd=upload&page=1813
  ...
  

• そしてこれとは別件なんだけど， 敵国語版 Wiki を， という雑用圧力を二方向から受けつつあるわけだが …… ふーむ， UTF8 版 PukiWiki でどうにかなる問題でもないのかねえ (検討中) ……
• このあとの鍋嶋さん質問にそなえて 時系列データ解析の勉強． お， Wikipedia の カルマンフィルター 項がさらに充実しているな …… 再帰ベイズ推定との関係などが追加されている． EN 版 のほうは図つき． あ， Recursive Bayesian estimation の項が独立してたり， Particle filter の項なんかもあるな．
• 1100 すぎからお茶部屋で鍋嶋さんと年輪データ解析の話 …… とはいえ， 本日はとりあえず年輪の X 線写真 (!) を拝見して， 全体的な解析の方針について話を． その後， お昼御飯にお供させていただく．
• 昼飯から帰って， 一時間ほど母子里論文原稿なおし． あとで学術交流会館で原さん・隅田さんとこのへんについて ちょっと話しあう機会あり．
• 1440 学術交流会館へ． 昨日は COE しんぽとかいうやつだったけど， 本日は看板だけかけかえて GLP シンポジウム． 1450 から一時間ほど Heike Lischke さんの keynote speech． 森林動態モデル (gap → 地理的スケイル) の upscaling のやりかたについて． かかる面倒な内容を楽しくてぎわよくまとめた講演． 最後のほうは， adaptive gridding だのモーメント近似のハナシだのがでたり． とくに後者は全聴衆の中で私ぐらいしかわかってないんだろうなあ， と思ってみたり．
• そのあと， もと職場 (通称「独房群」) 地球フロンティアから甲山さんに呼ばれた佐藤さん SEIB-DGVM と伊藤さん (現在は環境研) Sim-Cycle ハナシ． すでに何度も聞いてるものだけど， 研究のハナシというのはオチがわかっていても 意外と何度でも楽しめるものだねえ．
• 最後は環境研の竹中さんのサトり (もしくはその逆) 気味なコメントで終了． 1800 ごろ研究室もどる． 矢澤さんと Pasoh 話を少々． 1830 ごろ研究室でて JR タワー 36F スカイバンケット で代替営業してるレストラン (ブッフェ) で， Heike Lischke さんを囲む晩餐会． 女王とその下僕 5 名 (甲山・高田・竹中・佐藤・久保) というかんぢ． 私以外は喋り好きなヒトたちばかりで， 楽しいひとときとなる．
• あ， そうだ． 10 分間ほどテイブル横で本職のマジシャンが手品やってくれた． ギミックというかシカケものも面白かったけど， トランプ 3 枚だけ + 指先わざだけで至近距離から見ている こちらの目をあっさりあざむけるのはスゴい， と素直に感心しました．
• 2100 ごろおひらき． Lischke さんを Aspen hotel までお送りしたあと， そのまま竹中さんの宿である北大前東横インロビーでモデリング雑談． キーワードは懐疑的・不可知論・中立モデル・べいぢあん …… うーん， ひねくれものの対話， というべきか．
• 竹中さんと話しているうちに， 生態学データ解析における多変量事前分布の使いみちがひとつわかった! 種差 (や個体差) のパラメーターが種間 (や個体間) 相互作用に依存してるかもしれない場合の推定だ …… いや， しかし分散共分散行列とかでは これを十分に表現できない場合もありそう?
• 2310 帰宅． 今日はかなりヒトの話を聞く日， になったような． 明朝は 0930 から二次試験監督説明会． 1030 からの研究室セミナーまでに終ってくれるかしらん?
• [今日の運動]
  • うんどう休養日， か
• [今日の食卓]
  • 朝 (0830): もりもとバタール． ニラ卵．
  • 昼 (1250): アスベスト工事避難のためクラーク会館 3F で営業中のきゃら亭． A ランチ 700 円， ホタテ & サーモンフライ．
  • 晩 (1900): JR タワー 36F スカイバンケットのブッフェ， 和洋中各種料理． 4000 円!