95% prediction interval for the VAR model can be obtained this way (vars library):
library(vars) data(Canada) series_mod = wind