ぎょーむ日誌 2007-03-15

2007 年 03 月 15 日 (木)

• 0900 起床． うう， 9 時間も寝てしまった …… しかしまだアタマがややへろへろしている． コーヒー． 朝飯． 1020 自宅発． 晴． 1035 研究室着．
• アタマがへろへろしてる理由: (昨日の作業である) LaTeX ポスター作りはある意味ラクでいいんだけど， これってアタマの中でレイアウトの座標計算やりながら 作成することになるんで， このあたりが意外と (自分でも気づかぬうちに) 負担になっているんだろうな ……
• 窓の外はときどき吹雪． 三月になってもこんなに降るとはなかなかのモノだな．
• 午前中は何だかメイル書いたり， 屋久島 web page や自由集会サイトまわりでふらふらしていたり …… いまやありがたいことにネット上をさまようだけで 愛媛大学→「大街道」周辺の地理がアタマに入ったりするんだけど， 自由集会後のうちあげ会場はいつものことながら何とも決めがたい． 人数が把握できないから， という理由もあるんだけど．
• 気分転換にちょっと昼飯調達にでる． 研究室にもどって昼飯．
• 1400 から 2 時間弱ほど石沢さんと豊平川アミノ酸の 階層ベイズモデルのハナシ． じつはこれは (来年度あれこれやらねばならぬ) 群集生態学ハナシのモデリングとも関連ありそうだ …… といった統計モデリングおたくな話はともかく， このアミノ酸うんぬんとサケ母川回帰のハナシがどう関連するか いまいちよくわからない． おそらく石狩川・豊平川溯上してるサケが 分岐点あたり で「えーと，豊平川は右か」とか判断するのに使うんだろう (??)けど， 札幌市豊平川さけ科学館 の豊平川のサケ標識放流調査ってのを見ても 「どれぐらい母川を識別してるか」 はわからない． 「まちがって左の石狩川とかに行ってしまいました」 データはどこにもないんで．
• 粕谷さんと 自由集会 メイルやりとり …… 粕谷さんの一般講演 (ポスター) は「主成分分析と検定」なるスゴそうなタイトルとわかったので， ネット上の要旨 を読んでしまった …… PCA と統計学的検定を連関させんのは意外と難しいような気もする． そりゃ， bootstrap とか使えばどういう場合でも何らかの P 値みたいなものは計算できるでしょうけどね ……
• 経験ベイズ法の例題としては， 縮小推定量 (Stein 推定量) のハナシとかわかりやすいかも …… たとえば， 「100 人のプロ野球選手それぞれ 5 打席ぶんのデータから 各打者の『真の打率』 (?) を推定しなさい」， みたいな． 打者ごとに割りざん推定してしまうと 「0 割打者」とか「10 割打者」とか ぞろぞろ出るんだけど， 階層ベイズモデル化 (打者 100 人の「個体差」考慮しつつ， 全員のデータから「打率の事前分布」を推定しつつ) → 経験ベイズ法で推定計算 (これは最尤推定) するともうちょっとマシな推定結果になる (5 打数 5 安打でも「10 割打者」とはならない)， といったハナシですね．
• 中間さん からのメイルを拝読してるうちに， 今さらながら気づいたんだけど …… 2 年ちょい前 のぎょーむ日誌で R の expression() で日本語文字があつかえない， とかぶちぶち文句を言ってたんだけど， なかまさんたちの 尽力 によってその後の version の R では ちゃんとあつかえるようになっている， と今ごろになって気づいた! device は pdf() や postscript() でも問題ナシ． R 開発者の皆さん， ありがとう ……

  x <- rnorm(100)
  y <- rnorm(100)
  plot(x, y, xlab = expression(Delta * beta["control"]^italic(Z) * " の値"))
  

• 自由集会準備は具体的にはあまりススまぬまま撤退． 1910 研究室発． 1930 帰宅． 体重 68.6kg． 晩飯の準備． 晩飯．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 10 × 3 回． スクワット 100 回．
• [今日の食卓]
  • 朝 (0930): 米麦 0.5 合． ワカメスープ．
  • 昼 (1320): 研究室． 北大生協にぎり飯．
  • 晩 (2220): 米麦 0.7 合． タマネギ・ブナシメジ・ホタテのカレー． ダイコンのゴマサラダ．