階層線形モデル(HLM)の構造と予測値計算
階層線形モデル(Hierarchical Linear Model, HLM)は、複数レベルの構造を持つデータを分析するための統計モデルです。教育、医学、組織研究など様々な分野で使用されます。
2レベル階層線形モデルの基本構造
レベル1(個体レベル):
$y_{ij} = \beta_{0j} + \beta_{1j}x_{ij} + r_{ij}$
レベル2(グループレベル):
$\beta_{0j} = \gamma_{00} + \gamma_{01}W_j + u_{0j}$
$\beta_{1j} = \gamma_{10} + \gamma_{11}W_j + u_{1j}$
Step 1: モデル成分の説明
レベル1モデルの構成要素:
- $y_{ij}$:グループ$j$の個体$i$の応答変数
- $\beta_{0j}$:グループ$j$の切片
- $\beta_{1j}$:グループ$j$の傾き
- $x_{ij}$:個体レベルの説明変数
- $r_{ij}$:個体レベルの誤差項
レベル2モデルの構成要素:
- $\gamma_{00}$:全体の切片
- $\gamma_{01}$:グループレベル変数$W_j$の切片への効果
- $\gamma_{10}$:全体の傾き
- $\gamma_{11}$:グループレベル変数$W_j$の傾きへの効果(本問では0)
- $u_{0j}, u_{1j}$:グループレベルのランダム効果
Step 2: 問題設定の確認
与えられた条件:
- $\gamma_{00} = 10$(全体の切片)
- $\gamma_{01} = 2$($W_j$の切片への効果)
- $\gamma_{10} = 3$(全体の傾き)
- $W_j = 1$(グループレベル変数)
- $x_{ij} = 2$(個体レベル変数)
本問のレベル2モデル:
$\beta_{0j} = \gamma_{00} + \gamma_{01}W_j + u_{0j} = 10 + 2W_j + u_{0j}$
$\beta_{1j} = \gamma_{10} + u_{1j} = 3 + u_{1j}$
Step 3: グループ$j$の回帰係数計算
$W_j = 1$のとき:
$\beta_{0j} = 10 + 2 \times 1 + u_{0j} = 12 + u_{0j}$
$\beta_{1j} = 3 + u_{1j}$
予測値を求める際は、ランダム効果の期待値$E[u_{0j}] = E[u_{1j}] = 0$を使用:
$E[\beta_{0j}] = 12$
$E[\beta_{1j}] = 3$
Step 4: 予測値の計算
レベル1モデルに代入:
$E[y_{ij}] = E[\beta_{0j}] + E[\beta_{1j}] \times x_{ij}$
$= 12 + 3 \times 2 = 12 + 6 = 18$
統合モデル(Combined Model)
レベル1とレベル2を統合すると:
$y_{ij} = \gamma_{00} + \gamma_{01}W_j + \gamma_{10}x_{ij} + u_{0j} + u_{1j}x_{ij} + r_{ij}$
固定効果部分:$\gamma_{00} + \gamma_{01}W_j + \gamma_{10}x_{ij}$
ランダム効果部分:$u_{0j} + u_{1j}x_{ij} + r_{ij}$
Step 5: 直接計算による確認
統合モデルを使用した計算:
$E[y_{ij}] = \gamma_{00} + \gamma_{01}W_j + \gamma_{10}x_{ij}$
$= 10 + 2 \times 1 + 3 \times 2 = 10 + 2 + 6 = 18$
Step 6: 階層線形モデルの特徴
HLMの主要な特徴
- 階層構造の考慮:データの階層性を明示的にモデル化
- グループ間変動:$u_{0j}, u_{1j}$によりグループ固有の効果を表現
- 交互作用効果:レベル間の交互作用(クロスレベル交互作用)を分析可能
- 分散成分分解:各レベルでの分散成分を推定
1. 固定効果と変量効果の解釈
- $\gamma_{00}$:全グループの平均切片
- $\gamma_{01}$:グループレベル変数が切片に与える効果
- $\gamma_{10}$:全グループの平均傾き
- $u_{0j}, u_{1j}$:グループ固有の偏差
2. 分散構造
$\text{Var}\begin{pmatrix} u_{0j} \\ u_{1j} \end{pmatrix} = \begin{pmatrix} \tau_{00} & \tau_{01} \\ \tau_{01} & \tau_{11} \end{pmatrix}$