lmer研究7 プログラム改訂
一年前、公開したlmer()自動Akaike weight計算プログラムですが、lmerの仕様が変わったため使えなくなっていた。
原因を探ったら、意外と簡単に解決できたので、新しいプログラムを載せておく。
「model_average_lmer_new.r」をダウンロード
ランダム要因が一つのもの用
「model_average_lmer2_new.r」をダウンロード
ランダム要因が二つ以上のもの用
これを作ってもまた仕様がかわるかもなあ。
説明変数シラミつぶし関数をつかったものに変えた方がいいのかも。
| 固定リンク
「統計/R」カテゴリの記事
- PDF出力の日本語文字(2009.07.13)
- RODBC(2009.07.09)
- GUIでsource()(2009.03.07)
- factor()の順番(2009.02.18)
- lmer研究7 プログラム改訂(2009.02.17)
この記事へのコメントは終了しました。
コメント
lmer()は、Random effectが複数入る上に、切片と傾き両方のケースがありえるから厄介ですよね。プログラム上でRandom effectとそうでないものをどう見分けて総当りのモデルをつくるか。
あれだけ批判されていても、やはりlmer()のユーザーは多いです。間違ったGLMMの使い方を要求してくる教員やReviewerが多いせいでしょう。こうした状況は後何年ぐらい続くのか......
投稿: Totoron | 2009年2月18日 (水) 09時07分
lmer()はなんだかんだで使い勝手がいいんですよ。
複数のランダム要因の指定およびいろいろな分布形にも対応しているし。
ベイズは敷居が高いですから。
自分もまだベイズは使えないし。
投稿: martesorex | 2009年2月18日 (水) 16時44分