data("SURdata", package="panelSUR")
## Data preparation
library(plm)
datap <- pdata.frame(data, index=c("IND", "TIME"))
## Equations specification
eq1<-Y1~X1+X2
eq2<-Y2~X1+X2+X3
eqlist<-c(eq1,eq2)
## System estimation
mod1<-SURest(eqlist=eqlist,method="1wayWB",data=datap)
## Summary of estimation results
printSUR(mod1)
Run the code above in your browser using DataLab