2.1. Approximating the Taylor Residuals via the Fourier Function for the ESTAR (1) Unit Root Test
This section will show that the Fourier function can be used to increase the power of an ESTAR type unit root test obtained by computing the Taylor expansion. Ref. [
11] mentioned the importance of the transition speed parameter while producing the ESTAR type unit root test. When the transition speed parameter is zero, there is no nonlinearity, and the unit root problem arises. In order to understand the situation better, it is helpful to summarize the study of [
11].
where
,
and
are unknown parameters.
is a stochastic process with zero-mean. Ref. [
11] uses the ESTAR function as a transition function.
In Equation (2), is the transition variable where and denote the speed of the mean reversion and the lag of the transition variable, respectively. Moreover, is a bounded function between 0 and 1.
We can define the model by inserting Equation (2) into Equation (1):
When we rewrite Equation (3) in terms of the traditional unit root testing methodology, we obtain Equation (4) as follows:
where
and
.
In Equation (3), if
is away from the attractor, it can converge to the attractor (i.e., mean of the nonlinear process). On the one hand, if
is sufficiently away from the attractor, it can converge to the mean of the nonlinear system. On the other hand, if
is not sufficiently away from the attractor, it exhibits some instabilities in the neighborhood of the attractor. Such behavior is also prevalent in the asset markets. If the difference between the risk-adjusted returns of two assets is large, then the profitability of arbitrage is higher than the case where the mentioned difference is small due to the transaction costs. Consequently, the speed of the mean reversion to the equilibrium changes depending on the size of this differential inversely proportionally. The second empirical application carried out by [
11] was motivated by Sercu et al., (1995) [
20] and Michael et al. (1997) [
21], which examine the nonlinearity included in the purchasing power parity (PPP) hypothesis. These studies point out that, the larger the deviation from the PPP, the faster the mean-reverting behavior to the equilibrium. Again, the speed of the mean reversion is determined by the distance to the attractor. Here, in the framework of Equation (4),
; that is, the outer regimes are convergent
and inner regime is divergent (
), finally
. It states that the proposed nonlinear stochastic process is locally unit root (explosive) globally stationary. In these circumstances, the transition variable stays in the middle regime for values of
whereas the transition variable proceeds to the outer regime when it has bigger values. Therefore, it has stable dynamics and so, it is geometrically ergodic. Ref. [
11] shows that these conditions are satisfied. This nonlinear behavior can also be seen in
Figure 1. In the null hypothesis, which indicates the unit root case,
and
whereas in the alternative hypothesis,
and
. In this case,
. So, the nonlinear ESTAR stationarity is satisfied when
and
.
In view of these circumstances,
and we follow [
11] by considering the case where
. Then the model can be expressed in the following form:
Since
, we focus on the value of
to test the unit root. Thus, the null hypothesis is as follows:
and the alternative hypothesis is:
We cannot directly test the unit root hypotheses because there is no
under the null hypothesis. Therefore, following Davies (1987) [
22] and Luukkonen, Saikkonen and Teräsvirta (1988) [
23], we apply the Taylor expansion to
under the null hypothesis.
We obtain the auxiliary regression as in [
11] when we take
.
After some algebraic operations, we obtain the following auxiliary regression:
where the null hypothesis is,
and the alternative hypothesis is
Then, the
-type test statistics can be calculated as follows:
In this study, we use the approximation of the nonlinear process instead of its actual values, where there is an information loss. The method that compensates the information loss between the approximated and actual nonlinear processes can provide testing procedures as effective as the procedures using the actual nonlinear process. This issue has already been pointed out in the literature, yet it is misinterpreted. To show the method proposed in this study, we prefer to show the residual structure of Taylor expansion.
In the above given equation,
represents the Taylor expansion residuals that are considered as either a modeling error or an omitted functional form (such as omitted variable bias or wrong functional form) in the literature. The simplest way to restore this
is with the help of Fourier functions. Fourier functions can converge to functions of unknown structures with ([
16]). The Fourier function along with the ESTAR test was first used in [
15] who consider the Fourier function as a smooth transition of multiple structural breaks in the FKSS test, as it is done in [
17,
18]. However, using simulation studies, [
16] show that the power of the integer Fourier function–ESTAR test is smaller than the integer Fourier function test. Following [
16], in this study, we will show that the use of fractional Fourier instead of integer one can improve this convergence.
If the Taylor expansion can be performed as in the method given in Equation (15), the KSS-FF test proposed in this study will be considered as a state-dependent nonlinear unit root test, rather than being a test combining state-dependent nonlinearity with a structural break.
To clarify this point, we run the following simulation.
- Step 1.
ESTAR (1) data are generated. Then, they are regressed on the series obtained with ESTAR (1) Taylor expansion.
- Step 2.
Use the same ESTAR (1) data in Step 1. Then, they are regressed on the series with augmenting fractional frequency Fourier function. Obtain the estimated fractional frequency.
- Step 3.
Residuals obtained from these two regressions are subtracted from originally obtained residuals and the differenced series thereby obtained.
- Step 4.
The distance between the original series and the cumulative summation of squared differenced residuals is calculated.
- Step 5.
In the case where the difference between KSS-FF residuals and original model error terms is less than the difference between KSS residuals and original model error terms, Equation (12) becomes valid.
As it is displayed in
Table 1 and
Figure 2, the KSS-FF model reintroduces Taylor residuals into the model in a way that satisfies Equation (12). It is remarkable that the results in which we are interested are coherent with the power analysis applied to ESTAR (1) DGP by [
16]. Here, the power of the KSS test is better than that of all the other tests previously used in the literature, naturally. Nonetheless, the power of the KSS-FF test is the best in almost all parameter regions, while the second-best is KSS test and the third-best is integer frequency Fourier unit root test (FADF) used in the study of [
18]. Ref. [
17] claimed that the KSS test is a rival to their FADF test, but they did not compare the power of the FADF test with that of the KSS test. In this sense, the FADF test is a competitor to KSS, which is also confirmed by [
16]. Nevertheless, we have shown that the KSS-FF test, which is proposed in this study, really increases the power by taking the Taylor expansion term in the KSS test, according to the simulation. Another point to emphasize is that the KSS-FF test proposed is not a test for detecting structural breaks, but a state-dependent nonlinear unit root test. The KSS-FF test is one type of the ESTAR (1) test, which exhibits information loss due to the Taylor expansion, but which regains the lost information by the fractional Fourier function. Accordingly, it would be a more appropriate terminology to describe the KSS-FF test as an ESTAR test with augmented fractional frequency Fourier trend. As it is displayed in
Figure 3, this can be seen more clearly by drawing the transition function including the Taylor residuals to the model with the Fourier function. For this purpose, we use the generated data and plot them with the transition functions. To do so, we take the
and
Fourier trends to the 0–1 range to bring them to the same scale.
2.2. ESTAR (1) Fractional Frequency Fourier Tests and Their Asymptotic Distributions
Let us consider the following data-generating process:
where
and the initial condition
is zero.
is a de-meaned and de-trended series with fractional frequency Fourier function. Equation (18) assumes that the adjustment speed is nonlinear and follows an ESTAR process developed in [
11]. In this study, we propose the two-step testing procedure of [
15]. In the first step, we obtain the de-meaned or de-trended series,
as follows:
where
,
,
and
are OLS estimators for de-meaned and de-trended cases, respectively. Next, we build the ESTAR (1) augmented fractional frequency Fourier function unit root test by using the de-meaned and de-trended series
in the second step. Then, by applying the steps from Equation (1) until Equation (11), the following test equation is obtained:
where
with
representing the Taylor expansion residuals and
. If residuals,
, are serially correlated, then Equation (21) can be augmented to correct the serial correlation as follows:
We can collect the t–statistics for
against
as follows:
where
is the OLS estimator of
and
is the standard error of
.
To obtain the asymptotic distribution of the test, we need the subsequent outcomes, where we let , , be an integer close to . During the derivation implies weak convergence as approaches .
Theorem 1. Under the null hypothesis of unit root, the asymptotic distribution of is as follows:where is the Wiener process as defined over the interval . for is a demeaned and de-trended Brownian motion, respectively. Apparently, the asymptotic distribution of the subsequent test statistics under the null depends on the fractional Fourier frequency , however invariant to other parameters in the model. Therefore, the importance of frequency estimation in ESTAR fractional frequency (1) test is shown explicitly by deriving the asymptotic distribution.
Asymptotic critical values of the KSS-FF statistics for the intercept only and the intercept and trend cases have been arranged via simulations with T = 2000 and for 100000 replications. Critical values of the mentioned cases are displayed in
Table 2 and
Table 3 (For critical values for the one step intercept only and the one step intercept and trend cases, see
Appendix B).