residuals from a VAR estimation
[resid]=vare(y,nlag,x)
* y = an (nobs x neqs) matrix of y-vectors
* nlag = the lag length
* x = optional matrix of variables (nobs x nx)
(NOTE: constant vector automatically included)
* resid = matrix of residuals (nobs x neqs)