混合効果モデルに関するベイツの本の 12ページで、彼はモデルを次のように説明しています。
スクリーンショットの終わり近くで、彼は
相対共分散係数 に応じて、分散成分パラメータ、θ
正確に関係を説明せずに。我々が与えられていると言うどのように我々は導き出すだろう、Λのθ、それから?
関連するノートでは、これはベイツの説明が少し詳細に欠けていると私が思う多くの例の1つです。パラメータ推定の最適化プロセスとテスト統計の分布の証明を実際に通過するより良いテキストはありますか?
混合効果モデルに関するベイツの本の 12ページで、彼はモデルを次のように説明しています。
スクリーンショットの終わり近くで、彼は
相対共分散係数 に応じて、分散成分パラメータ、θ
正確に関係を説明せずに。我々が与えられていると言うどのように我々は導き出すだろう、Λのθ、それから?
関連するノートでは、これはベイツの説明が少し詳細に欠けていると私が思う多くの例の1つです。パラメータ推定の最適化プロセスとテスト統計の分布の証明を実際に通過するより良いテキストはありますか?
回答:
それは階層的な推論です。線形モデルには、bのコンポーネントである多数のパラメーターがあります。純粋な固定効果モデルでは、これらの推定値を取得するだけで済みます。代わりに、bの値自体が、thetaによってパラメーター化された共分散行列を持つ多変量正規分布から抽出されることを想像してください。これは簡単な例です。10か所の異なる場所で5つの異なる期間の動物数を確認するとします。count〜time + factor(location)のような線形モデル(ここではRトークを使用しています)が得られるため、(この場合は)すべての回帰に共通の勾配(それぞれに1つ)が得られます場所)が、場所ごとに異なる切片。パントして固定効果モデルと呼び、すべての切片を推定することができます。しかしながら、多数の可能な場所から選択された10の場所である場合、特定の場所は気にしないかもしれません。したがって、切片に共分散モデルを配置します。たとえば、切片を多変量正規であり、共通分散sigma2で独立であることを宣言します。その場合、sigma2は「シータ」パラメーターです。これは、各場所での切片の母集団を特徴付けるためです(これはランダム効果です)。