混合モデル

混合モデル（こんごうモデル、）とは、固定効果（fixed effect）と変量効果（random effect）を共に含む（ゆえに混合効果と呼ばれる）統計学的モデルであり、医学・生物学・社会科学等の広い領域に用いられる。特に縦断研究においてある項目を繰り返し観察する等で有用である。欠測データの取り扱いに優れ、混合効果モデルは多くの場合、反復測定分散分析等の伝統的なアプローチよりも望ましい。
歴史と現状.
1918年、ロナルド・フィッシャーが関連する項目間の特性を変量効果モデルとして導き出した事に由来する。1950年代、チャールズ・ヘンダーソンが固定効果モデルである最良線形不偏推定量（BLUE）および変量効果モデルである（BLUP）を規定した。
その後、混合モデルは最尤推定量、非線形混合効果モデル、欠測のあるモデル、混合効果モデルのベイズ推定量の計算等に用いられるようになった。混合モデルは各測定点での値が相互に影響し合うケースに用いられ、現在ではヒトに対する臨床試験や動物実験で汎用されているほか、。
定義.
行列を用いて以下の様に記述する。
ここで、
推定.
formula_2 と formula_5 の結合密度関数は次の様に書ける: formula_18
formula_19 と formula_20 および formula_21 には正規分布を仮定し、formula_4 と formula_5 の同時密度関数（結合密度関数とも）を最大化すると、ヘンダーソンの“mixed model equations (MME)”が得られる。
このMMEを解く時、 formula_25 と formula_26 はそれぞれ、formula_4 と formula_5 の最良線形不偏推定量（BLUE）と最良線形不偏予測量（BLUP）である。これは、目的変数の条件付き分散が単位行列のスカラー倍にならない場合のガウス＝マルコフの定理の解である。条件付き分散が既知である時、逆分散の加重最小二乗推定値はBLUEであるが、条件付き分散が既知であることは稀である。従ってMMEを解く時は、分散と加重推定値を同時推定する必要がある。
この様な混合モデルに適用する方法の一つとして、EMアルゴリズムがある。EMアルゴリズムにおいては分散成分が結合尤度における未観測のとして扱われる。現在は、R言語（「nlme」ライブラリの「lme」関数）や（「proc mixed」プロシジャ）に実装されている。混合モデル式の解法として、誤差が正規分布する場合は最尤推定法を用いる。