ぎょーむ日誌 2010-07-25

2010 年 07 月 25 日 (日)

• 0820 起床． よくねてしまった …… 朝飯の準備． 朝飯．

[ネパールの青いケシ]

結婚おいわいということで， 卒業した院生から送っていただいた． 夏休みにネパールに行ってきたそうで …… どうもありがとうございます． 元の画像ファイルは こちら (0.7 MB)．

• 今日は日射が強いので (気温は低いけど)， 北大構内走はあとまわしにして， 洗濯． 怠業． 1030 自宅発． 晴． 1050 研究室着． コーヒー．
• 昨日の いまさらながら 対数尤度についていろいろ「発見」， まあ「自分がきちんと理解してなかったことを発見」 してしまったので …… いまさらながらに 「情報量統計学」(1983) と「情報量規準」(2004) をみなおしてみる．
• 「情報量統計学」には上下それぞれの「ずれ」 に関してカイ二乗分布と書いてある． なぜやたらとカイ二乗分布なんぞがもちだされるかというと， 最尤推定量 (estimator) が (ある条件下では) 漸近正規性があるから．
• 「情報量規準」ではカイ二乗分布とか登場しないんだよね． もっと一般化されてるのか? …… うーむ， そのかわりにというか， 多変量中心極限定理をもちだしてるな． いや， これは最尤推定量の漸近正規性を示すためか．
• 現代的には， カイ二乗分布はいらない? まあ， 「説明のひとつとして」 カイ二乗分布をもちだしてもいいんだろうけど．
• 昼飯．
• まあ， これ以上うだうだと考えてる時間もないし， えいやっと書いてしまうしかないな …… 岩波本作業のつづき．
• まずは作図プログラミング． 岩波本， 絵本化しつつあるな ……
• これが第 4 章 AIC モデル選択章の最後のほうの図． こんなンでいいのかな?

• えーと， この図からわかることは， 無意味なパラメーター追加によって AIC は平均 1 増加 (-2 をかけるので)． AIC が変化しないなら， 推定用データへの deviance 減少は 2 ぐらい． deviance 減少 1 との確率を近似計算すると:

  > pchisq(1:2, 1)
  [1] 0.68269 0.84270
  

• しかし， 赤池さんは紙とエンピツだけで， よくもまあこういうことを思いつくもんですなぁ …… たぶん， ここに到達する経路はぜんぜん別のものかもしれないけれど．
• 統計学の歴史上の偉人たち， 皆さんアタマおかしいとしか思えない． (おおくの場合) 計算機実験が不可能な状況のもとで， 想像困難な世界において正しい方向に空想をハシらせるものだと， いまさらながらに感心させられますよ． 人間の想像力はすごい．
• まあ， たしかにネストした統計モデル間の比較なら ……

• 1900 研究室発． 1910 帰宅． 1920 自宅発北大構内走． 2000 帰宅． 体重 71.4 kg． 晩飯の準備． 晩飯．
• おお， 南半球旅行中の伊庭さんにいろいろご心配をおかけしてるみたいだ ……

• うーむ， 「3 本の放物線」 というと， ブートストラップなハナシですかね ……? そっち方面のハナシはしない予定なのですが． ともあれ， 明日， 小西・北川本を確認しなくては．
• 夜中にまゆみさんが道東から帰ってきた．
• [今日の運動]
  • 北大構内走 40 分間．
• [今日の食卓]
  • 朝 (0910): 米 0.6 合． キャベツ・ニンジン・タマネギ・鶏ササミのスープ． トマト・卵の炒めもの．
  • 昼 (1300): 研究室． 米 0.8 合． レトルトパウチドのカレー． トマト．
  • 晩 (2130): 米 0.8 合． トマト・ネギ・卵炒飯． キャベツ・ニンジン・タマネギ・鶏ササミのスープ． 今日はトマトばかりくってたな．