Monte-Carlo Impulse Response function for VAR
[irf_low,irf_upp]=irf_mc1(res,mres,S,N,p,niter,siz)
results = results tlist returned by VAR
mres = decomposition method
N = # of endogenous variables in the VAR
S = # of periods
p = # of lags in the VAR
niter = # of draws in the simulation
siz = size of the confidence band
irf_low = ((S+1) x T) lower range of impulse response confidence band
irf_upp = ((S+1) x T) upper range of impulse response confidence band
[irf_low,irf_upp]=irf_mc1(results,mres,S,N,p,niter,siz) // Example taken from function irf(). This is equivalent to the simpler call: irf(results,S,'mres=name of mres','size=siz','meth=mc1','niter=value of niter')).