ぎょーむ日誌 2010-07-16

2010 年 07 月 16 日 (金)

• 0720 起床． 朝飯の準備． 朝飯． コーヒー． 0900 自宅発． 晴． 0920 研究室着．
• あれこれののちに， 岩波本作業のつづき． さてさて， 今日から第 4 章のモデル選択・検定章の書きなおしなのだが …… どーにもこーにもモデル選択まわりの説明がひどすぎる． というか， 説明にもなってないような．
• いやべつに AIC の導出なんぞを岩波本で示す必要はないんですよ． そんなのは， 小西・北川「情報量規準」本とか読めばいいんだから．
• とはいえ， 逆に， ありがちな「2 かけるパラメーター数はぺなるちーで」 といった 「なんちゃって美しきぱーしもにー実現する厳罰主義」 みたいな安易な説明なんかも通俗宗教的なかんぢで (職業的うそつきとしては) おもしろくないな， と． 今までさんざん使ってきましたが．
• そこで， とりあえず， 「とにかく何でも具体的な数値例・図解」 な岩波本のぽりしー (いつのまにやら?) の旗をかかげつつ ……
• 構成変更ぷらんを検討してみる． これまでの第 4 章を分断して，
  • 新第 4 章: 旧第 4 章アタマのポアソン回帰の AIC モデル選択はそのまま， 図解・平均対数尤度 (これは従来の教科書ではなかったもの)， 旧第 4 章末尾につけてた FAQ モデル選択 のちょー圧縮版 (つづきは web 版みろ的な)
  • 新第 5 章: 旧第 4 章の尤度比検定パラメトリックブートストラップ版そのまま， よけーな記述は削除する
  といった構成を検討中． 新第 4, 5 章どちらも短い． 現在の構成は章の長さのばらつき大きいのだけど， 主題によって章の長さがかわるのはとうぜん， などとひらきなおっている今日このごろ．
• さてさて， 平均対数尤度 & 最大化対数尤度のバイアス補正 なる概念は下のよーな図で説明されることが多く， これがわかりやすいような気もするのだが ……

• えーと上の図は， 川崎能典さん の AIC 導出の説明 スライド の一部を利用させてもらいました．
• R をひねくって数値例をいろいろ生成してみると …… 実際のところは， たとえばこんな描像だ．

• 自分のど忘れ用にメモしとくと， 青曲線が平均対数尤度， 赤点が最尤推定値・最大化対数尤度， 赤破線がその標本平均， グレイ曲線がある推定における対数尤度関数．
• 1925 研究室発． 買いもの． 2000 帰宅． 晩飯の準備． 晩飯．

[100 均の工具]

いつものごとく北 12 生協で晩飯の買いもの後に帰宅中， 「いまや北大周辺で (自転車修理に使える) 工具うってる店はないな」 (昼飯後にちょっとヨドバシにいったんだけど， 工具はまったく置いてなかった) …… 「そうだ北 12 地下鉄駅そばの 100 円均一店はけっこう工具あった!」 と寄り道してみた． たしかにいろいろあるんだけど大口径棒レンチはなかった．

• [今日の運動]
  • うんどう休養日， つづく ……
• [今日の食卓]
  • 朝 (0740): 米 0.6 合． 油揚・ワカメの味噌汁． コマイ．
  • 昼 (1300): 研究室． シリアル．
  • 晩 (2200): 米 0.8 合． ホウレンソウ． タマネギ・ズッキーニ・ニンジン・ピーマン・トマトの煮物． ニラ・ショウガ・卵・豚肉の炒めもの …… らでぃっしゅぼーやで野菜がたくさんとどいたので， 野菜な晩飯に．