Calculating Joint Confidence Bands for Impulse Response Functions using Highest Density Regions