的 GLS 估计量才是最佳线性无偏估计(BLUE)。 由计算分块矩阵逆的引理1容易证明，
( ) ( ) 1引理：对于矩阵 A、B 和 D， A + BDB' −1 = A−1 − A−1B B' A−1B + D−1 −1 B' A−1 .
面板数据计量分析 白仲林
( ) X
' i
ωi−1
X
i
=
⎡ ⎢⎣
ε m = ( IN ⊗ιT )ξ m + um
E
⎛ ⎜ ⎝
ξ u
m m
⎞ ⎟ ⎠
(ξ
l
'
ul
')
=
⎛ ⎜
σ
2 ξ ,ml
I
N
⎝
⎞
σ I 2 u,ml NT
⎟ ⎠
，
m，l = 1, 2, …, M
ξ ~（0，∑ ξ ⊗ IN）
u ~（0，∑u ⊗ INT），
( ) ( ) 其中，∑ ξ =
σ2 ξ ,ml
"
xKi2
⎥ ⎥
# " #⎥
⎡X1
⎢ ， X =⎢⎢
X2 %
⎤ ⎥ ⎥， ⎥
⎢⎥ ⎣YN ⎦ NT×1
⎢⎥ ⎣ X N ⎦ NT×K
⎢ ⎣x1iT
x2iT
"
⎥ x ⎦ KiT T×K
⎢ ⎣
⎥ XN⎦
⎡ ξ1 ⎤
⎡ξ1i ⎤
⎡ u1 ⎤
ξ
=
⎢ ⎢
ξ
2
⎢#
⎥ ⎥ ⎥
， ξi
=
⎢⎢ξ2i ⎢#
⎥ ⎥ ⎥
sureg (I1 = F1 C1) (I2 = F2 C2[, noconstant]) (I3 = F3 C3) (I4 = F4 C4) (I5 = F5 C5)[ , isure]
. sureg (I1 = F1 C1) (I2 = F2 C2) (I3 = F3 C3) (I4 = F4 C4) (I5 = F5 C5)
(占总资产)比率是影响工业电力需求的重要因素。但是，考虑到宏观或产业政策的相互影响， 使得这三个模型的误差项具有相关性，因此需要将它们联立估计。
这时，面板数据模型应该设定为
K
∑ y
m it
=
β
m k
x
m kit
+
εm it
i = 1,2," , N ； t = 1,2," ,T
k =1
m = 1，2，…, M
.0059292 .2278877
C2
.4503213 .1218427
3.70 0.000
.211514 .6891286
_c ons
47.17259 114.8141
0.41 0.681
-177.859 272.2041
I3
F3
.0352793 .0127776
2.76 0.006
.0102357 .0603228
46.73 0.0000
I4
20
2 12.34291 0.9122
210.78 0.0000
I5
20
2 8.391486 0.6778
39.19 0.0000
迭代估计法选项
Coef. Std. Err.
z
P>|Leabharlann z|[95% Conf. Interval]
I1
F1
.1288887 .0212979
6.05 0.000
假设 3：对每个个体 i，误差向量 Ui 是均值为零、具有协方差矩阵为 σi2 IT 的独立同分布
( ) 随机向量，即， E (ui ) = 0 ， E
ui u'j
=
⎧σ ⎨
2 i
IT
⎩0
(i = j) (i ≠ j) .
假设 4：模型（5.3）的系数向量 βi 是均值 β 和协方差矩阵 Σ 的独立同分布随机向量，
面板数据计量分析 白仲林
( ) ( ) 类似地，设 y' = y1' " yM ' ， ym = y1m1 " y1mT " yNm1 " yMmT '
⎡ X1
⎢ X =⎢
⎢
X2 %
⎤
⎡
Xm 1
⎤
⎥ ⎥，
其中， X m
=
⎢ ⎢
Xm 2
⎥ ⎥
⎥
⎢⎥
Xm i
=
⎛ ⎜ ⎜ ⎜
xm 1i1
xm 1i 2 #
（5.4）
ωi = E ⎡⎣(Xiξi + ui )(Xiξi + ui )' ⎤⎦ = E ⎡⎣(Xiξi )(Xiξi )' + uiui' ⎤⎦ ，
ωi
=
X
i
ΣX
' i
+
σ
2 i
IT
i = 1,2,3," , N .
5.2 Swamy 随机系数模型的估计
由于模型（5.3）的合并随机项 Xξ +U 存在异方差和序列相关性，所以，模型（5.3）
随机向量，并且， β 和 Σ 分别是不随时间变化的向量和矩阵，则面板数据随机系数模型可
写为矩阵形式
Y = Xβ + Xξ + U
（5.3）
⎡ Y1 ⎤
其中，Y
=
⎢ ⎢
Y2
⎢#
⎥ ⎥ ⎥
⎡ X1 ⎤
，
X
=
⎢ ⎢ ⎢
X2 #
⎥ ⎥ ⎥
⎡ x1i1
，
Xi
=
⎢⎢x1i ⎢#
2
x2i1 " xKi1 ⎤
x2i2
( ) 即，
βi
=
β
+ξ i
， E (ξi )
=
0；E
ξi
ξ
' j
⎧Σ
=
⎨ ⎩
0
(i (i
= ≠
j) j)
；
假设 5：对任意的 i 和 j，误差向量 ui 与系数向量 β j 独立。 于是，模型（5.3）的合并随机项 Xξ +U 的协方差矩阵
( )( ) Ω
=
E
⎡ ⎢⎣
Xξ +U
Xξ +U
⎡⎡
.3141704 .0260555 12.06 0.000
.2631025 .3652383
_c ons
.6961959 11.57599
0.06 0.952 -21.99233 23.38473
I5
F5
.1444101 .0501274
2.88 0.004
.0461623
.242658
C5
.0069288 .0192621
i =1
（5.5）
∑ ( ) ( ) ( ) 其中，Wi
=
⎧ ⎨ ⎩
N j =1
⎡ ⎣
Σ
+σ 2 j
X
' j
X
j
−1
⎤ −1 ⎦
⎫−1 ⎬ ⎭
⎡ ⎣
Σ
+σ 2 i
X
' i
X
i
⎤ −1 −1 ⎦
，
βˆ i
=
X
' i
X
i
−1
X
Y'
ii
是第
i
个个体模型系数向量 βi 的 OLS 估计。
由（5.5）式可见，模型（5.3）的 GLS 估计 βˆ 是 βi 的 OLS 估计 βˆ i 的矩阵加权平均。并 且， βˆ 的协方差矩阵
"
β
m K
' ，m = 1, 2, …, M
显然，模型（4.3）是标准的线性 SUR 模型（Zellner，1962）。从而，可以基于 SUR
模型的两种估计方法（FGLS 估计量和迭代估计法 ITERZEF）估计模型（4.1）。
类似地，也可以讨论双因素误差的模型（4.1），即，
ε
m it
= ξim
+ λtm
面板数据计量分析 白仲林
第三讲 面板数据变系数回归模型
1 确定性变系数模型——SUR 模型
单因素面板数据 SUR 模型 前面所讨论的面板数据模型属于面板数据单方程模型，常常也需要建立多方程的面板数
据模型。这里只讨论一种最简单的情形，即假设在 M 个方程中，每个单方程模型的解释变 量各不相同。对于更一般的情况，请参考 Baltagi(2008，P121-141).
+ uimt 的情形。
SUR 模型的检验
Breusch 和 Pagan（1980）基于 Lagrange 乘数（Lagrange multiplier）方法提出了检验零
假设
H0： Ω 是对角矩阵
的 LM 统计量。
不含截距选
SUR 模型的 Stata 估计
以 Grunfeld(1958)数据的前 5 家公司数据为例。 Stata 命令：
0.36 0.719 -.0308242 .0446818
_c ons
25.00319 6.239317
4.01 0.000
12.77435 37.23202
面板数据计量分析 白仲林
2 面板数据随机系数模型
自 Swamy（1970、1973 和 1978 等）应用面板数据的随机系数模型研究美国各州汽油需 求函数等问题以来，面板数据的随机系数模型得到了一些应用。然而，由于该类模型的参数 估计计算比较复杂，制约了它的广泛应用，经验研究主要集中于随机效应模型的使用。但是 这并不意味着随机系数模型不重要，实际上，在研究经济增长收敛理论（Durlauf，2001） 等许多经济问题时，建立面板数据随机系数模型是解决问题的合理方法（Canova，1999）。 本节主要介绍两种面板数据随机系数模型，一种是 Swamy 随机系数模型，另一种是 Hsiao 随机系数模型。