ネストされていないモデルとAICの比較


19

GLMMにしなければならないと言う

mod1 <- glmer(y ~ x + A + (1|g), data = dat)
mod2 <- glmer(y ~ x + B + (1|g), data = dat)

これらのモデルは、通常の意味ではネストされていません。

a <- glmer(y ~ x + A + (1|g),     data = dat)
b <- glmer(y ~ x + A + B + (1|g), data = dat)

私たちは行うことはできませんanova(mod1, mod2)我々はと同じようにanova(a ,b)

代わりにAICを使用して、どちらが最適なモデルであるかを指定できますか?

回答:


24

AICは、ネストされていないモデルに適用できます。実際、これはAICについて最も拡張された神話(誤解?)の1つです。見る:

注意する必要があることの1つは、すべての正規化定数を含めることです。これらは、異なる(ネストされていない)モデルによって異なるためです。

こちらもご覧ください:

GLMMのコンテキストでは、この種のモデルを比較するためのAICの信頼性がより繊細な問題になります(@BenBolkerも参照)。AICの他のバージョンについては、次のペーパーで説明および比較しています。


2
異なるモデルを比較しようとすると、限界対条件付きAICの区別が最も重要であることに注意してくださいランダム効果の彼らのセットでの
ベンBolker

@ChandelierとBen Bolkerは、両方の答えをありがとう。このようにAICを使用するための議論について、どちらかがより正式なリファレンスを持っていますか?
user1322296 14

2
@ user1322296ルートに行くことをお勧めします。これは赤池の論文です。AICは、モデルと「真のモデル」との相違の推定量として取得されます。したがって、ネストは想定されず、いくつかの規則性条件のみが想定されます。
シャンデリア

lm1 = x〜A + BのAICを比較することは有効ですか Cとlm2 = x〜D Cか?おかげで
crazjo

AICの使用が適切でないネストされていないモデルがあるようです。以下に2つの例を示します。12です。ネストされていないモデルの選択が機能する条件をいくつか教えてください。
カール

10

参考のため、反論:ブライアン・リプリーは「モデルの大きなクラスの中から選択する」 pp。6-7

重要な仮定 ...モデルはネストされています(脚注:赤池(1973)の再版の615ページの下部を参照)。– AICが広く使用されていない場合

関連する文節(赤池の別の再版のp。204)は、「統計モデルの同定の問題はの選択の問題として定式化されることが多い」というフレーズから始まると思うf(x|kθ

Ripley、BD2004。「モデルの大きなクラスの中からの選択」。N。Adams、M。Crowder、D。J Hand、およびD. Stephensが編集した統計の方法とモデル、155〜70。イギリス、ロンドン:Imperial College Press。

赤池H.(1973)情報理論と最尤法の拡張。で情報理論に関する第2回国際シンポジウム(編BNペトロフとF.なCaSki)、頁267から281まで、ブダペスト。アカデミアカイドー。Breakthroughs in Statistics 、eds Kotz、S。に転載 &Johnson、NL(1992)、volume I、pp。599–624。ニューヨーク:スプリンガー。


3

赤池は、AICはネストされていないモデルを比較するための有用なツールであると考えたようです。

「AICに関する重要な観察事項の1つは、真のモデル[f(x |kθ)]への特定の参照なしで定義されていることです。したがって、任意の有限数のパラメトリックモデルについて、 [f(x |kθ)]これは、少なくとも原則として、AICがネストされていないモデルの比較、つまり、従来の対数尤度比検定が適用できない状況の比較に役立つことを示唆しています。

(赤池1985、ページ399)

赤池、広t。「予測とエントロピー」。赤池弘tの論文 スプリンガー、ニューヨーク、ニューヨーク、1985年。387-410。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.