Consistency of the Model Order Change-Point Estimator for GARCH Models

GARCH models have been commonly used to capture volatility dynamics in financial time series. A key assumption utilized is that the series is stationary as this allows for model identifiability. This however violates the volatility clustering property exhibited by financial returns series. Existing methods attribute this phenomenon to parameter change. However, the assumption of fixed model order is too restrictive for long time series. This paper proposes a change-point estimator based on Manhattan distance. The estimator is applicable to GARCH model order change-point detection. Procedures are based on the sample autocorrelation function of squared series. The asymptotic consistency of the estimator is proven theoretically.


Introduction
Modelling volatility of financial asset returns is particularly an important area in Finance.This is because volatility is considered to be a measure of risk when pricing financial instruments.The series particularly is characterized by the property of volatility clustering and thus can be considered to display a stationary behavior for some time then suddenly the variability changes, it stays constant for some time at this new value until another change occurs.This therefore suggests that the financial returns series is non-stationary and can be looked at as a union of several stationary series.GARCH models have been commonly used to capture volatility dynamics in financial time series particularly in model-ling of stock market volatility as seen in [1] [2] [3] and derivative market volatility as utilized by [4] [5] and [6].
Given the changing pace of the underlying economic mechanism and technological progress, modeling economic processes over a long time horizon, it is possible that structural changes may occur.This can cause the time series to deviate from stationarity and result to volatility clustering.The detection of these structural change points is therefore vital to various players in a given economy to ensure timeliness of decisions.A fundamental problem in financial trading is the correct and timely identification of turning points in stock value series.This detection enables one to make profitable investment decisions, such as buyingat-low and selling-at-high, hence traders require early identification of local troughs and peaks of stock values.In macroeconomics, knowing the beginning of a recession leads to an increase of government expenditure or an expansion of money supply.
A key assumption of the GARCH models used is that the process is stationary as this allows for model identifiability.However, this violates the volatility clustering property exhibited by the financial returns series.This phenomenon is manifested by the fact that the absolute value of returns or their squares display a positive, significant and slowly decaying autocorrelation function despite the fact that the returns are uncorrelated.This indicates that modeling financial returns series over long time horizons deviates from the stationarity assumption suggesting the existence of a change-point in the series.A modification of the GARCH model, specifically the IGARCH model, has been proposed to model the persistent changes in volatility as the stationarity assumption is relaxed.
However the IGARCH model is prone to some shortcomings.[7] showed that the behavior of an IGARCH process depends on the intercept, such that, if the intercept is positive then the unconditional variance of the process grows linearly with time.In practice this means that the amplitude of the clusters of volatility to be parametrized by the model on the average increases over time.The rate of increase needs not, however, be particularly rapid.If the intercept is zero in the IGARCH model, the realizations from the process collapse to zero almost surely.
However, a potentially disturbing fact is that the model assumes that the unconditional variance of the process to be modeled does not exist in that the variance may be infinite [8] and [9].
It is argued that in applications, the assumption of parameter constancy in GARCH models may not be appropriate especially when the series to be modeled are long [10].To overcome this problem of modeling financial time series in the presence of structural changes, the duo suggests that one option is to assume that the parameters change at specific points of time, divide the series into sub-series according to the location of the change-points and fit separate GARCH models to the sub-series.This brings about the challenge of determining the number of change-points and their location because they are normally not known in advance.This proposition has been adopted by various researchers who have utilized different methodologies to be able to locate the change-points attributed to change in parameter specification.The use of squared model residuals and likelihood ratio to detect parameter changes is proposed by [11] while [12] proposes the use of Markov-switching GARCH models estimated through Markov Chain Monte Carlo simulation methods.Modeling of equity volatilities as a combination of macroeconomic effects and time series dynamics by combining exponential splines and GARCH models is utilized by [13].An alternative approach is to use smooth transition GARCH model.This can be achieved by defining a transition function where the coefficients are expressed as a function of time as in [9] or by employing non-linear functions that are lagged for the squared observations [14] and [15] or lagging the conditional variance [16].
CUSUM tests have also been proposed as suitable methods of identifying change-points by establishing breaks in moments of the time series.The use of the unconditional variance is proposed by [17] while [18] [19] use the mean.
However, these methods are mainly based on the assumption that change-points occur solely due to change in parameter specification.The approach presented here seeks to identify change-points attributed to change in model order specification.
This paper is organized as follows.Section 2 gives an overview of GARCH model specification with corresponding assumptions utilized in the proof of the main result.Section 3 presents the proposed change-point estimator for change attributed to in the model orders p and q in GARCH models.The estimator is based on the Manhattan distance of sample autocorrelations of a squared returns series.Section 4 provides proof of the consistency of the proposed change-point estimator.

GARCH Model
Assume that the data t X , for t ∈  , are independent and sampled at equi- spaced points.The series ( ) t t X ∈ describe a financial returns time series mod- eled using ( )

GARCH ,
p q model specified as: The sequence of innovations ( )


is an independent and identically distributed (iid) sequence with mean zero and unit variance.( ) σ ∈ is the volatil- ity sequence of the GARCH model.Assume that 0 p q α β ≠ and that all coeffi- cients i α and j β are non-negative to avoid ambiguity with regards to orders ( ) , p q .Since we are not interested in the trivial solution 0 t X ≡ to (1), further assume that 0 0 α > .
Let p q = and is a sequence of independent and identically distributed random variables such that , i t c is independent of t σ .This allows us to rewrite (1) as ( ) Model ( 2) is utilized in the proof of consistency of the proposed change-point estimator.
The GARCH (p,q) model ( 1) can also be represented as an ARMA(max(p,q),q) as showed by [20] as where and ( ) u ∈ is white noise.
This representation of the GARCH model follows the standard ARMA form for the squared series, therefore, conventional methods used to identify ARMA processes may be used to determine the presence of GARCH.Of keen interest is the use of the sample autocorrelation function (SACF) and partial autocorrelation functions (PACF).Specifically, the orders p and q are drawn from the autocorrelation function and partial autocorrelation function respectively.Empirically, these orders are chosen such that the SACF cuts off after lag p and the PACF decline exponentially to zero after lag q for which they are significant.In light of this it can be asserted that the SACF and PACF can be used to distinguish GARCH model with different model orders specifications.
The following assumptions are necessary to prove the subsequent theoretical results.According to [21] the existence of a unique strictly stationary solution to (1) is the negativity of the top Lyapunov exponent.This however cannot be calculated explicitly but a sufficient condition for this is given by Assumption 3 (Ergodic Process) According to [22] standard ergodic theory yields that ( ) t X is an ergodic process.Thus its properties can be deduced from a single sufficiently large random sample of the sample.

Change-Point Estimator
Assume that the data { } 1 n t t X = describes a financial returns time series modeled using ( ) GARCH , p q process.A single change-point testing problem is first considered where it is assumed that a change-point can happen only at time k where 1 1 k n < < − .The hypotheses to be investigated are assumed to follow the following definition: where , \ 0 GARCH , for 1, , : Let I =  be a finite index sequence and ( ) ( ) The autocovariance and autocorrelation functions can be expressed in terms of the inner product as ( ) sd Y represents the standard deviation of X and Y respectively which represents an 2 L distance from the mean.
By the Assumption 3 that the series ( ) X ∈ is ergodic, then it is implied that the sample moments converge in probability to the population moments.It therefore follows that the sample autocovariance and autocorrelation converge in probability to the population autocovariance and autocorrelation respectively.
Theorem 1. (Holder's Inequality) Let I be a finite or countable index set.Given , where ( ) ( ) in the Holders Inequality Theorem 1 we obtain Thus, applying the result in ( 9) to ( 6) and ( 7) yields Following (11) define sequences of autocorrelation functions 1, i j ρ + where for fixed 0 i = , 1 1 j n ≤ ≤ − and for fixed j n = , 1 1 i n ≤ ≤ − to be such that we Journal of Mathematical Finance have two subsequences ( ) ( )  ) where for 0 , 0 Specifically, assume the case when 1 p = in (12) resulting into a weighted Manhattan distance and by linearity and absolute value of inequalities of the expectation operator results into ( ) To facilitate the construction of the proposed estimator the lower bound of the divergence measure ( 13) is assumed.Further assume that the autocorrelation function is calculated at lag : 0 h h n < < .The proposed change-point estimator is thus developed from the process generated by this measure as follows: ( ) ( ) From ( 14) it can be seen that the proposed test is a weighted difference between the sample autocorrelation functions 1,k ρ and The weight k w is a measurable function that depends on the sample size n and change-point .It is arbitrarily chosen such that it satisfies the condition that Equating ( 14) and ( 15) determines the weight k w as follows: The resultant process is obtained from ( 14) and ( 16) and defined as The change-point estimator k of a change point * k is the point at which there is maximal sample evidence for a break in the sample autocorrelation function of the squared returns process.It is therefore estimated as the least value of k that maximizes the value of

Simulation Study
The performance of the proposed estimator is examined by considering the effects of the change in sample size.Assume that { } The following table gives the parameter estimates used in the simulation. in Table 1 for changes in order q.
The results for the change in order q in Table 1 show that as the sample size increases, the similarity index given by ARI generally increases.

Consistency of the Change-Point Estimator
Consider a sample , , , n X X X  satisfying ( 2) and ( 5) and the change-point estimator k given by (18).If the sequences { } then for where C is a positive constant.
Proof Suppose that { }

GARCH , p q
sequences as defined in model (2).Further suppose that a sample where * k is the unknown change point.More specifically assume that the two sequences have different model order specification such that  ( ) Equation (24) shows is that ( )

p p p t i t i t i i t i t i l t l m t m t l t m
To establish ( ) is determined using the following Theorems as proved by [23].
The expected value of the sample autocorrelation function, ( ) k E φ , is first evaluated using (29) and (30).
Further assuming that ( 22) and ( 23) are satisfied, evaluate (31) for 1  1 0 0 For GARCH (1,1)  k k Equation (34) shows that, in the presence of a change-point, the expected value of the sample autocorrelation function before and after the true change-point * k is not equal.We consider a special case of change from GARCH(1,1) to GARCH(2,2) where we evaluate (31) for 2 2 2 p q = = as follows: 1 0 0 From ( 34) and (37), it can be seen that in the presence of a change-point then To proof consistency of the estimator we need to show that as the sample size n increases ( ) ( ) From ( 39) and (40) it follows that We also have that Thus from (41) and (50) as well as replacing τ with τ in (41) we have that Consider k n D as given in (17), the estimate ( ) ( ) , 1 1 As n → ∞ from (48) we can see that using a divergence measure motivated by the weighted p L distance, with k denoting the change-point.For 0 p > define

(
. The single change-point estimation problem is considered where the change-point k is fixed at the plots for the location of the change-point estimator as estimated by the proposed estimator(18) for various sample sizes.The hypothesis considered here is when change occurs in model order q, described as; 500, 1000 and 2000 respectively as displayed in Figures1-3.The performance of the estimators is evaluated using the Adjusted Rand Index (ARI) which compares the segmentation created by the change-point estimator and the true segmentation.The Adjusted Rand Index lies between 0 and 1.When the two partitions agree perfectly, the ARI is 1.The results of the ARI are provided I. W. Irungu et al.DOI: 10.4236/jmf.2018.82018273 Journal of Mathematical Finance

Figure 1 .Figure 2 .Figure 3 .
Figure 1.Single Change-Point for Stationary Series GARCH t X for 500.n = then in the presence of the change-point, the sequence { } 2 k X is no longer stationary.The foundation of this proof is based on the second and fourth moments of { } t X which will first be derived.Assume that the GARCH(p,p) model (2) has a finite fourth moment and let

Table 1 .
Adjusted rand index given changes in order q.