3.2.2 Random E ff ect Model (
ランダム効果モデル) Model:
y
it= X
itβ + v
i+ u
it, i = 1 , 2 , · · · , n , t = 1 , 2 , · · · , T where i indicates individual and t denotes time.
The assumptions on the error terms v
iand u
itare:
E(v
i| X) = E(u
it| X) = 0 for all i ,
V(v
i| X) = σ
2vfor all i , V(u
it| X) = σ
2ufor all i and t ,
Cov(v
i, v
j| X) = 0 for i , j , Cov(u
it, u
js| X) = 0 for i , j and t , s , Cov(v
i, u
jt| X) = 0 for all i, j and t .
Note that X includes X
itfor i = 1 , 2 , · · · , n and t = 1 , 2 , · · · , T .
In a matrix form with respect to t = 1 , 2 , · · · , T , we have the following:
y
i= X
iβ + v
i1
T+ u
i, i = 1 , 2 , · · · , n ,
where y
i=
y
i1y
i2...
y
iT
, X
i=
X
i1X
i2...
X
iT
and u
i=
u
i1u
i2...
u
iT
are T × 1, T × k and T × 1, respectively.
u
i∼ N(0 , σ
2uI
T) and v
i1
T∼ N(0 , σ
2v) = ⇒ v
i1
T+ u
i∼ N(0 , σ
2v1
T1
0T+ σ
2uI
T) . Again, in a matrix form with respect to i, we have the following:
y = X β + v + u ,
where y =
y
1y
2...
, X =
X
1X
2...
, v =
v
11
Tv
21
T...
and u =
u
1u
2...
are nT × 1, nT × k, nT × 1 and
nT × 1, respectively.
The distribution of u + v is given by:
v + u ∼ N (
0 , I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T) ) The likelihood function is given by:
L( β, σ
2v, σ
2u) = (2 π )
−nT/2I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T)
−1/2× exp (
− 1
2 (y − X β )
0(
I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T) )
−1(y − X β ) ) . Remember that f (x) = (2 π )
−k/2|Σ|
−1/2exp (
−
12(x − µ )
0Σ
−1(x − µ ) )
when X ∼ N( µ, Σ ),
where X denotes a k-variate random variable.
The estimators of β , σ
2vand σ
2uare given by maximizing the following log-likelihood function:
log L( β, σ
2v, σ
2u) = − nT
2 log(2 π ) − 1
2 log I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T)
− 1
2 (y − X β )
0(
I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T) )
−1(y − X β ) . MLE of β , denoted by ˜ β , is given by:
β ˜ = ( X
0(
I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T) )
−1X )
−1X
0(
I
n⊗ ( σ
2v1
T1
0T+ σ
2uI
T) )
−1y
= ( ∑
ni=1
X
i0( σ
2v1
T1
0T+ σ
2uI
T)
−1X
i)
−1( ∑
ni=1