library(Ecdat)
data(Grunfeld)
pdata.frame(Grunfeld,"firm","year")
wi <- within(Grunfeld$inv)
w2 <- within(Grunfeld$inv,effect="twoways")
X <- model.matrix(~value+capital-1,Grunfeld)
Xt <- within(X,Grunfeld$firm)
Run the code above in your browser using DataLab