Calculating joint confidence bands for impulse response functions using highest density regions