generate a VAR model and provide impulse response functions, variance decomposition and granger causality test

We have 6 "agents" with 4 variables all being quarterly observed We need to set up a VAR model and test for stationarity and heteroskedasticity etc, and then produce impulse response functions. We also need to test for Granger causality and proceed with variance decomposition.