# Anderson and Hsiao estimator
library(Ecdat)
data(Grunfeld)
pdata.frame(Grunfeld,"firm","year")
form <- diff(inv)~diff(inv,1)+diff(value)+diff(capital)
hs <- plm(form,data=Grunfeld,instruments=~lag(inv,2),endog=~diff(inv,1),model="pooling")
Run the code above in your browser using DataLab