Joint Prediction Regions for time-series models