Wavelet Transform as an Alternative to Power Transformation in Time Series Analysis

C. Ogbonna, C. Nweke, Eleazer C. Nwogu, I. Iwueze
{"title":"Wavelet Transform as an Alternative to Power Transformation in Time Series Analysis","authors":"C. Ogbonna, C. Nweke, Eleazer C. Nwogu, I. Iwueze","doi":"10.18052/WWW.SCIPRESS.COM/BMSA.17.57","DOIUrl":null,"url":null,"abstract":"This study examines the discrete wavelet transform as a transformation technique in the analysis of non-stationary time series while comparing it with power transformation. A test for constant variance and choice of appropriate transformation is made using Bartlett’s test for constant variance while the Daubechies 4 (D4) Maximal Overlap Discrete Wavelet Transform (DWT) is used for wavelet transform. The stationarity of the transformed (power and wavelet) series is examined with Augmented Dickey-Fuller Unit Root Test (ADF). The stationary series is modeled with Autoregressive Moving Average (ARMA) Model technique. The model precision in terms of goodness of fit is ascertained using information criteria (AIC, BIC and SBC) while the forecast performance is evaluated with RMSE, MAD, and MAPE. The study data are the Nigeria Exchange Rate (2004-2014) and the Nigeria External Reserve (1995-2010). The results of the analysis show that the power transformed series of the exchange rate data admits a random walk (ARIMA (0, 1, 0)) model while its wavelet equivalent is adequately fitted to ARIMA (1,1,0). Similarly, the power transformed version of the External Reserve is adequately fitted to ARIMA (3, 1, 0) while its wavelet transform equivalent is adequately fitted to ARIMA (0, 1, 3). In terms of model precision (goodness of fit), the model for the power transformed series is found to have better fit for exchange rate data while model for wavelet transformed series is found to have better fit for external reserve data. In forecast performance, the model for wavelet transformed series outperformed the model for power transformed series. Therefore, we recommend that wavelet transform be used when time series data is non-stationary in variance and our interest is majorly on forecast. 1.0 Introduction In several organizations, managerial decisions are largely based on the available information of the past and present observations and possibly on the process that generate such observations. A time series data provides such information. Time series is used to represent the characterized time course of behavior of wide range of several systems which could be biological, physical or economical. The utility of the time series data lies in the result of the time series analysis. Such analysis will be helpful in achieving the aim for collection of such data which could be for description (exposing the main properties of a series), explanation (revealing the relationship between variables of a series especially when observations are taken on two or more variables), forecasting (prediction of the future values of a series) and control (taking appropriate corrective actions) [1]. To analyze any time series data, time series analysis techniques are adopted. The commonly used techniques are: descriptive technique, probability models technique and spectral density analysis technique. The inference based on the descriptive method and probability models is often referred to as analysis in time domain while inference based on spectral density function is referred to analysis in frequency domain [2,3,4,5]. Bulletin of Mathematical Sciences and Applications Submitted: 2016-09-12 ISSN: 2278-9634, Vol. 17, pp 57-74 Revised: 2016-09-25 doi:10.18052/www.scipress.com/BMSA.17.57 Accepted: 2016-10-24 2016 SciPress Ltd, Switzerland Online: 2016-11-01 SciPress applies the CC-BY 4.0 license to works we publish: https://creativecommons.org/licenses/by/4.0/ All these models assume that the error component of the series (et) is normally distributed with zero mean and constant variance   2  or that the series is normally distributed with constant mean (μ) and constant variance   2  . When any study data violates any or all of these assumptions, the series is subjected to transformation. Transformation helps to (i) stabilize the variance of a series, (ii) make the seasonal effect when present additive and (iii) make the data normally distributed [1]. One of the transformations commonly used is the power transformation developed by [6]. [7] noted that the power transformation (i) changes the scale of the original series, (ii) may introduce bias in the forecast especially when data have to be transformed back to its original scale and (iii) often the transformed series have no physical interpretation. [1] argued that transformation alone may not be helpful when variance changes through time in the absence of trend. In such case (i.e., when variance changes through time in the absence of trend), he recommended that a model that allows for changes in variance should be considered. Wavelet method is one such method that allows for changes in variance which has been found to be useful in time series analysis. It involves decomposition, de-nosing and reconstruction of series. Decomposition involves breaking down time series into two main components namely the detail and the smooth components, de-noising deals with the removal of the non significant components of the series, while reconstruction involves recovering of the original series devoid of noise. Unlike the power transformation, wavelet method does not change the scale of the series, poses no problem in its interpretation and does not rely on the assumption of any underlying distribution of the study data. Additionally, wavelet transformation method allow for decomposition of a series without knowing the underlying functional form of the series [8]. Could this lead to an improved model and forecast performance from those based on power transformation? This and other related questions are what this study intends to address. Therefore, the objective of this study is to examine the precision (in terms of goodness-of-fit) and forecast performances of the models for wavelet transformed series while comparing it with models for power transformed time series. 2.0 Literature Review Various transformations exist, but the power transformation developed by [6] is often used. This transformation requires a correct choice of the transformation parameter often denoted as (λ). [4] suggested using a maximum likelihood value for the choice of the value of λ that results in the smallest residual sum of squares. [9] proposed a Bayesian method to choose the value of λ for a given model structure. The correct choice of the transformation parameter (λ), the simultaneous transformation and fitting of the model of a given series are the noticed limitations in the use of [6] power transformation. To remedy these limitations, [10] have shown how to apply Bartlett’s transformation to time series data. Accordingly, they regress the natural logarithms of the group standard deviation ) ,..., 2 , 1 , ˆ ( m i i   against the natural logarithms of the group means ( , 1,2, , ) i X i m  of time series data arranged chronologically in equal groups and determine the slope (β) of the relationship. [11] derived a confidence interval for the index of a power transformation that stabilizes the variance of a time series. They claimed that the confidence interval for the minimum coefficient of variation can also be used to construct confidence interval for any coefficient of variation. [12] used Box-Cox transformation approach to transform a streamflow time series data to turn the non-Guassian heavy tailed distribution to a nearly Gaussian series. [13] applied log-transformation in time series modeling of US macroeconomic data. He demonstrated that the claim previously made concerning improvement in forecast accuracy following bias correction for the transformed data were not generally well founded. [14] in his work applied differencing and Box-Cox transformation on a non-stationary series prior to application of ARMA model. He discovered that the technique was not optimal for forecasting of the study data. Therefore, he recommended a method for modeling time series data with unstable roots and changing parameters. The appropriateness in context of out-of-sample forecast of vector autoregressive models (VAR) when the series is log-transformed was examined by [15]. The 58 BMSA Volume 17","PeriodicalId":252632,"journal":{"name":"Bulletin of Mathematical Sciences and Applications","volume":"33 1","pages":"0"},"PeriodicalIF":0.0000,"publicationDate":"2016-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Bulletin of Mathematical Sciences and Applications","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.18052/WWW.SCIPRESS.COM/BMSA.17.57","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

Abstract

This study examines the discrete wavelet transform as a transformation technique in the analysis of non-stationary time series while comparing it with power transformation. A test for constant variance and choice of appropriate transformation is made using Bartlett’s test for constant variance while the Daubechies 4 (D4) Maximal Overlap Discrete Wavelet Transform (DWT) is used for wavelet transform. The stationarity of the transformed (power and wavelet) series is examined with Augmented Dickey-Fuller Unit Root Test (ADF). The stationary series is modeled with Autoregressive Moving Average (ARMA) Model technique. The model precision in terms of goodness of fit is ascertained using information criteria (AIC, BIC and SBC) while the forecast performance is evaluated with RMSE, MAD, and MAPE. The study data are the Nigeria Exchange Rate (2004-2014) and the Nigeria External Reserve (1995-2010). The results of the analysis show that the power transformed series of the exchange rate data admits a random walk (ARIMA (0, 1, 0)) model while its wavelet equivalent is adequately fitted to ARIMA (1,1,0). Similarly, the power transformed version of the External Reserve is adequately fitted to ARIMA (3, 1, 0) while its wavelet transform equivalent is adequately fitted to ARIMA (0, 1, 3). In terms of model precision (goodness of fit), the model for the power transformed series is found to have better fit for exchange rate data while model for wavelet transformed series is found to have better fit for external reserve data. In forecast performance, the model for wavelet transformed series outperformed the model for power transformed series. Therefore, we recommend that wavelet transform be used when time series data is non-stationary in variance and our interest is majorly on forecast. 1.0 Introduction In several organizations, managerial decisions are largely based on the available information of the past and present observations and possibly on the process that generate such observations. A time series data provides such information. Time series is used to represent the characterized time course of behavior of wide range of several systems which could be biological, physical or economical. The utility of the time series data lies in the result of the time series analysis. Such analysis will be helpful in achieving the aim for collection of such data which could be for description (exposing the main properties of a series), explanation (revealing the relationship between variables of a series especially when observations are taken on two or more variables), forecasting (prediction of the future values of a series) and control (taking appropriate corrective actions) [1]. To analyze any time series data, time series analysis techniques are adopted. The commonly used techniques are: descriptive technique, probability models technique and spectral density analysis technique. The inference based on the descriptive method and probability models is often referred to as analysis in time domain while inference based on spectral density function is referred to analysis in frequency domain [2,3,4,5]. Bulletin of Mathematical Sciences and Applications Submitted: 2016-09-12 ISSN: 2278-9634, Vol. 17, pp 57-74 Revised: 2016-09-25 doi:10.18052/www.scipress.com/BMSA.17.57 Accepted: 2016-10-24 2016 SciPress Ltd, Switzerland Online: 2016-11-01 SciPress applies the CC-BY 4.0 license to works we publish: https://creativecommons.org/licenses/by/4.0/ All these models assume that the error component of the series (et) is normally distributed with zero mean and constant variance   2  or that the series is normally distributed with constant mean (μ) and constant variance   2  . When any study data violates any or all of these assumptions, the series is subjected to transformation. Transformation helps to (i) stabilize the variance of a series, (ii) make the seasonal effect when present additive and (iii) make the data normally distributed [1]. One of the transformations commonly used is the power transformation developed by [6]. [7] noted that the power transformation (i) changes the scale of the original series, (ii) may introduce bias in the forecast especially when data have to be transformed back to its original scale and (iii) often the transformed series have no physical interpretation. [1] argued that transformation alone may not be helpful when variance changes through time in the absence of trend. In such case (i.e., when variance changes through time in the absence of trend), he recommended that a model that allows for changes in variance should be considered. Wavelet method is one such method that allows for changes in variance which has been found to be useful in time series analysis. It involves decomposition, de-nosing and reconstruction of series. Decomposition involves breaking down time series into two main components namely the detail and the smooth components, de-noising deals with the removal of the non significant components of the series, while reconstruction involves recovering of the original series devoid of noise. Unlike the power transformation, wavelet method does not change the scale of the series, poses no problem in its interpretation and does not rely on the assumption of any underlying distribution of the study data. Additionally, wavelet transformation method allow for decomposition of a series without knowing the underlying functional form of the series [8]. Could this lead to an improved model and forecast performance from those based on power transformation? This and other related questions are what this study intends to address. Therefore, the objective of this study is to examine the precision (in terms of goodness-of-fit) and forecast performances of the models for wavelet transformed series while comparing it with models for power transformed time series. 2.0 Literature Review Various transformations exist, but the power transformation developed by [6] is often used. This transformation requires a correct choice of the transformation parameter often denoted as (λ). [4] suggested using a maximum likelihood value for the choice of the value of λ that results in the smallest residual sum of squares. [9] proposed a Bayesian method to choose the value of λ for a given model structure. The correct choice of the transformation parameter (λ), the simultaneous transformation and fitting of the model of a given series are the noticed limitations in the use of [6] power transformation. To remedy these limitations, [10] have shown how to apply Bartlett’s transformation to time series data. Accordingly, they regress the natural logarithms of the group standard deviation ) ,..., 2 , 1 , ˆ ( m i i   against the natural logarithms of the group means ( , 1,2, , ) i X i m  of time series data arranged chronologically in equal groups and determine the slope (β) of the relationship. [11] derived a confidence interval for the index of a power transformation that stabilizes the variance of a time series. They claimed that the confidence interval for the minimum coefficient of variation can also be used to construct confidence interval for any coefficient of variation. [12] used Box-Cox transformation approach to transform a streamflow time series data to turn the non-Guassian heavy tailed distribution to a nearly Gaussian series. [13] applied log-transformation in time series modeling of US macroeconomic data. He demonstrated that the claim previously made concerning improvement in forecast accuracy following bias correction for the transformed data were not generally well founded. [14] in his work applied differencing and Box-Cox transformation on a non-stationary series prior to application of ARMA model. He discovered that the technique was not optimal for forecasting of the study data. Therefore, he recommended a method for modeling time series data with unstable roots and changing parameters. The appropriateness in context of out-of-sample forecast of vector autoregressive models (VAR) when the series is log-transformed was examined by [15]. The 58 BMSA Volume 17
小波变换在时间序列分析中的替代作用
本文研究了离散小波变换作为一种变换技术在非平稳时间序列分析中的应用,并将其与幂变换进行了比较。用巴特利特恒方差检验进行恒方差检验和适当变换的选择,用Daubechies 4 (D4)最大重叠离散小波变换(DWT)进行小波变换。利用增广Dickey-Fuller单位根检验(ADF)检验了变换后的幂级数和小波级数的平稳性。采用自回归移动平均(ARMA)模型技术对平稳序列进行建模。使用信息标准(AIC、BIC和SBC)确定拟合优度方面的模型精度,同时使用RMSE、MAD和MAPE评估预测性能。研究数据是尼日利亚汇率(2004-2014)和尼日利亚外汇储备(1995-2010)。分析结果表明,汇率数据的幂变换序列允许随机游走(ARIMA(0,1,0))模型,而其小波当量则充分拟合到ARIMA(1,1,0)。同样,外部储备的幂变换版本充分拟合ARIMA(3,1,0),其小波变换当量充分拟合ARIMA(0,1,3)。在模型精度(拟合好度)方面,幂变换系列的模型更适合汇率数据,而小波变换系列的模型更适合外部储备数据。在预测性能上,小波变换序列模型优于幂变换序列模型。因此,我们建议在时间序列数据方差是非平稳且我们主要关注预测的情况下使用小波变换。在一些组织中,管理决策在很大程度上是基于过去和现在的观察所得的信息,可能还基于产生这些观察的过程。时间序列数据提供了这样的信息。时间序列用来表示生物、物理或经济等多种系统行为的特征时间过程。时间序列数据的效用在于时间序列分析的结果。这种分析将有助于实现收集这些数据的目的,这些数据可以用于描述(揭示一系列的主要属性),解释(揭示一系列变量之间的关系,特别是当对两个或多个变量进行观察时),预测(预测一系列的未来值)和控制(采取适当的纠正措施)[1]。为了分析任何时间序列数据,都采用时间序列分析技术。常用的技术有:描述技术、概率模型技术和谱密度分析技术。基于描述方法和概率模型的推理通常称为时域分析,而基于谱密度函数的推理通常称为频域分析[2,3,4,5]。数学科学与应用通报提交日期:2016-09-12 ISSN: 278-9634 Vol. 17 pp . 57-74修回日期:2016-09-25 doi:10.18052/www.scipress.com/BMSA.17.57接受日期:2016-10-24 2016 SciPress Ltd, Switzerland在线日期:2016-11-01 SciPress对我们发表的作品采用CC-BY 4.0许可协议:https://creativecommons.org/licenses/by/4.0/这些模型都假定该序列的误差分量(et)是均值为零、方差为常数的正态分布2或均值(μ)和方差为常数的正态分布2。当任何研究数据违反任何或所有这些假设时,该系列将受到转换。变换有助于(i)稳定序列的方差,(ii)使季节效应存在加性时,(iii)使数据正态分布[1]。常用的变换之一是b[6]开发的功率变换。[7]指出,功率转换(i)改变了原始序列的尺度,(ii)可能在预测中引入偏差,特别是当数据必须转换回原始尺度时,(iii)转换后的序列通常没有物理解释。[1]认为,在没有趋势的情况下,当方差随时间变化时,仅靠转换可能没有帮助。在这种情况下(即,当方差随时间变化而没有趋势时),他建议应该考虑允许方差变化的模型。小波方法是一种允许方差变化的方法,在时间序列分析中被发现是有用的。它涉及到序列的分解、去噪和重构。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
自引率
0.00%
发文量
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信