ぎょーむ日誌 2010-04-26

2010 年 04 月 26 日 (月)

• 0730 起床． ねむい． 朝飯． コーヒー． 0830 自宅発． 晴． 0840 研究室着．
• メイルかき & 研究室セミナーでの自己紹介の準備 ……
• 1030 より 研究室セミナー， 今日も 7-8 人の 5 分間研究自己紹介．
• 自分の自己紹介． やはりウケたのはここでした (アフリカンな大学院生たちも笑ってくれた!)． こういうヘンな画像をネット上でさがすのに時間かかりました (いちおう出典サイト URL も併記)．

• 蛇足ながら， "data cooking" にはデータ捏造といった悪い意味で使われることもあるけど， ここではとりあえずにゅーとらるな表現ということで． マトメはこんなかんぢで．

• そのあとは保要さんのサロベツなモウセンゴケ研究計画説明．
• 昼飯．
• 甲山さん指令で新しい液晶プロジェクターを手配することに． といっても， 北大生協の外売部 (地環研担当は今年度から八木さんにかわった) に依頼しただけなんだけど ……

  研究室で購入する液晶プロジェクターを探しています．
  条件としては
  
  ・十分に明るい
   - 40 人程度入れる教室で使う
   - 投影スクリーンは 3m 四方 ぐらい
   - スクリーンまでの距離は 5-7 m ぐらい
  ・静音性ができるだけ良いもの
  
  その他の条件はとくにありません．
  

  このあといただいたメイルによると， 北大生協 外売店専用WEB受注サイト なるものを作ったそうで …… うーむ， これはまだ発展途上かな?
• で， 他のことは放置して岩波本作業にとりくんでみる ……
• なかなかススまない …… 生協のヒトがやってきて， セミナー室であれこれ相談． じつは上記のスクリーンサイズはかなりいいかげんで， 実測は幅 1.7 m・高さ 1.3 m だった． 甲山さんがやってきて， 白い壁に直接投影するプランもあるかも， と提案． 4000 ルーメンぐらいの機器は 20 万円ぐらいするので， 後日またここでテストしましょう， という生協のヒトの提案で終了．
• 本日の岩波本作業， 7.8 節「GLM のベイズモデル化」 の図の準備とか．

[二項分布の確率 q の事前分布]

logit(q) = a としたときに， (A) a の事前分布が平均ゼロ SD 100 ぐらいのいわゆる無情報事前分布としたときに， (B) q = logistic(a) がどうなるかを図示したもの． この (B) の図示では， もとの logit な世界のスケイルの確率密度を正確に変換するあたりで， ちょっとアタマを使う必要がある．

• ところで， 先日かいたつじつまあわせの部分では， 確率 q の事前分布を一様分布とゆーことにしてハナシを 強行していったわけだが ……

• こういう二項分布の確率 q では以下の三つの 「無情報」事前分布が設定されることが多い:
  • 一様分布 U(0, 1) (== ベータ分布 Beta(1, 1))
  • Jeffery の事前分布: Beta(0.5, 0.5) ⇔ 一次元のランダムウォーク → arcsin 関数: 0.5 π arcsin(q^0.5) (これは岩波「科学」 2008 年 4 月号の p.402)
  • Prob(logit(q)) = 一定とすると， q の事前分布は improper な Beta(0, 0) となる
  といったことが Gelam et al. (2004) Bayesian Data Analysis の p.63 あたりにでてる …… で，じつよー的には Beta(0, 0) から Beta(1, 1) までのあいだの どれでもいーんじゃない，といった記述なんかも．
• このあたり， 林さんのメモ も参照． ここで粕谷さんがコメントしてる Lee (2004) Bayesian Statistics: an introduction (3rd ed.) だと p.85 あたり．
• やはり Jeffery な事前分布の図も追加するかな …… まあ， これは明日にしておこう． 1920 研究室発． 買いもの． 1940 帰宅． 晩飯の準備． 晩飯．
• ばててたのか， 晩飯くったらいきなり寝てしまった …… なぜか子供のころに住んでたあたりの道路がどこも断崖絶壁になっていて， 昔はするする登っていたんだけど， 今はカラダが重くてなかなか登攀できない …… といった悪夢にうなされて真夜中めざめてしまった． うーむ ……
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0800): 米麦 0.6 合． ワカメ・豆腐の味噌汁．
  • 昼 (1300): 米麦 0.8 合． キャベツ・タマネギ・カツオブシのサラダ．
  • 晩 (2030): 米麦 0.8 合． タマネギ・ネギ・豚肉炒飯．