vare — 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)