ぎょーむ日誌 2010-01-23
2010 年 01 月 23 日 (土)
-
0830 起床.
よくねた
……
朝飯.
コーヒー.
洗濯.
怠業
……
ばかりもしてられないんだよね.
-
せとぎわ感あるメイルがきてるので,
データ解析こんさるにとりくむ.
ヒヨコ実験データ解析の卒業研究
……
はい,
理学部からです.
-
卒研生に階層ベイズモデルってのもアレなので,
R だけで簡単にできるところを
……
交互作用がちょっと独特なので,
x1 + x2 + x3 + x1:x3 + x2:x3
みたいに *
演算子を使わずにぽちぽちと指定していく
……
-
stepAIC()
なモデル選択やらせると,
この関数に独特なルール
「x1:x3
項が残ると,
x1
と x3
は必ず残る」
が適用されてみとおしの悪い結果になった
……
とゆーかみとおし悪くなった責任は,
stepAIC()
ではなく,
ありがちな線形予測子を指定している当方にあり,
だな.
-
ここは工夫が必要とされるところだ.
この実験データでは
x3
(って仮名ですよ)
は他とは異なる独特なものであり,
これを as.numeric()
などして「因子型変数の交互作用と同等な数値型変数」
の列をどんどん追加していき
(じつは x3
は二値型なので簡単にできる)
……
これでパラメーター推定 &
モデル選択するとすっきりとわかりやすく,
ついでに AIC も低くなるモデルが得られた!
報告かき.
メイル送信.
-
……
いやー,
これってすごくおもしろいヒヨコちゃんうろうろ実験の結果ではなかろーか,
などと考えつつ遅い昼飯.
-
1450 自宅発.
曇.
1510 研究室着.
R 棟の廊下が暖かいような気がしたので,
「今日は暖房運転してるのかな?」
と思ったけどじつはいつもの土日と同じく最小出力運転だった.
防寒態勢をとる.
-
さっさと 2/9 の公開講座の準備せんといかんのだが
……
なかなかススまず
……
-
室内もかなり寒くなってきたので,
1820 研究室発.
買いもの.
1900 帰宅.
晩飯.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0900):
米麦 0.6 合.
ナメコ・ネギの味噌汁.
- 昼 (1400):
パン.
ホタテスープ.
- 晩 (2040):
米麦 0.8 合.
コマツナ・ブナシメジの味噌汁.
ジャガイモ・ニンジン・タマネギ・しらたき・牛肉の煮物.