Econometric Analysis of Panel Data: William Greene Department of Economics Stern School of Business
Econometric Analysis of Panel Data: William Greene Department of Economics Stern School of Business
http://people.stern.nyu.edu/wgreene/CumulantInstruments-Racicot-AE(2014)_46(10).pdf
The NYU
No Action
Letter
g
( )= i1z i (yi x
i1z i i
i )
N
N
1 Ni1i2zizi
1 Ni1ei2zi zi
Asy.Var[g
( )]
, estimated with
N
N
N
N
i ) '
i1zi (yi x
N
N
N
N
2
i1 i i i
1 N
z
(y
x
)
i .
N i1 i i
y1= x
1 1 1
2
y2= x
...
M
yM= x
x G ] 0
y1= x
1 1 1 , E[1 |x 1 ,x 2 ,...
2
y2= x
...
G
yM= x
2 ,...
G ,...
yi
, Xi
...
... ...
yiG
0 0
K1 K2
yi X+
i
i
0
1
... 0
, = 2 , i =
...
... ...
iG
... x
G
... KG columns
...
i1
i2
...
iG
Instruments
zi1 0
0 zi2
Zi
... ...
0 0
L1
L2
... 0
... 0
, G rows (1 for each equation)
... ...
... xiG
... LG columns
Such that
E zi1i1
zi1,1i1
0
zi1,2i1
0
...
...
i1,L1 i1
0
Moment Equations
zi1i1
z
i ] E i2 i2
E[Z
...
ziGiG
0 L1 rows
L rows
0 2
, for observation i
...
...
0 LG rows
Summing over i gives the orthogonality condition,
1 N
i E
i=1Z
zi1i1
0 L1
L
z
0 2
1 N i2 i2
i=1
...
...
N
0 LG
ziGiG
rows
rows
...
rows
Estimation-1
ig yig x
ig
M
m=1
(1/N)
N
i=1 ig,m
)
(yi xi
gg(g )'gg(g )
Leads to 2SLS
For all equations at the same time
= the minimizer of
G
g=1
g
'g
) g(
g( g
G
g=1
M
m=1
(1/N)
N
i=1 ig,m
(yi xi )
Estimation-2
Assuming ig are all uncorrelated, equation by equation GMM
1
ig
qg Ni1zig(yig x
N
)
g '
1 e z z 1 N
z
(y
) .
i1 ig
ig
ig g
N
N
2
i1 ig ig ig
Estimation-3
Combining GMM criteria
G
g 1
1 N
1 N
1
1 N 2
z
(y
x
)
'
z
z
iz1
ig g
i 1 ig ig ig
N i1 ig ig
N
N
Ni1zi1 (yi1 x
i1 1 ) Ni1i12 zi1zi1
0
N
N
2
z
(y
x
)
0
i 1 i2
i2
i2 2
i 1 i2zi2 zi2
q
'
...
...
...
N
z
(y
x
)
0
0
iG G
i1 iG iG
Ni1zi1 (yi1 x
i1 1 )
N
z
(y
x
)
i2 2
i1 i2 i2
...
N
z
(y
x
)
iG G
i1 iG iG
(y
ig
ig
ig g
)
...
0
...
0
...
...
... Ni1iG
ziGziG
Estimation-4
If disturbances are correlated across equations,
i1 1 )
Ni1zi1(yi1 x
z
(y
x
)
1
i1 i2
i2
i2 2
q
'
N
...
z
(y
x
)
i1 iG
iG
iG G
... Ni1i1iGzi1ziG
... Ni1i2iGzi2ziG
...
...
2
... Ni1iG
ziGziG
i1 1 )
Ni1zi1(yi1 x
N
i2 2 )
1 i1zi2 (yi2 x
N
...
iG G )
Ni1ziG (yiG x
Ni1i12 zi1zi1
Ni1i1i2zi1zi2
N
2
Ni1i2
zi2zi2
1 i1i2i1zi2zi1
N2
...
...
Ni1iGi1ziGzi1 Ni1iGi1ziGzi1
Estimation-5
If disturbances are correlated across equations,
G
G
(1/ N)Niz1
ig g) W gh
q g1 h1 (1/ N)Ni1zig (yig x
N 2
i1 i1 i1 i1
N
i1 i1 i2 i1 i2
N 2
i1 i2 i2 i2
... z z
... z z
...
...
N 2
... i1iGziGziG
N
i1 i1 iG i1 iG
N
i1 i2 iG i2 iG
ih
(yih x ih
h )
x2 it
i
z1
i
z2
it ui
NT
NT
NT
NT
1
1
1
XMD, so plim
X*'w plim
XMD 0
NT
NT
NT
1
1
plim
X*'X plim
X'MDX within groups sums of squares 0.
NT
NT
X* is a valid instrument.
plim b*=plim X* ' X
X
* y =
x2 it
i
z1
i
z2
it ui
zi1 zi2
zi1 xi1,1
L1 K1 columns
M
M
zi1 xi1,T i
estimates
and
(4) Transform
W*
i = [xit1 , xit2 ,zi1 ,zi2 ] - i[xi1 , xi2 ,zi1 ,zi2 ]
and
x1i
Application
Bias in Conventional Estimation
Development of Consistent Estimators
Efficient GMM Estimators
E[ci | Xi ] g( Xi )
No correlation across individuals
LSDV is Inconsistent
[(Steven) Nickell Bias]
yi,t yi ( xi,t xi )'+(yi,t1 yi ) (i,t i )
2 (T 1) T T
Cov[(yi,t1 yi ),(i,t i )] 2
T
(1 )2
Large when T is moderate or small.
Proportional bias for conventional T (5 - 15), is
on the order of 15% - 60% .
Instrumental Variables
Predetermined variables
, xi,2
yi,1 , xi,1
Zi
...
, xi,2
, xi,3
...
yi,1 , yi,2 , xi,1
...
0
...
...
, xi,2
,...xi, T 1
... yi,1 , yi,2 ,..., yi, T 2 , xi,1
...
0
(T rows)
Zi
...
0
...
, xi,2
,...xi, T 1 ...
yi,1 , yi,2 , xi,1
0
...
0
...
...
, xi,2
,...xi, T 1
... yi,1 , yi,2 ,..., yi, T 2 , xi,1
(T rows)
Simple IV Estimation
= N XZ N ZZ 1 N ZX
i=1 i i
i=1 i i
i=1 i i
N XZ N ZZ 1 N Zy
i=1 i i
i=1 i i
i=1 i i
N
N
]=
Est.Asy.Var[
i=1 XiZi i=1ZiZi
N
i=1
ZiXi
(y y )]2
Ni1 tT 3[( yi, t yi, t 1 ) ( xi, t xi, t 1 )'
i, t
i, t 1
Ni1 (Ti 2)
]= N XZ N Zv
Est.Asy.Var[
i=1 i i vi Zi
i=1 i i
N
i=1
ZiXi
Arellano/Bond
First Difference Formulation
yi, t 1 it
yit x
it
Parameters : =
[ , ]
The data
xi3
yi3
xi4
i4
yi
, Xi
...
M
iT i
xiT
K
yi,2 yi,1
yi,3 yi,2
, Ti -2 rows
yi,T yi, T1
1
columns
Arellano/Bond - GLS
yi,t yi,t1 ( xi,t xi,t1)'+(yi,t1 yi,t2 ) (i,t i )
i,3 i,2
Cov
2 1 0 ...
1 2 1 ...
0 1 2 ...
i,4 i,3
i,5 i,4
...
i,T i,T 1
...
...
0
0
0 2 i
1 ... 1
... 1 2
Ni=1ZiXi
Ni=1Z iy i
= N XZ N Z Z
i=1 i i i
i=1 i i
N XZ N Z
i iZ i
i=1
i
i
i=1
Z X
= XZ ZZ
X Z Z Z Z y
GMM Estimator
yi,t xi,t'+yi,t1 i,t
We make no assumptions about the disturbance. In first differences
yi,t yi,t1 ( xi,t xi,t1 )'+(yi,t1 yi,t2 ) (i,t i,t1)
(1) Two stage least squares
= N XZ N ZZ
i=1 i i
i=1 i i
Zi Xi
N
i=1
N
i=1
XiZi Ni=1ZiZi
1 Ni=1Ziv
1 Ni=1ZiXi
N
N
i=1
1 Ni=1Ziy i
XiZi W
1 Ni=1Zi Xi
N
Est.Asy.Var[
GMM ] i=1 XiZi W
Zi y i
N
i=1
Arellano/Bond/Bovers Formulation
Start with H&T
it
yit x1
x2 it
i
z1
i
z2
it ui
x1i
Let vit it ui
Let zit [( x1it - x1i )',( x2it - x2i )',z1i , x1']
Then E[zit vit ] 0
We formulate this for the Ti observations in group i.
Arellano/Bond/Bovers Formulation
Dynamic Model
it
yit yi,t1+x1
x2 it
i
z1
i
z2
it ui
Parameters : = [, 1
, 2
, 1
, 2']
The data
yi,2
yi,1 x1i2 x2i2 z1i z2 i
y
y
x1
x2
z1
z2
i3
i3
i
i,3
i,2
i
yi
, Xi
, Ti -1 rows
yi,T i
yi,T-1 x1iTi x2iTi z1i z2 i
1 K1
K2
L1 L2 columns
Arellano/Bond/Bovers Formulation
it
yit yi,t1 + x1
x2 it
i
z1
i
z2
it ui
Arellano/Bond/Bovers Formulation
zit [yi,1 yi ,..., yi,t1 yi ,( x1it - x1i )',( x2it - x2i )',z1i , x1']
Then E[zit vit ] 0
We formulate this for the last Ti -1 observations in group i.
(yi,1 , x1i2 ,x2i2 ,z1i )
(0,0,0)
(0,0,0)
...
(0,0,0)
(0,0,0)
(0,0,0)
...
(0,0,0)
Zi
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
(0,0,0)
...
(0,0,0)
1/(Ti 1)
1/(Ti 1)
, where MiD,(T -1) MDi without the last column.
i
...
1/(Ti 1)
Arellano/Bond/Bovers Formulation
For unbalanced panels the number of
columns for Zi varies. Given the form of
Zi, the number of columns depends on Ti.
We need all Zi to have the same number
of columns. For matrices with less
columns than the largest one, extra
columns of zeros are added.
Arellano/Bond/Bovers Formulation
The covariance matrix defines the model:
i=2 I - Classical (pooled) regression model (no effects)
i=2 I + u2ii' - Random effects model
i= A positive definite TxT
i
i matrix - GR model
Arellano/Bond/Bover Estimator
= N X H Z N Z H
H Z
i=1
i
i
i
i=1
i
i
i
i
i
N
N
H Z
i Hi
X
H
Z
Z
i=1
i
i
i
i=1
i
i
i
Zi Hi Xi
Ni=1 Zi Hi yi
N
i=1
= I . Compute residuals v
i yi Xi
(1) Use
i
= 1 N Hv
Then HH
i i i
i1 i iv iHi
N
.
(2) Recompute
]= N X H Z N Z H
H Z
Est.Asy.Var[
i=1
i
i
i
i=1
i
i
i
i
i
Zi Hi Xi
N
i=1
GMM Criterion
The GMM criterion which produces this estimator is
H Z
i i Zi Hi
iHZ
q= v
i
i
i
N
i1
N
i=1
i i i
ZHv
Application: Maquiladora
Maquiladora
Side Issue
How does y(t) = 1.220175 y(t-1) - 0.262198 y(t-2) + a behave?
y(t) = 1.220175 y(t-1) + a
is obviously explosive.
1.220175 0.262198
How to tell: A =
1
0
Postscript
x2 it
i
z1
i
z2
it ui