ぎょーむ日誌 2012-06-08

2012 年 06 月 08 日 (金)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0750 自宅発． 0820 研究室着．
• なぜかしら， なかなか進捗しない 統計モデリング入門 のちょくはんサイト準備 ……

[今日の店頭写真]

丸善丸の内本店 に陳列されている写真を送っていただきました． ありがとうございます．

• 昼飯．
• 午後はまず，雑用 ……
• オジロワシ繁殖の統計モデリングこんさる．
  1. データが少ないわりには， GIS な説明変数が過剰というありがちな状況 ……
  2. glm() + stepAIC() などやると， どう見ても無意味な (まったく解釈のできない) 説明変数を結合したようなモデルが AIC 最良となる (呪われな GIS 的説明変数!)
  3. とりあえず， データのとりかたをよく検討し， 「地域差」という効果があることがわかったので glmmML() を使うことに
  4. しかし現時点の最新の glmmML() …… これはまたしても以前のように stepAIC() できなくなってしまっている!
  5. しょーがないので 2003 年ごろ につくっていた自作の stepAIC.glmmML() 関数を復活させることに …… これ は少し修正しただけで使えた! (後記: この関数はかなりヘンなもので， stepAIC.glmmML("y ~ x1 + x2 + x3", cluster = d$id, data = d, ...) といった引数指定をします． 数式が formula ではなく文字列!!)
  6. しかし呪われな GIS 的説明変数の個数はなかなか減らない
  7. 呪われな GIS 的説明変数どうしの正負の相関を pairs() で図示して確認
  8. 論文の method section における「これこれこういう理由で」 という部分の正当化を考え， 「少数個の説明変数を使ったモデルたち」 という範囲で AIC によるモデル選択， これにて一件落着
  というかんじで， とりあえず何とかなった． しかし， また後日， これまた呪われぎみになりがちな「気象要因」データを 「どう捨てるか」について研究しなければならぬことに．
• 16 時 30 分ごろから 1.5 時間ほど， 北大祭の展示雑用 …… ぐったり．
• 1900 研究室発． 買いもの． 1930 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1220): 研究室． 米 0.7 合． 興部ハーブソーセイジ． アオノリ天． キュウリのピクルス． ミニトマト． ミズナ．
  • 晩 (2000): 米 0.7 合． 冷凍保存しといたチキン野菜カレー．