ぎょーむ日誌 2006-07-18

2006 年 07 月 18 日 (火)

• 0750 起床． 朝飯． コーヒー．
• 今朝もまた (やや現実逃避的と評価せざるをえない) 光計算問題にとらっぷされてしまった． 昨日すでに 「計算量を激減できる画期的な方法はこの天地のあいだには (あったとしても) 見つからない」 とゆーサトリの境地に到達した …… というかすでに数日前に， (光遮断オブジェクト数) × (見る方向) 個ぶんの計算はおそらく回避不可能で， 最悪の場合には (光遮断オブジェクト数) × (光遮断オブジェクト数) × (見る方向) ぐらいのケタになる， などと自分で書いてるじゃない． 極端にいえば， この範囲こそが 「人間の努力の限界」 ⇔ 「愚かしさ」について規定しているのでは， と思えてくる．
• で， 昨日から水平面を分割することばかり考えてたわけだが， 垂直面や天球のうまい分割ができれば何かよいことないのか， と思ったりしたわけだが …… うーむ． 天球の区画化は案外重要かも． しかしよくわからん． 0955 自宅発． 曇． 1010 研究室着．
• ネット雑用いくつか．
• そしてまた中国奥地に anova.glm() の結果を送らねばならなかったり．
• RjpWiki トップペイジに朝倉書店の新刊本のハナシがでてるな．
  • 数理統計学ハンドブック : Hogg の Introduction to Mathematical Statistics の翻訳
  • ビジネスへの統計モデルアプローチ : 「びぢねす」とあるけど， 目次をみると意外とおもしろそう ……

    3. 芝居をプロデュースする−統計モデルのプランニング
      3.1 芝居の目的と配役の検討−勧善懲悪ドラマと推理ドラマ

    おそらく「勧善懲悪ドラマ」とは Neyman-Pearson わくぐみでの (つまりふつーの) 「検定」のことなんだろうなぁ …… 最初から「とうぜん到達するべき結末」 が決まってる点において「勧善懲悪ドラマ」だよね
  (後記: 実際に同書を買ってみたら
  1. 勧善懲悪ドラマ: 入力変量と出力変量の因果構造を あきらかにする構造モデル
  2. 推理ドラマ: 出力の背後にある「概念」を探り， その level を推定する測定モデル
  という何だかよくわかぬものだった)
• つまり「ゆーい差決戦主義」が成立するためには， Neyman-Pearson lemma が必要不可欠である， と． そのわりには主義者どもは最強力検定たる尤度比検定とかやらないわけだが． 自分でもよくど忘れするんで注記しとくと， 尤度比検定が最強力なのは単純仮説の比較において …… 複合仮説の場合はどう考えたらよいものやら， 私はまだよくわかってない．
• ホーエル訳本「入門数理統計学 (4th)」 の尤度比検定の説明を (自分につごうよい方向で) 解釈してみると，
  1. 対立仮説が複合仮説である場合には， 複合仮説内のうちもっとも尤度が大きい仮説で その複合仮説全体を代表させる
  2. これで対立仮説は (最尤推定値という一点の値をつかった) 単純仮説と帰無仮説 (これはふつーは単純仮説) の尤度比検定は， 単純仮説どうしの比較になっているので Neyman-Peason が使える
  ってことかしらん? つまり最尤推定値を使った尤度比検定ってのは， 1. のぽりしーを認めるなら， つねに最強力検定であると考えてよい， ということだろうか．
• で， Neyman-Pearson で思い出してしまったので， また 山口さん ところにコメントなど書いてみる． うーむ， あまり役にたたんこと書いてるような気がする．
• するとどんどんこちらの深みにはまりこんでしまった． 岩波の統計科学のフロンティア 「モデル選択」 の 下平-長谷川検定 のところをながめる …… これすごくややこしいんだけど (AIC + ブートストラップ法 + 多重比較)， 「SH 検定って何するものなの?」 と疑問もつヒトたちをケムにまける程度の理解はできた．
• たとえば AIC のばらつきに関して，
  • AIC のばらつきを評価する bootstrap 法において， 「AIC そのもの」 のばらつきを見ても意味がなく， 「モデル間の AIC の差」 が bootstrap データでどう変わるかをみる
  • しかし単純な bootstrap 法ではだめで， bias 補正を考慮した下平-長谷川検定でやっているような multiscale bootstrap 法が必要である
  • なぜ multiscale bootstrap が必要かといえば， bias 補正のためである …… このあたりどういう理由で bias が発生して， それをどう補正すればいいのか， そのあたり私にはまだよくわかっていない
  などなどは知っておいてソンではありませんよ， とか何とか．
• ということで， また山口さんのところでコメント． うーむ， この午後は統計学の基礎 (すなわちめんどうな問題) に関する表面的な理解を進展させるひととき， になりつつあるなあ．
• 環境生物科学部門 における今年度の教育経費なる数字が．

  今年度の教育経費は、学院拠出分１３千円をマイナスして
  修士：７６千円
  博士：１１４千円
  です。
  計画的な費消をお願いいたします。
  

• 1815 研究室発． 1830 帰宅． 体重 74.6kg． 昼飯くわないぐらいではなかなか体重おちないね． この時刻なので NHK ラジオ夕刊 が放送されており， metabolic syndrome についてあれこれと． ウェスト 85cm 以上で「アウト」ですか．
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0830): 米麦 0.7 合． ネギ・コンニャク・豆腐の味噌汁．
  • 昼 : 食ってない
  • 晩 (1900): 米麦 0.8 合． ニラ卵炒飯． ネギ・ワカメ・豆腐の味噌汁．