const month dbh.log dbh.log.sq ...
1 TRUE TRUE TRUE TRUE
2 TRUE FALSE TRUE TRUE
3 TRUE TRUE FALSE TRUE
4 TRUE FALSE FALSE TRUE
5 TRUE TRUE TRUE FALSE
6 TRUE FALSE TRUE FALSE
7 TRUE TRUE FALSE FALSE
8 TRUE FALSE FALSE FALSE
...
……
こうなってしまうモデル選択の問題はおいといて，
full model 推定やらんといかんな．
full model からしてすでに面倒そう．
これで計算おそかったら
R
で計算させるのやめよかな．
exp()，
忘れじの log()」
で 90 分つぶせそうだな．
lme4
ライブラリに新型の一般化線形混合モデル推定関数
GLMM()
がふくまれてるらしい．
推定計算法は罰則つき尤度法 (penralized quasi-likelihood; PQL)
でいまいち感あるけど
……
おお，
R-help ml とか検索すると，
またしても Ripley 先生が
「glmmPQL(MASS)
のどこがイカンのじゃ」
(むろん表現はこれよりずっと英国的)
とけんか腰ですぞ．
sapply()
のたぐいが．
これいれると計算すごく遅くなるのに，
ついつい sapply()
してしまう，
と．
vector の計算はほどかずに vector のままで，
というふうに手直しする．