Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter July 25, 2017

Generating prediction bands for path forecasts from SETAR models

  • Daniel Grabowski , Anna Staszewska-Bystrova and Peter Winker ORCID logo EMAIL logo

Abstract

Prediction bands for time series are usually generated point-wise by bootstrap methods. Such bands only convey the prediction uncertainty for each horizon separately. The joint distribution is not taken into account. To represent the forecast uncertainty over the entire horizon, methods for constructing joint prediction bands for path forecasts from SETAR models are proposed. This class of nonlinear models is increasingly used in time series analysis and forecasting as it is useful for capturing nonlinear dynamics. Approaches based on statistical theory and explicit sequential and global optimization methods are both considered. Monte Carlo simulation is used to assess the performance of the proposed methods. The comparison is done with regard to the actual coverage of the constructed prediction bands for full path forecasts as well as with regard to the width of the bands. An empirical application demonstrates the relevance of the choice of bands for indicating the uncertainty of path forecasts in nonlinear models.

Acknowledgements

Support from the MNiSW/DAAD PPP Grant (56268818) and the National Science Center, Poland (NCN) through HARMONIA 6:UMO-2014/14/M/HS4/00901 are gratefully acknowledged.

A Bootstrap procedure for residual rescaling

Algorithm 3:

Algorithm 3: Bootstrap procedure for residual rescaling

1: Estimate SETAR parameters β^ij(r) and γ^ij(r) and obtain residuals ε^t(r1);

2: Compute the standard deviation σ(r1) of the residuals ε^t(r1);

3: Re-center ε^t(r1);

4: forc = 1 to Cdo

5:  Construct bootstrapped series ytr similarly to Eq. (4), using β^ijr, γ^r and a draw from the re-centered ε^t(r1) ;

6:  Estimate the SETAR parameters for ytr and obtain the residuals ε^t(r2);

7:  Compute the standard deviation σc(r2) of ε^t(r2);

8: end for

9: The estimated rescaling factor is given by r=σ(r1)(1Cσc(r2))1;

B Simulation results for heteroscedastic error terms

Table 6:

Estimated coverage probabilities and band width and respective standard deviations for the joint coverage of the prediction bands for DGPs with regime dependent standard deviation. The effect of different standard deviations is shown.

β11 = 0.8, β21 = 0.8, T = 60
σε1 = 0.9, σε2 = 1.1σε1 = 0.8, σε2 = 1.2σε1 = 0.5, σε2 = 1.5
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI82.9664.7982.4665.1480.0464.68
(11.22)(20.50)(12.19)(21.32)(16.22)(26.27)
Bt96.39122.6695.82123.7093.06144.06
(4.29)(141.22)(5.23)(101.13)(9.63)(205.52)
Badj94.3189.1793.5090.0190.1692.02
(5.83)(29.63)(7.12)(31.19)(11.43)(40.80)
TA94.9984.1994.2184.7990.7685.57
(5.88)(22.29)(7.12)(23.48)(11.41)(29.94)
  1. BPI, Bt, Badj and TA stand for the bootstrap interval method, the traditional Bonferroni method, the adjusted Bonferroni method and the method using Threshold Accepting, respectively. The results were obtained using 2000 Monte Carlo replications and 2000 bootstrap replications.

Table 7:

Estimated coverage probabilities and band width and respective standard deviations for the joint coverage of the prediction bands for DGPs with regime dependent standard deviation. The effect of different sample sizes is shown.

β11= 0.8, β21 = 0.8, σε1 = 0.8, σε2 = 1.2
T = 60T = 120T = 240
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI82.4665.1480.5352.9878.6248.78
(12.19)(21.32)(10.58)(8.68)(9.63)(4.71)
Bt95.82123.7095.8378.2095.7270.37
(5.23)(101.13)(4.38)(13.80)(3.58)(6.95)
Badj93.5090.0193.8870.8293.7765.42
(7.12)(31.19)(5.37)(10.91)(4.36)(5.94)
TA94.2184.7993.9770.6993.7865.42
(7.12)(23.48)(5.35)(10.68)(4.36)(5.94)
  1. BPI, Bt, Badj and TA stand for the bootstrap interval method, the traditional Bonferroni method, the adjusted Bonferroni method and the method using Threshold Accepting, respectively. The results were obtained using 2000 Monte Carlo replications and 2000 bootstrap replications.

Table 8:

Estimated coverage probabilities and band width and respective standard deviations for the joint coverage of the prediction bands for DGPs with regime dependent standard deviation. The effect of different autoregressive coefficients is shown.

β21 = 0.8, σε1 = 0.8, σε2 = 1.2, T = 60
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI78.5955.1182.4665.1488.45355.69
(13.31)(39.93)(12.19)(21.32)(18.75)(5131.30)
Bt94.4589.7695.82123.7097.271062.70
(7.11)(82.15)(5.23)(101.13)(9.91)(15064.0)
Badj92.6274.7793.5090.0194.57472.79
(7.92)(53.25)(7.12)(31.19)(14.49)(6260.80)
TA92.8471.7994.2184.7995.33338.84
(7.99)(36.91)(7.12)(23.48)(13.99)(4453.30)
  1. BPI, Bt, Badj and TA stand for the bootstrap interval method, the traditional Bonferroni method, the adjusted Bonferroni method and the method using Threshold Accepting, respectively. The results were obtained using 2000 Monte Carlo replications and 2000 bootstrap replications.

C Simulation results for alternative threshold values

Table 9:

Estimated coverage probabilities and band width and respective standard deviations for the joint coverage of the prediction bands for DGPs with threshold value γ = 0.5.

sample size: T=60
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI78.9949.4683.264.9185.69428.98
(11.83)(17.41)(11.13)(19.85)(20.58)(3286.50)
Bt95.3778.8596.51122.5397.031391.4
(5.14)(39.35)(4.17)(105.18)(10.85)(11734.0)
Badj93.7566.68294.5489.7093.48638.86
(6.09)(23.59)(5.58)(29.71)(16.33)(5169.60)
TA93.9964.87795.2384.1394.05394.41
(6.11)(18.88)(5.66)(21.59)(16.32)(2541.60)
sample size: T=120
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI76.0443.0680.9952.4386.17128.03
(8.47)(6.27)(8.12)(7.89)(17.11 )(159.95)
Bt95.4960.9496.4976.1797.82217.35
(3.44)(9.51)(3.08)(12.03)(7.36)(317.82)
Badj94.2357.5694.7969.2594.52162.93
(4.00)(7.89)(3.78)(9.46)(12.75)(209.89)
TA94.2557.5194.8569.1494.62147.67
(4.00)(7.66)(3.81)(9.59)(13.25)(148.2)
sample size: T=240
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI75.0141.6278.3948.1584.46100.42
(6.34)(3.78)(6.14)(4.02)(16.52 )(62.80)
Bt95.8258.5996.2768.2497.40159.74
(2.43)(5.49)(2.26)(5.87)(8.73)(104.24)
Badj94.6356.0094.4663.7093.94125.50
(2.81)(4.87)(2.73)(5.05)(12.50)(71.99)
TA94.6355.9994.4663.6993.75119.65
(2.81)(4.81)(2.74)(5.05)(12.28)(63.47)
  1. BPI, Bt, Badj and TA stand for the bootstrap interval method, the traditional Bonferroni method, the adjusted Bonferroni method and the method using Threshold Accepting, respectively. The results were obtained using 2000 Monte Carlo replications and 2000 bootstrap replications.

Table 10:

Estimated coverage probabilities and band width and respective standard deviations for the joint coverage of the prediction bands for DGPs with threshold value γ = 1.

sample size: T=60
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI79.8551.0883.9268.6787.39253.05
(12.16)(18.70)(11.54)(25.17)(16.90)(1477.50)
Bt95.5883.7496.66135.8297.57752.55
(5.13)(50.97)(4.31)(180.87)(7.74)(5380.20)
Badj94.1069.2594.8495.0794.99363.50
(6.04)(26.86)(5.67)(37.88)(11.43)(2346.40)
TA94.3766.9995.6088.4495.48254.88
(6.06)(19.86)(5.52)(28.04)(11.8)(1195.00)
sample size: T=120
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI76.4844.2180.8652.7285.78120.74
(8.64)(7.03)(8.33)(8.21)(17.116 )(124.88)
Bt95.6562.3796.4876.5297.77200.56
(3.45)(10.16)(3.08)(12.62)(7.13)(245.08)
Badj94.5059.0194.7469.6394.43152.38
(3.92)(8.67)(3.91)(10.16)(12.14)(152.09)
TA94.5258.9594.8269.5294.52140.90
(3.92)(8.34)(3.93)(10.01)(12.87)(119.55)
sample size: T=240
β11 = 0.2β11 = 0.8β11 = 1
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
Coverage

(SD)
Width

(SD)
BPI74.8942.3178.8048.3684.8097.56
(6.14)(3.58)(6.14)(4.00)(15.88 )(59.68)
Bt95.9059.4696.3868.5097.75153.98
(2.31)(5.18)(2.17)(5.77)(8.31)(101.97)
Badj94.8056.8994.6463.9594.08121.83
(2.61)(4.65)(2.63)(4.98)(12.18)(68.20)
TA94.8056.8994.6463.9494.01116.77
(2.61)(4.65)(2.63)(4.98)(12.76)(60.15)
  1. BPI, Bt, Badj and TA stand for the bootstrap interval method, the traditional Bonferroni method, the adjusted Bonferroni method and the method using Threshold Accepting, respectively. The results were obtained using 2000 Monte Carlo replications and 2000 bootstrap replications.

References

Arora, S., M. A. Little, and P. E. McSharry. 2013. “Nonlinear and Nonparametric Modeling Approaches for Probabilistic Forecasting of the US Gross National Product.” Studies in Nonlinear Dynamics and Econometrics 17: 395–420.10.1515/snde-2012-0029Search in Google Scholar

Bai, J.. 1997. “Estimating Multiple Breaks One at a Time.” Econometric Theory 13: 315–352.10.1017/S0266466600005831Search in Google Scholar

Bec, F., O. Bouabdallah, and L. Ferrara. 2014. “The Way Out of Recessions: A Forecasting Analysis for Some Euro Area Countries.” International Journal of Forecasting 30: 539–549.10.1016/j.ijforecast.2013.01.004Search in Google Scholar

Boero, G., and F. Lampis. 2017. “The Forecasting Performance of SETAR Models: An Empirical Application.” Bulletin of Economic Research 69: (3): 216–228. .10.1111/boer.12068Search in Google Scholar

Chan, K.. 1993. “Consistency and Limiting Distribution of the Least Squares Estimator of a Threshold Autoregressive Model.” The Annals of Statistics 21: 520–533.10.1214/aos/1176349040Search in Google Scholar

Claveria, O., and S. Torra. 2014. “Forecasting Tourism Demand to Catalonia: Neural Networks vs. Time Series Models.” Economic Modelling 36: 220–228.10.1016/j.econmod.2013.09.024Search in Google Scholar

Clements, M., and D. Hendry. 1998. Forecasting Economic Time Series Cambridge: Cambridge University Press.10.1017/CBO9780511599286Search in Google Scholar

Clements, M. P., P. H. Franses, J. Smith, and D. Van Dijk. 2003. “On SETAR Non-linearity and Forecasting.” Journal of Forecasting 22: 359–375.10.1002/for.863Search in Google Scholar

Clements, M. P., and J. Smith. 1997. “The Performance of Alternative Forecasting Methods for SETAR Models.” International Journal of Forecasting 13: 463–475.10.1016/S0169-2070(97)00017-4Search in Google Scholar

Clements, M. P., and J. Smith. 1999. “A Monte Carlo Study of the Forecasting Performance of Empirical SETAR Models.” Journal of Applied Econometrics 14: 124–141.10.1002/(SICI)1099-1255(199903/04)14:2<123::AID-JAE493>3.0.CO;2-KSearch in Google Scholar

Clements, M. P., and J. Smith. 2000. “Evaluating the Forecast Densities of Linear and Non-linear Models: Applications to Output Growth and Unemployment.” Journal of Forecasting 19: 255–276.10.1002/1099-131X(200007)19:4<255::AID-FOR773>3.0.CO;2-GSearch in Google Scholar

Cryer, J., and K.-S. Chan. 2008. Time Series Analysis With Applications in R New York, NY: Springer.10.1007/978-0-387-75959-3Search in Google Scholar

Dueck, G., and T. Scheuer. 1990. “Threshold Accepting: A General Purpose Algorithm Appearing Superior to Simulated Annealing.” Journal of Computational Physics 90: 161–175.10.1016/0021-9991(90)90201-BSearch in Google Scholar

Feng, H.. 2011. “Forecasting Comparison Between Two Nonlinear Models: Fuzzy Regression Versus SETAR.” Applied Economics Letters 18: 1623–1627.10.1080/13504851.2011.554369Search in Google Scholar

Feng, H., and J. Liu. 2003. “A SETAR Model for Canadian GDP: Non-linearities and Forecast Comparisons.” Applied Economics 35: 1957–1964.10.1080/0003684032000160674Search in Google Scholar

Fenga, L., and D. N. Politis. 2015. “Bootstrap Order Selection for SETAR Models.” Journal of Statistical Computation and Simulation 85: 235–250.10.1080/00949655.2013.811507Search in Google Scholar

Kapetanios, G., and Y. Shin. 2006. “Unit Root Tests in Three-Regime SETAR Models.” Econometrics Journal 9: 252–278.10.1111/j.1368-423X.2006.00184.xSearch in Google Scholar

Kilian, L.. 1998. “Small-Sample Confidence Intervals for Impulse Response Functions.” Review of Economics and Statistics 80: 218–230.10.1162/003465398557465Search in Google Scholar

Li, J.. 2006. “Testing Granger Causality in the Presence of Threshold Effects.” International Journal of Forecasting 22: 771–780.10.1016/j.ijforecast.2006.01.003Search in Google Scholar

Li, J. 2011. “Bootstrap Prediction Intervals for SETAR Models.” International Journal of Forecasting 27: 320–332.10.1016/j.ijforecast.2010.01.013Search in Google Scholar

Lütkepohl, H., A. Staszewska-Bystrova, and P. Winker. 2015a. “Comparison of Methods for Constructing Joint Confidence Bands for Impulse Response Functions.” International Journal of Forecasting 31: 782–798.10.1016/j.ijforecast.2013.08.003Search in Google Scholar

Lütkepohl, H., A. Staszewska-Bystrova, and P. Winker. 2015b. “Confidence Bands for Impulse Responses: Bonferroni vs. Wald.” Oxford Bulletin of Economics and Statistics 77: 800–821.10.1111/obes.12114Search in Google Scholar

Petruccelli, J.. 1986. “On the Consistency of Least Squares Estimators of a Threshold AR(1) Model.” Journal of Time Series Analysis 7: 269–278.10.1111/j.1467-9892.1986.tb00494.xSearch in Google Scholar

Petruccelli, J. D., and S. W. Woolford. 1984. “A Threshold AR(1) Model.” Journal of Applied Probability 21: 270–286.10.1017/S0021900200024670Search in Google Scholar

Potter, S. M.. 1995. “A Nonlinear Approach to US GNP.” Journal of Applied Econometrics 10: 109–125.10.1002/jae.3950100203Search in Google Scholar

Shaman, P., and R. A. Stine. 1988. “The Bias of Autoregressive Coefficient Estimators.” Journal of the American Statistical Association 83: 842–848.10.1080/01621459.1988.10478672Search in Google Scholar

Staszewska-Bystrova, A.. 2011. “Bootstrap Prediction Bands for Forecast Paths from Vector Autoregressive Models.” Journal of Forecasting 30: 721–735.10.1002/for.1205Search in Google Scholar

Staszewska-Bystrova, A., and P. Winker. 2013. “Constructing Narrowest Pathwise Bootstrap Prediction Bands Using Threshold Accepting.” International Journal of Forecasting 29: 221–233.10.1016/j.ijforecast.2012.09.004Search in Google Scholar

Staszewska-Bystrova, A., and P. Winker. 2016. “Improved Bootstrap Prediction Intervals for SETAR Models.” Statistical Papers 57: 89–98.10.1007/s00362-014-0643-1Search in Google Scholar

Stine, R. A.. 1987. “Estimating Properties of Autoregressive Forecasts.” Journal of the American Statistical Association 82: 1072–1078.10.1080/01621459.1987.10478542Search in Google Scholar

Teräsvirta, T.. 2006.. “Univariate Nonlinear Time Series Models.” In: edited by T. Mills, and K. Patterson (Eds.),. Palgrave Handbook of Econometrics 1, 396–424. Basingstoke: Palgrave Macmillan.Search in Google Scholar

Thombs, L. A., and W. R. Schucany. 1990. “Bootstrap Prediction Intervals for Autoregression.” Journal of the American Statistical Association 85: 486–492.10.1080/01621459.1990.10476225Search in Google Scholar

Tiao, G. C., and R. S. Tsay. 1994. “Some Advances in Non-linear and Adaptive Modelling in Time-Series.” Journal of Forecasting 13: 109–131.10.1002/for.3980130206Search in Google Scholar

Tong, H.. 1990. Non-Linear Time Series: A Dynamical System Approach Oxford: Oxford University Press.Search in Google Scholar

Winker, P.. 2001. Optimization Heuristics in Econometrics: Applications of Threshold Accepting Chichester: Wiley.Search in Google Scholar

Wolf, M., and D. Wunderli. 2015. “Bootstrap Joint Prediction Regions.” Journal of Time Series Analysis 36: 352–376.10.1111/jtsa.12099Search in Google Scholar


Supplemental Material

The online version of this article offers supplementary material (DOI:https://doi.org/10.1515/snde-2016-0066).


Published Online: 2017-7-25

©2017 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 7.12.2023 from https://www.degruyter.com/document/doi/10.1515/snde-2016-0066/html
Scroll to top button