I estimated VAR (4 variables with 2 lags) in EViews and then in R using library vars:
mydata = data.frame(lgdp_