Next Article in Journal
A Complete Characterization of Bipartite Graphs with Given Diameter in Terms of the Inverse Sum Indeg Index
Next Article in Special Issue
Representative Points Based on Power Exponential Kernel Discrepancy
Previous Article in Journal / Special Issue
Confidence Intervals for the Ratio of Variances of Delta-Gamma Distributions with Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Unified Test for the AR Error Structure of an Autoregressive Model

1
School of Statistics, Jiangxi University of Finance and Economics, Nanchang 330013, China
2
Key Laboratory of Data Science in Finance and Economics, Jiangxi University of Finance and Economics, Nanchang 330013, China
*
Author to whom correspondence should be addressed.
Submission received: 25 October 2022 / Revised: 20 November 2022 / Accepted: 28 November 2022 / Published: 1 December 2022
(This article belongs to the Special Issue Computational Statistics & Data Analysis)

Abstract

:
A direct application of autoregressive (AR) models with independent and identically distributed (iid) errors is sometimes inadequate to fit the time series data well. A natural alternative is further to assume the model errors following an AR process, whose structure however has essential impacts on the statistical inferences related to the autoregressive models. In this paper, we construct a new unified test for checking the AR error structure based on the empirical likelihood method. The proposed test is desirable because its limit distribution is always chi-squared regardless of whether the autoregressive model is stationary or non-stationary, with or without an intercept term. Some simulations are also provided to illustrate the finite sample performance of this test. Finally, we apply the proposed test to a financial real data set.
MSC:
62M10; 91B84; 01-08; 65C60

1. Introduction

When auxiliary variables are not available, autoregressive models are widely used to model this kind of time series data. Typically, the response is often assumed to depend linearly on its previous values. Among all autoregressive models, the autoregressive model of the first order, i.e., AR(1), is the simplest, which takes the following form:
X t = μ + ϕ X t 1 + ε t , t = 1 , 2 , , n ,
where μ and ϕ are unknown parameters with μ being the intercept item and ϕ the autoregression coefficient, and { ε t } denotes the sequence of random errors or innovations having means of zero.
In many previous studies, a considerable amount of work has been provided on statistical inferences [1,2,3,4,5,6] and related applications [7,8,9] for AR models. In terms of practical applications, AR models are commonly used to describe the behavior of inflation or logarithmic exchange rate, where people are interested in whether there is a unit root or persistence of related variables. However, a precondition for an accurate unit root test or persistence test is that the model is properly fitted so that the parameters can be reasonably estimated. To guarantee this, it is important to perform predefined tests, e.g., the unit root test and serial correlation test, on the rationality of using the AR model before conducting a relevant economic analysis.
Among them, the unit root test is the most commonly mentioned. Note that the limit distributions of the estimators of μ and ϕ depend on whether the process { X t } is stationary or non-stationary, i.e., Case (i) | ϕ | < 1 (stationary), Case (ii) μ = 0 and ϕ = 1 + c n for some nonzero constant c (nearly integrated if c 0 , and unit root if c = 0 ), and Case (iii) μ 0 and ϕ = 1 + c n for some nonzero constant c (nearly integrated if c 0 ). It is well known that when the AR process has a unit root, its many statistical procedures have quite complex limit distributions, differing from that for the stationary case. Hence, various testing methods have been developed in the past decades to address the issue of unit root, including the augmented Dickey–Fuller (ADF) test [10], the Phillips–Perron (PP) test [11], the DF–GLS test [3], and the KPSS test [12], etc.
It is worth mentioning that if the true underlying innovations are correlated, the finite sample performance of the tests above may be greatly affected. To improve the efficiency of the estimation, a natural idea is to take into account the special structure of the errors if available. Note that it is common to assume that the errors further follow an AR process when they are correlated, while the performance of some testing procedures can be greatly improved once the AR structure has been addressed sufficiently, as shown in [13,14].
In detail, Ref. [13] considered the following autoregressive model with AR errors:
X t = μ + ϕ X t 1 + ε t , e t = ε t + i = 1 p ψ i ε t i .
where ψ = ( ψ 1 , ψ 2 , , ψ p ) denotes the vector of unknown parameters involved in the AR errors, and e t denotes the random error involved in ε t . Compared to Model (1), Ref. [13] here further assumed that ε t follows an AR process. Note that (2) implies ε t = i = 1 p ( ψ i ) ε t i + e t . A unified unit root test was developed by considering the special structure in { e t } . Their test has been shown to have desirable properties, as the related statistic converges in distribution to a standard chi-squared distributed variable. However, their test depends on preconditions such that the AR structure of { ε t } has been well specified, and p is properly predefined. The violation of these conditions may result in power loss in this method, as shown in our simulations.
To this end, we are interested in producing statistics to test whether ψ is equal to some given constant vector ψ 0 under Cases (i)–(iii), which has not been considered in the literature to the best of our knowledge. Note that although many tests have been developed for testing the possible serial correlation in { ε t } , including the Lagrangian multiplier (LM) test [15], Box–Pierce (BP) test [16], and Ljung–Box (LB) test [17], etc., they cannot be used directly to test the hypothesis above. In view of this, we propose an empirical likelihood-based statistic for testing this issue by taking into account the AR structure. Note that the setting in Case (ii) causes issues in the derivation of the asymptotic distribution, as well as the related applications. A new data-splitting idea is also employed in order to unify Cases (i)–(iii). It turns out that the proposed statistic converges in distribution to a standard chi-squared distributed variable regardless of { X t } being stationary or non-stationary, due to the special block structure of the asymptotic covariance matrix. The simulations show that our method has a good size, as well as nontrivial power performance in finite sample cases.
As a nonparametric method, empirical likelihood (EL) was firstly proposed by [18]. Because of its many excellent properties, i.e., no need to assume the parameter distribution in advance, it has been widely used in the literature when parametric methods do not work well to produce satisfactory results. Many authors have devoted themselves to extending this method. To name but a few, Ref. [19] obtained confidence regions for vector-valued statistical functions, which is a multivariate generalization of the work of [18]. Refs. [20,21] extended the empirical likelihood method to the setting of regression models and general estimation equations, respectively. Recently, Ref. [22] discussed the possibility of constructing unified tests by using empirical likelihood based on a weighted technique for time series models. Ref. [13] further extended this weighted technique to AR models with AR errors. Further, Ref. [23] applied the empirical likelihood method to test the heteroscedasticity for errors of single-index model. Ref. [6] developed a unified empirical likelihood inference method to test the predictability regardless of the properties of the predicting variable. Ref. [24] considered the unified test problem in a predictive regression model. To move the effect of the possible existence of an intercept, the idea of data-splitting has also been developed in [24]. The literature above inspired the current research.
We organize the rest of this paper as follows. Section 2 develops the unified test for the AR structure of the AR models. Section 3 reports the finite-sample simulation results. Section 4 applies the proposed test to the exchange rates between the U.S. dollar and eight countries. Section 5 concludes this paper. The detailed proof of the main theorem is specified in Appendix A.

2. Methodologies and Asymptotic Results

Supposing the random observations { X t } t = 1 n are generated from the model (2) with possible AR errors. Formulate ψ = ( ψ 1 , ψ 2 ) and let ψ 0 = ( ψ 1 , 0 , ψ 2 , 0 ) be its true value.
Note that when θ = θ 0 , { e t ( θ ) } is a sequence of iid variables, it is more efficient to construct a statistical procedure on t = p + 1 n e t 2 ( θ ) than on t = p + 1 n ( X t μ ϕ X t 1 ) 2 , as discussed in [13], where
e t ( θ ) = X t μ ϕ X t 1 + i = 1 p ψ i ( X t i μ ϕ X t i 1 ) , and ε t ( θ ) = X t μ ϕ X t 1 ,
for a given θ = ( μ , ϕ , ψ ) . However, their method depends on an assumption that the structure of the AR errors has been correctly specified, which needs to be pretested in practice. This motivates us to consider the following hypothesis:
H 0 : ψ 2 = ψ 2 , 0 v e r s u s H 1 : ψ 2 ψ 2 , 0 .
Remarkably, when ψ 2 , 0 = 0 , { ε t } is a sequence of iid errors.
Note that when θ takes the true underlying value θ 0 , we have
E ( Z t * ( θ ) | F t 1 ) = 0 , for t = p + 1 , , n ,
where Z t * ( θ ) = ( Z t , 1 * ( θ ) , , Z t , 2 + p * ( θ ) ) , F t denotes the sigma field generated by { e s : t t } , and
Z t , 1 * ( θ ) = e t ( θ ) Z t , 2 * ( θ ) = e t ( θ ) ( X t 1 + i = 1 p ψ i X t + m i 1 ) Z t , 2 + i * ( θ ) = e t ( θ ) ( X t i μ ϕ X t i 1 ) , i = 1 , , p ,
which can be obtained by taking the partial differential to the sum of least squares, i.e.,
t = p + 1 n X t μ ϕ X t 1 + j = 1 p ψ j ( X t j μ ϕ X t j 1 ) 2 ,
with respect to ψ . Then, similar to [21], one can use the profile empirical likelihood method to construct a test for hypothesis H 0 based on { Z t * ( θ ) } .
However, following [22], it is easy to verify that the resulting test does not converge in distribution to a standard chi-squared variable because the quantity 1 n t = p + 1 n Z t * ( θ 0 ) Z t * ( θ 0 ) does not converge in probability for Case (ii), i.e., μ = 0 and ϕ = 1 + c n for some nonzero constant c (nearly integrated if c 0 , and unit root if c = 0 ). As an improvement, one may use the weighted technique developed in [22] to construct a weighted empirical likelihood-based test. Unfortunately, the resulting testing statistic still faces a similar problem in the optimization step during the process of profiling the redundant parameters; see a similar discussion in [25].
To overcome this problem, we propose the construction of the following empirical likelihood function for θ :
L ( θ ) = sup t = p + 1 m m δ t : δ 1 0 , , δ m 0 , t = p + 1 m δ t = 1 , t = p + 1 m δ t Z t ( θ ) = 0 ,
based on the data-splitting idea, where Z t ( θ ) = ( Z t , 1 ( θ ) , , Z t , 2 + p ( θ ) ) with
Z t , 1 ( θ ) = e t ( θ ) Z t , 2 ( θ ) = e t + m ( θ ) X t + m 1 1 + X t + m p 1 2 + i = 1 p ψ i X t + m i 1 1 + X t + m p i 1 Z t , 2 + i ( θ ) = e t ( θ ) ( X t i μ ϕ X t i 1 ) , i = 1 , , p ,
where m = [ n / 2 ] with [ · ] is the floor function. That is, we use the second half of the data to handle ϕ , and the first half of the data to handle the rest of the parameters. Here, 1 + X t + m p i 1 is mainly used for technical consideration, which can relieve the correlation among { Z t ( θ ) } , and consequently improve the finite sample performance of the EL test.
Since our aim is to test H 0 related to ψ , we are only interested in the parameter ψ . To this end, we treat the other parameters as redundant parameters, as in [21], and obtain the profile empirical likelihood ratio as p ( ψ ) :   = min μ , ϕ ( μ , ϕ , ψ ) .
To derive the asymptotic result for p ( ψ ) , we need the following regular conditions:
  • (C1) Suppose { X t } follows one of the following cases:
    -
    (i) (Stationary) | ϕ | < 1 , independent of n;
    -
    (ii) (Non-stationary without an intercept) ϕ = 1 c n for some constant c independent of n with μ = 0 ;
    -
    (iii) (Non-stationary with an intercept) ϕ = 1 c n for some constant c independent of n with μ 0 ;
  • (C2) ψ ( z ) = 1 j = 1 p ψ j z j 0 when | z | < 1 , and ψ ( z ) has no common root with ψ p 0 .
  • (C3) { e t } are iid random errors, and satisfy E ( | e t | 2 + δ ) < for some constant δ > 0 .
These conditions are quite common, and can be found in studies such as [13]. Here, (C2) is assumed to guarantee the stationarity of { ε t } .
Under these conditions, we have the following result.
Theorem 1.
Suppose Conditions (C1)–(C3) hold. Then, under the null hypothesis H 0 ,
p ( ψ 0 ) d χ p 2 ,
as n , where χ p 2 denotes a chi-squared random variable with p degrees of freedom, and ‘ d ’ denotes the convergence in distribution.
Remark 1.
Using a similar proof to that of Theorem 1, we can show that
˜ p ( ψ 0 , 2 ) d χ r 2 , as n ,
where ˜ p ( ψ 0 , 2 ) = min μ , ϕ , ψ 1 ( μ , ϕ , ( ψ 1 , ψ 0 , 2 ) ) with r being the dimension of ψ 0 , 2 , which is the true value.
Theorem 1 is desirable because it shows that the proposed test has a standard chi-squared distribution asymptotically, regardless of which one of the Cases (i)–(iii) is followed by { X t } . Based on Theorem 1, we may reject the null hypothesis H 0 once p ( ψ 0 ) > χ r 2 ( 1 a ) at the significance level a ( 0 , 1 ) , where χ r 2 ( 1 a ) denotes the ( 1 a ) -th quantile of χ r 2 .

3. Simulation Results

In this section, we conduct some simulations to investigate the finite sample performance of the proposed test in terms of both size and power. The simulations consist of three parts. In the first part, we investigate the finite sample performance of the proposed profile empirical likelihood, and compare it with a combination of the LB test and the Akaike information criterion (AIC), i.e., using firstly the LB test to detect whether there exists a serial correlation in the residuals, and then by employing the AIC to determine the order of the AR structure in residuals. In the second part, we investigate the possibility of using the proposed method to test whether or not ψ is equal to some given ψ 0 , which may be useful when verifying the extent of the stationarity of the AR errors. Note that the combination of the LB and AIC cannot be used to fulfill this type of task. In the last part, we study the impact of misdetermining the AR structure of the errors on the finite sample performance of the unit root test developed in [13]. The LB test is computed with the R function Box.test.R, while for the computing of the profile empirical likelihood, we first use R package emplik to obtain the log-empirical likelihood ratio, and then optimize this log ratio by using the nlm.R function. All of these R functions are well-documented, and are currently available from the CRAN of the R-project.
In the first part, the random observations { X t } are generated from the model (2) with μ { 0 , 0.01 } , which indicates that the model has no intercept and an intercept item, respectively. We take ϕ from { 0.5 , 1 , 1 1 n } , where 0.5 indicates that X t is a stationary process, and 1 indicates that it is a unit root process, while 1 1 n indicates a near unit root process. { e t } is a sequence of iid random variables with means of zero and variances of one. { ε t } follows the three different scenarios listed below.
  • S1: The null hypothesis H 0 ( 1 ) : ψ = ( 0 , 0 , 0 ) , i.e., ε t has no serial correlation. The local alternative hypothesis: ψ = ψ 0 = ( d / n , 0 , 0 ) for some d > 0 .
  • S2: The null hypothesis H 0 ( 2 ) : ψ = ( 0.1 , 0 , 0 ) , i.e., ε t has first-order serial correlation. The local alternative hypothesis: ψ = ψ 0 = ( 0.1 , d / n , 0 ) for some d > 0 .
  • S3: The null hypothesis H 0 ( 3 ) : ψ = ( 0.1 , 0.1 , 0 ) , i.e., ε t has second-order serial correlation. The local alternative hypothesis: ψ = ψ 0 = ( 0.1 , 0.1 , d / n ) for some d > 0 .
In all Scenarios S1S3, d is taken from { 1 , 3 , 5 , 7 } . All computations are carried out 10,000 times with n ranging from 300 to 1200.
Table 1 reports the size performance of the proposed method with different settings at the significance levels τ = 0.05 . We also report the ratios of determining the order of the AR error incorrectly by using the AIC of Scenarios S1S3 under the condition of H 0 for comparison. The EL method has a good performance in all Scenarios S1S3. The results show that the size values of the EL method gradually converge to the significance level as the sample size n increases, no matter whether X t is a stationary process, a near unit root process, or a unit root process, and regardless of whether μ is 0 or not. Conversely, for the AIC method, when X t follows a stationary process, the ratios of determining the order of the AR error incorrectly are only closer to 5% in S1. Note that it performs poorly for the rest of the settings, meaning that it is affected greatly by the stationarity of { X t } .
Figure 1 shows the power performance of the EL method. We can see that in S1 and S2, when X t follows a stationary process, the convergence rate is the slowest. When X t follows a near unit root process or a unit root process, as the value of d increases, the power converges quickly to 1. In S3, when X t follows a near unit root process or a unit root process, as the value of d increases to 7, the power values have a slightly descending tendency. This implies that although the stationarity of { X t } does not impact on the order of the local alternative hypothesis, it does affect the power function of the EL method.
In the second part, we consider testing whether or not ψ is equal to some given ψ 0 . We simulate two settings, i.e.,
  • (I): The null hypothesis H ˜ 0 ( 1 ) : ψ = ( ψ 1 , ψ 2 ) = ( 0.1 , 0.3 ) against the local alternative hypothesis: ψ = ψ 0 = ( 0.1 + d n , 0.3 + d n ) , for some d > 0 .
  • (II): The null hypothesis H ˜ 0 ( 2 ) : ψ 2 = 0.3 , against the local alternative hypothesis: ψ 2 = 0.3 + d n , for some d > 0 .
The other parameters are the same as those in the first part. The size ( d = 0 ) and power ( d { 3 , 5 , 10 } ) performances are shown in Figure 2. As expected, similar observations can be found in Figure 3 as in the first part of simulations.
The simulation results in the first and second parts show that the proposed EL method has a good performance in specifying the AR error structure and testing whether or not ψ is equal to some given ψ 0 , thereby confirming the theoretical result obtained in Theorem 1. It is worth noting that when taking the AR error structure into account, accurate identification is crucial, because it will affect the unit root test of the AR model. Therefore, in the third part, we conduct the following simulation to show the benefit of conducting a predefined test.
Step 1: We generate an AR model with an AR(1) error structure, and the parameters are μ = 0 , ϕ = 1 1 n , ψ = 10 n . Then, we use the EL and AIC methods to determine the order of the AR error. We consider sample sizes of 600 and 1200, repeat the tests 10,000 times, and record the order determination counts under the two methods. The results are shown in Figure 3. The abscissa represents the order of the AR error, and the ordinate is the number of each order. It can be seen from Figure 3 that under all sample sizes, the two methods show that the residuals have a serial correlation. For the EL test, in 10,000 experiments, 9398 of them are correctly ordered, and the error rate is only 6.02%. When the sample size increases to 1200, the error rate decreases to 5.84%. For the AIC, when the sample size is 600, the error rate is 34.21%, and when the sample size is 1200, the error rate is 22.75%. It is obvious that compared with the AIC method, the EL test has advantages in identifying the order of the correlated errors, which is consistent with the above simulation results.
Step 2: We use the method proposed in [13] to test the unit root of an AR (1) model when the AR error order is correctly and incorrectly determined. Table 2 records the probability of identifying a unit root when the real data are a near unit root. The results show that when the true underlying structure of the AR error is incorrectly specified, the power of the test proposed in [13] suffers from a loss compared to the case when the true underlying structure of the AR error is correctly specified. This shows the necessity of correctly testing the AR error structure before conducting the unit root test if one wants to obtain a more reliable unit root test result.
To summarize, the EL method proposed in this paper has obvious advantages in identifying the AR error structure, and these two methods are crucial in the subsequent real data analysis.

4. A Financial Real Data Application

In this section, we provide a real financial data example. The purpose of this section is to explore the error structure of different exchange rate markets. We collected the exchange rates of eight countries, including developed and developing countries, against the U.S. dollar. Currencies from developed countries include the Canadian dollar (CAD), Norwegian Kroner (NKR), Singapore dollar (SGD), Swedish Kronor (SKR) and Japanese yen (JPY). Currencies from developing countries include Chinese yuan (CNY), Thai baht (THB) and Sri Lanka Rupees (SRE). All data are downloaded from FRED database (fred.stlouisfed.org). The sample period is the daily data from 2 January 2017 to 31 December 2020 ( n = 1044 ). Their time series graphs are provide in Figure 4.
We report the least squares estimation of the unknown parameters μ and ϕ , and the testing results of the EL, LB, and AIC methods, where the LB and AIC tests were conducted on residuals obtained from the least squares method. All results are listed in Table 3, in which the second and third columns are the estimated intercept and autoregressive coefficients, respectively; the fourth column is the order determination result of the EL test; the fifth column is the p-values of the EL test; and the last two columns are the p-values of LB test and the order determination result of AIC, respectively. The AIC test shows that most sequences have a serial correlation, except for CNY and JPY, while the EL method indicates that only one country’s data has an AR error of up to an order of 2. Note that the AIC tends to determine the correlated errors with a higher order than for Cases (ii)–(iii), while for most cases, the estimated ϕ ^ is very close to 1, i.e., a near unit root. It seems that the testing results for this dataset coincide roughly with the observations in the simulations.

5. Conclusions

The AR model is widely used in time series data modeling. However, the direct application of AR models with iid errors is inadequate sometimes. A common practice is to further assume that AR models have errors of AR structure. Note that the relevant structure of the error affects the statistical inference of the AR model. Therefore, it is important to test the error structure of the model in both theoretical and practical analyses, which is not considered in the literature to the best of our knowledge. Motivated by this, this paper proposed a consistency empirical likelihood test method based on the idea of data splitting. The limit distribution of the EL statistic was proved to be chi-squared asymptotically regardless of the process { X t } being stationary or non-stationary, and with or without an intercept term. The proof is challenging and different from that of the traditional profile empirical likelihood in [21], as the quantity 1 n t = p + 1 n Z t ( θ 0 ) does not converge in distribution to a normally distributed vector. Fortunately, the limit distribution of the profile empirical likelihood-based test is still chi-squared because of the special block structure of the asymptotic covariance matrix. The simulation results illustrated that the proposed method could not only have a good performance in specifying the AR error structure, but also could sufficiently test whether the coefficients of the error item are equal to some given values, which can not be achieved by some existing serial correlation tests in the literature. The technique in the proof of Theorem 1 is challenging as in Case (ii) the theoretical proof involves handling the convergence in space, and the special block structure of the asymptotic covariance matrix. Hopefully, it is of potential usage in practice as it is difficult to detect whether the process { X t } is a unit root or near unit root process. Note that it is not necessary to make clear which case of (i)–(iii) the { X t } follows when using our proposed test in practice.
As noted by an an anonymous reviewer, an issue of interest is whether the current result can be extended to the case when ε t follows a autoregressive moving average model. Note that the discussion in this paper involves nonstationary Cases (ii) and (iii); it seems challenging to derive the related theoretical results. We will further consider this in the future.

Author Contributions

Conceptualization, X.L.; Methodology, X.L., X.W. and Y.F.; Software, X.W.; Validtion, X.L., Q.L. and L.T.; Formal analysis, X.W. and Y.F.; Investigation, Q.L. and L.T.; Resources, X.L.; Data curation, X.W.; Writing—original draft preparation, X.W., X.L., Y.F., L.T. and Q.L.; Writing—review and editing, X.L. and X.W.; Visualization, Q.L. and L.T.; Supervision, Q.L. and L.T.; Project administration, X.L.; Funding acquisition, X.L. All authors have read and agreed to the published version of the published version of the manuscript.

Funding

Xiaohui Liu’s research was funded by NSF of China grant number 11971208, National Social Science Foundation of China grant number 21&ZD152 and the Outstanding Youth Fund Project of the Science and Technology Department of Jiangxi Province grant number 20224ACB211003. Li Tan’s research was funded by Natural Science Foundation of China grant number 12061034 and China Postdoctoral Science Foundation grant number 2022M711424. Qing Liu’s research was funded by NSF of Jiangxi Province grant number 20192BAB201005, China Post doctoral Science Foundation grant number 2020M671961, and Post doctoral Science Foundation of Jiangxi Province grant number 2019KY47 and The APC was funded by Xiaohui Liu.

Data Availability Statement

The real data used in the artical can be found at fred.stlouisfed.org.

Acknowledgments

We thank the three anonymous reviewers, the Associate Editor, and the Guest Editor of this Special Issue for their insightful comments, which have led to many improvements in this paper. Xiaohui Liu’s research is supported by NSF of China (Grant No. 11971208), National Social Science Foundation of China (21&ZD152) and the Outstanding Youth Fund Project of the Science and Technology Department of Jiangxi Province (No. 20224ACB211003). Li Tan’s research is supported by Natural Science Foundation of China (No. 12061034) and China Postdoctoral Science Foundation (No. 2022M711424). Qing Liu’s research is supported by NSF of Jiangxi Province (No. 20192BAB201005), China Post doctoral Science Foundation (No. 2020M671961), and Post doctoral Science Foundation of Jiangxi Province (No. 2019KY47).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of the Main Result

In this appendix, we provide the detailed proofs for the main results. Before proceeding further, we need to first provide some necessary lemmas. For convenience, denote μ 0 , ϕ 0 , and ψ 0 : = ( ψ 1 , 0 , ψ 2 , 0 , , ψ p , 0 ) as the true values of μ , ϕ , and ψ : = ( ψ 1 , ψ 2 , , ψ p ) , respectively. Write θ 0 : = ( μ 0 , ϕ 0 , ψ 0 ) , and let F t be the sigma field generated by { e s : 1 s t , m + 1 s m + t } .
For convenience, write S t : = ( S t , 1 , S ˜ t ) = ( S t , 1 , S t , 2 , , S t , p + 1 ) , where
S t , 1 = 1 m i = 1 t e i , S t , 2 = 1 m i = 1 t e i ε i 1 , , S t , p + 1 = 1 m i = 1 t e i ε i p ,
for t = 1 , 2 , , n . By following [2], it is easy to check for any s ( 0 , 2 ] that
S [ n s ] W ( s ) : = ( W e ( s ) , W ˜ ( s ) ) = ( W e ( s ) , W ˜ 1 ( s ) , , W ˜ p ( s ) ) ,
under Case (ii) as n , where ‘⇒’ denotes the convergence in space D ( 0 , 2 ] which is the space of real-valued functions of the interval ( 0 , 2 ] that are right continuous and have finite left limits, [ · ] denotes the floor function, and W ( s ) is a vector of Gaussian processes with covariance matrix diag { σ e 2 , Σ 22 } with σ e 2 = E ( e 1 2 ) and
Σ 22 = σ e 2 E ( ε 1 2 ) σ e 2 E ( ε 1 ε p ) σ e 2 E ( ε p ε 1 ) σ e 2 E ( ε p 2 ) .
Lemma A1.
Under the same conditions of Theorem 1, as n , we obtain
  • For Case (i),
    1 m t = p + 1 m Z t ( θ 0 ) d N ( 0 , Σ ) ,
    where ‘ p ’ denotes the convergence in probability, and Σ = diag { Σ 11 , Σ 22 } and
    Σ 11 = σ e 2 0 0 σ e 2 · lim t E X t + m 1 1 + X t + m p 1 2 + i = 1 p ψ i X t + m i 1 1 + X t + m p i 1 2 .
  • For Case (ii),
    1 m t = p + 1 m Z t ( θ 0 ) = ( W e ( 1 ) , ω 2 m * W e ( 2 ) ω m * W e ( 1 ) , W ˜ 1 ( 1 ) , , W ˜ p ( 1 ) ) + o p ( 1 ) ,
    where, for k = m , 2 m ,
    ω k * = X k p 1 1 + X k p 1 2 + i = 1 p ψ i X k p i 1 1 + X k p i 1 .
  • For Case (iii),
    1 m t = p + 1 m Z t ( θ 0 ) d N ( 0 , Σ ˜ ) ,
    where Σ ˜ = diag { σ e 2 , σ e 2 , Σ 22 } .
    for Case (i), or Σ 11 = diag { σ e 2 , σ e 2 } for Cases (ii) and (iii). Note that σ e 2 = E ( e t 2 ) .
Proof of Lemma A1.
Put e t = e t ( θ 0 ) and ε t = ( X t i μ 0 ϕ 0 X t i 1 ) for i = 1 , 2 , , p when θ = θ 0 . Then, it is easy to check that { Z t ( θ 0 ) } is a martingale difference sequence (MDS) with respect to the filter { F t } .
Next, for Case (i), under the conditions of Theorem 1, both { X t } and { ε t } are strictly stationary. Hence, for any a : = ( a 1 , a 2 , , a p + 2 ) R p + 2 , we obtain
1 m t = p + 1 m E ( a Z t ( θ 0 ) ) 2 | F t 1 = a 1 m t = p + 1 m E Z t ( θ 0 ) Z t ( θ 0 ) | F t 1 a p a Σ a ,
as n by using the law of large number for MDS [26]. Similarly, for any arbitrarily small ϵ > 0 , we obtain
1 m t = p + 1 m E | a Z t ( θ 0 ) | 2 I ( | a Z t ( θ 0 ) | ϵ m ) | F t 1 1 ϵ δ 1 n δ 1 / 2 1 m t = p + 1 m E | a Z t ( θ 0 ) | 2 + δ 1 | F t 1 ( p + 2 ) 1 + δ 1 a 2 + δ 1 ϵ δ 1 n δ 1 / 2 1 m t = p + 1 m E i = 1 p + 2 | Z t , i ( θ 0 ) | 2 + δ 1 | F t 1 p 0 ,
by noting that E ( | e t | 2 + δ 1 ) < , which implies E ( | ε t | 2 + δ 1 ) < and in turn E ( | X t | 2 + δ 1 ) < when { ε t } and { X t } are strictly stationary. (A5) and (A6) together show the normality for Case (i) by using the central limit theorem for MDS [26].
For Case (ii), note that
X t + m i 1 = k = 1 p ϕ k 1 ε t + m i k + 1 + ϕ p X t + m i p 1 , i = 1 , 2 , , p ,
and by [2,27], it holds that for any s ( 0 , 1 ]
X [ n s ] n J c ( s ) : = 0 s e c ( s r ) d W ( r ) in the space D ( ( 0 , 1 ] ) ,
as n . Using these, it is easy to check that, as n ,
1 m t = p + 1 m e t + m X t + m 1 1 + X t + m p 1 2 + i = 1 p ψ i X t + m i 1 1 + X t + m p i 1 = 1 m t = p + 1 m e t + m X t + m p 1 1 + X t + m p 1 2 + i = 1 p ψ i X t + m p i 1 1 + X t + m p i 1 + o p ( 1 ) .
Note that
1 m t = p + 1 m e t + m X t + m p 1 1 + X t + m p 1 2 + i = 1 p ψ i X t + m p i 1 1 + X t + m p i 1 : = 1 m t = p + 1 m e t + m ω t + m * = t = p + 1 m ( S t + m , 1 S t + m 1 , 1 ) ω t + m * = S 2 m , 1 ω 2 m * S p + m , 1 ω p + m + 1 * + t = p + 1 m 1 S t + m , 1 ( ω t + m * ω t + m + 1 * ) .
Using (A7), we obtain w 2 m p sgn ( J c ( 2 ) ) , where sgn ( · ) denotes the sign function. As X [ n s ] = O p ( m ) and n for Case (ii), it is easy to check that there exists some d ( 0 , 1 2 ) , for i = 0 , 1 , , p , such that
t = p + 1 m 1 S t + m , 1 X t + m p i 1 1 + X t + m p i 1 X t + m p i 1 + X t + m p i = t = p + 1 m 1 S t + m , 1 X t + m p i X t + m p i 1 ( 1 + ξ t , i , * 2 ) 3 2 O p ( m d ) × 1 m t = p + 1 m 1 { | S t + m , 1 | | X t + m p i X t + m p i 1 | } = o p ( 1 ) ,
where ξ t , i , * lies between X t + m p i and X t + m p i 1 . This shows
1 m t = p + 1 m Z t ( θ 0 ) = ( S m , 1 , S 2 m , 1 ω 2 m * S p + m , 1 ω p + m + 1 * , S m , 2 , , S m , p + 1 ) + o p ( 1 ) , as n .
Then, the asymptotic result for Case (ii) follows immediately based on (A1).
Case (iii) can be proved similarly as Cases (i) and (ii). We omit the details. □
Lemma A2.
Under the same conditions of Theorem 1, as n , we find that
  • (a) 1 m t = p + 1 m Z t ( θ * ) = 1 m t = p + 1 m Z t ( θ 0 ) + O p ( 1 ) , uniformly for ( μ , ϕ ) B ,
  • (b) 1 m t = p + 1 m Z t ( θ * ) Z t ( θ * ) = 1 m t = p + 1 m Z t ( θ 0 ) Z t ( θ 0 ) + o p ( 1 ) = Σ 0 + o p ( 1 ) , uniformly for ( μ , ϕ ) B ,
  • (c) max p + 1 t m sup B Z t ( θ * ) = o p ( m ) ,
    where θ * = ( μ , ϕ , ψ 0 ) ,
    B = { ( μ , ϕ ) : | μ μ 0 | + | ϕ ϕ 0 | < C / m } for Case ( i ) , { ( μ , ϕ ) : | μ μ 0 | + m | ϕ ϕ 0 | < C / m } for Case ( ii ) , { ( μ , ϕ ) : | μ μ 0 | + m | ϕ ϕ 0 | < C / m } for Case ( iii ) ,
    for some positive constant C, Σ 0 = d i a g { Σ 11 , Σ 22 } for Case (i), and d i a g { σ e 2 , σ e 2 , Σ 22 } for Cases (ii) and (iii).
Proof of Lemma A2.
We only prove Parts (a) and (c), as the proof of (b) is trivial based on those of (a) and (c).
For Part (a), note that
e t ( θ * ) e t = X t μ ϕ X t 1 + i = 1 p ψ i , 0 ( X t i μ ϕ X t i 1 ) { X t μ 0 ϕ 0 X t 1 + i = 1 p ψ i , 0 ( X t i μ 0 ϕ 0 X t i 1 ) } = ( μ μ 0 ) ( 1 + i = 1 p ψ i , 0 ) ( ϕ ϕ 0 ) { X t 1 + i = 1 p ψ i , 0 X t i 1 } ,
and ( X t i μ ϕ X t i 1 ) ε t i = ( μ μ 0 ) + ( ϕ ϕ 0 ) X t i 1 for any t = p + 1 , , 2 m and i = 0 , 1 , , p . Hence, we have
sup B 1 m t = p + 1 m ( Z t , 1 ( θ * ) Z t , 1 ( θ 0 ) ) sup B m ( μ μ 0 ) ( 1 + i = 1 p ψ i , 0 ) sup B m ( ϕ ϕ 0 ) 1 m m t = p + 1 m { X t 1 + i = 1 p ψ i , 0 X t i 1 } = C · ( 1 + i = 1 p ψ i , 0 ) · 1 + 0 1 J c ( s ) d s = O p ( 1 ) , as n .
The proofs of sup B 1 m t = p + 1 m ( Z t , k ( θ * ) Z t , k ( θ 0 ) ) , for k = 2 , , p + 2 , follow a similar fashion. This shows Part (a).
For Part (c), based on the decomposition of e t ( θ * ) given in Part (a), we similarly have
max p + 1 t m sup B | Z t , 1 ( θ * ) | max p + 1 t m | e t | + | 1 + i = 1 p ψ i , 0 | · sup B | μ μ 0 | + sup B | ϕ ϕ 0 | · max p + 1 t m | X t 1 + i = 1 p ψ i , 0 X t i 1 | = o p ( m ) ,
by using the Markov inequality based on the conditions of Theorem 1 as n . max p + 1 t m sup B | Z t , k ( θ * ) | = o p ( m ) , k = 2 , , p + 2 , can be proved similarly. We omit the details. This shows Part (c). □
Proof of Theorem 1.
In the following, we only prove Case (ii), as Cases (i) and (iii) follow a similar fashion.
Based on Lemmas A1 and A2, we can show by using similar techniques as in Theorem 1 of [28] that
( μ , ϕ , ψ 0 ) = 1 m t = p + 1 m Z t ( θ * ) Σ 0 1 1 m t = p + 1 m Z t ( θ * ) + o p ( 1 ) ,
uniformly for ( μ , ϕ ) B . Note that ( μ 0 , ϕ 0 ) B . Trivially, with (A9), it follows
( μ 0 , ϕ 0 , ψ 0 ) = 1 m t = p + 1 m Z t ( θ 0 ) Σ 0 1 1 m t = p + 1 m Z t ( θ 0 ) + o p ( 1 ) ,
as n .
ε t ( μ , ϕ ) = X t μ ϕ X t 1 , X ˜ t 1 = X t 1 + i = 1 p ψ i , 0 X t i 1 and γ 0 = 1 + i = 1 p ψ i , 0 . Note that, with (A7), we have
1 m t = p + 1 m { 1 m X ˜ t 1 } d γ 0 0 1 J c ( s ) d s , 1 m t = p + 1 m { 1 m X ˜ t + m 1 ω t + m * } d γ 0 1 2 sgn ( J c ( s ) ) d s , as n .
Next, since ε t = Φ 1 ( B ) e t , where Φ ( B ) = 1 + i = 1 p ψ i , 0 B i with B being the lag operator satisfying B j e t = e t j , which is a linear process of { e t } , we may show that
1 m t = p + 1 m { 1 m X t 1 e t } = O p ( 1 ) , and 1 m t = p + 1 m { 1 m X t 1 ε t } = O p ( 1 ) , as n ,
based on the martingale decomposition as those maintained in [29,30]; see the proof of Theorem 3.1 of [31] for similar discussions. Then, we have
1 m t = p + 1 m { Z t ( θ * ) Z t ( θ 0 ) } = 1 m t = p + 1 m ( μ μ 0 ) γ 0 + ( ϕ ϕ 0 ) X ˜ t 1 { ( μ μ 0 ) γ 0 + ( ϕ ϕ 0 ) X ˜ t + m 1 } ω t + m * { ( μ μ 0 ) γ 0 + ( ϕ ϕ 0 ) X ˜ t 1 } ε t 1 ( μ , ϕ ) { ( μ μ 0 ) γ 0 + ( ϕ ϕ 0 ) X ˜ t 1 } ε t p ( μ , ϕ ) 1 m t = p + 1 m 0 0 e t ( ( μ μ 0 ) + ( ϕ ϕ 0 ) X t 2 ) e t ( ( μ μ 0 ) + ( ϕ ϕ 0 ) X t p 1 ) = γ 0 γ 0 0 1 J c ( s ) d s γ 0 γ 0 1 2 sgn ( J c ( s ) ) d s 0 0 0 0 m ( μ μ 0 ) m ( ϕ ϕ 0 ) + o p ( 1 ) : = Γ m ( μ μ 0 ) m ( ϕ ϕ 0 ) + o p ( 1 ) ,
uniformly for ( μ , ϕ ) B as n .
Based on (A10), it is then easy to check that the minimizer, say ( μ ^ , ϕ ^ ) , of ( μ , ϕ , ψ 0 ) ( μ 0 , ϕ 0 , ψ 0 ) must be in B , and satisfies
m ( μ ^ μ 0 ) m ( ϕ ^ ϕ 0 ) = ( Γ Σ 0 1 Γ ) 1 Γ Σ 0 1 · 1 m t = p + 1 m Z t ( θ 0 ) + o p ( 1 ) , as n .
Then, it follows
( μ ^ , ϕ ^ , ψ 0 ) = 1 m t = p + 1 m Z t ( θ 0 ) Σ 0 1 Σ 0 1 Γ ( Γ Σ 0 1 Γ ) 1 Γ Σ 0 1 1 m t = p + 1 m Z t ( θ 0 ) + o p ( 1 ) .
Further note that
Γ ( Γ Σ 0 1 Γ ) 1 Γ = Γ ( Γ 1 Σ 11 1 Γ 1 ) 1 Γ = Σ 11 1 0 ,
where Γ 1 = γ 0 γ 0 0 1 J c ( s ) d s γ 0 γ 0 1 2 sgn ( J c ( s ) ) d s . Then, we have
Σ 0 1 Σ 0 1 Γ ( Γ Σ 0 1 Γ ) 1 Γ Σ 0 1 = Σ 0 1 Σ 0 1 Σ 11 1 0 Σ 0 1 = 0 Σ 22 1 .
Hence,
p ( ψ 0 ) = ( μ ^ , ϕ ^ , ψ 0 ) = 1 m t = p + 1 m Z t ( θ 0 ) 0 Σ 22 1 1 m t = p + 1 m Z t ( θ 0 ) + o p ( 1 ) = W ˜ ( 1 ) Σ 22 1 W ˜ ( 1 ) + o p ( 1 ) d χ p 2 , as n .
This completes the proof of this Theorem. □

References

  1. Chan, N.H.; Wei, C.Z. Asymptotic inference for nearly nonstationary AR (1) processes. Ann. Stat. 1987, 15, 1050–1063. [Google Scholar] [CrossRef]
  2. Phillips, P.C. Towards a unified asymptotic theory for autoregression. Biometrika 1987, 74, 535–547. [Google Scholar] [CrossRef]
  3. Elliott, G.; Rothenberg, T.J.; Stock, J.H. Efficient Tests for an Autoregressive Unit Root; National Bureau of Economic Research: Cambridge, MA, USA, 1992. [Google Scholar]
  4. Mikusheva, A. Uniform inference in autoregressive models. Econometrica 2007, 75, 1411–1452. [Google Scholar] [CrossRef]
  5. Zhang, R.; Li, C.; Peng, L. Inference for the tail index of a GARCH(1,1) model and an AR(1) model with ARCH(1) errors. Econom. Rev. 2019, 38, 151–169. [Google Scholar] [CrossRef]
  6. Liu, X.; Yang, B.; Cai, Z.; Peng, L. A unified test for predictability of asset returns regardless of properties of predicting variables. J. Econom. 2019, 208, 141–159. [Google Scholar] [CrossRef]
  7. Murray, C.J.; Papell, D.H. The purchasing power parity persistence paradigm. J. Int. Econ. 2002, 56, 1–19. [Google Scholar] [CrossRef] [Green Version]
  8. Rapach, D.E.; Wohar, M.E. The persistence in international real interest rates. Int. J. Financ. Econ. 2004, 9, 339–346. [Google Scholar] [CrossRef]
  9. Datteo, A.; Luca, F.; Busca, G. Statistical pattern recognition approach for long-time monitoring of the G.Meazza stadium by means of AR models and PCA. Eng. Struct. 2017, 153, 317–333. [Google Scholar] [CrossRef]
  10. Said, S.E.; Dickey, D.A. Testing for unit roots in autoregressive-moving average models of unknown order. Biometrika 1984, 71, 599–607. [Google Scholar] [CrossRef]
  11. Phillips, P.C.; Perron, P. Testing for a unit root in time series regression. Biometrika 1988, 75, 335–346. [Google Scholar] [CrossRef]
  12. Kwiatkowski, D.; Phillips, P.C.; Schmidt, P.; Shin, Y. Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root? J. Econom. 1992, 54, 159–178. [Google Scholar] [CrossRef]
  13. Hill, J.; Li, D.; Peng, L. Uniform interval estimation for an AR (1) process with AR errors. Stat. Sin. 2016, 26, 119–136. [Google Scholar] [CrossRef]
  14. Li, C.; Li, D.; Peng, L. Uniform test for predictive regression with AR errors. J. Bus. Econ. Stat. 2017, 35, 29–39. [Google Scholar] [CrossRef]
  15. Silvey, S.D. The Lagrangian multiplier test. Ann. Math. Stat. 1959, 30, 389–407. [Google Scholar] [CrossRef]
  16. Box, G.E.; Pierce, D.A. Distribution of residual autocorrelations in autoregressive-integrated moving average time series models. J. Am. Stat. Assoc. 1970, 65, 1509–1526. [Google Scholar] [CrossRef]
  17. Ljung, G.M.; Box, G.E. On a measure of lack of fit in time series models. Biometrika 1978, 65, 297–303. [Google Scholar] [CrossRef]
  18. Owen, A.B. Empirical likelihood ratio confidence intervals for a single functional. Biometrika 1988, 75, 237–249. [Google Scholar] [CrossRef]
  19. Owen, A. Empirical likelihood ratio confidence regions. Ann. Stat. 1990, 18, 90–120. [Google Scholar] [CrossRef]
  20. Owen, A. Empirical likelihood for linear models. Ann. Stat. 1991, 19, 1725–1747. [Google Scholar] [CrossRef]
  21. Qin, J.; Lawless, J. Empirical likelihood and general estimating equations. Ann. Stat. 1994, 22, 300–325. [Google Scholar] [CrossRef]
  22. Chan, N.H.; Li, D.; Peng, L. Toward a unified interval estimation of autoregressions. Econom. Theory 2012, 28, 705–717. [Google Scholar] [CrossRef]
  23. Liu, F.; Wang, P.F.; Chun-Yan, L.I.; Kang, X.M.; Statistics, D.O. Empirical Likelihood Based Diagnostics for Heteroscedasticity in Single-index Model. J. Chongqing Technol. Bus. Univ. (Nat. Sci. Ed.) 2018, 81, 255–281. [Google Scholar]
  24. Zhu, F.; Cai, Z.; Peng, L. Predictive regressions for macroeconomic data. Ann. Appl. Stat. 2014, 8, 577–594. [Google Scholar] [CrossRef] [Green Version]
  25. Liu, X.; Liu, Y.; Lu, F. Empirical likelihood-based unified confidence region for a predictive regression model. Commun. Stat. Simul. Comput. 2019, 51, 2122–2139. [Google Scholar] [CrossRef]
  26. Hall, P.; Heyde, C.C. Martingale Limit Theory and Its Application; Academic Press: Cambridge, MA, USA, 1981. [Google Scholar]
  27. Phillips, P.C. Time series regression with a unit root. Econom. J. Econom. Soc. 1987, 55, 277–301. [Google Scholar] [CrossRef] [Green Version]
  28. Ma, Y.; Zhou, M.; Peng, L.; Zhang, R. Test for zero median of errors in an arma–garch model. Econom. Theory 2022, 38, 536–561. [Google Scholar] [CrossRef]
  29. Phillips, P.C.; Solo, V. Asymptotics for linear processes. Ann. Stat. 1992, 20, 971–1001. [Google Scholar] [CrossRef]
  30. Phillips, P.C.; Magdalinos, T. Limit theory for moderate deviations from a unit root. J. Econom. 2007, 136, 115–130. [Google Scholar] [CrossRef] [Green Version]
  31. Guo, G.; Sun, Y.; Wang, S. Testing for moderate explosiveness. Econom. J. 2019, 22, 73–95. [Google Scholar] [CrossRef]
Figure 1. The power performance of EL method at τ = 0.05 .
Figure 1. The power performance of EL method at τ = 0.05 .
Axioms 11 00690 g001
Figure 2. Empirical reject probabilities at τ = 0.05 .
Figure 2. Empirical reject probabilities at τ = 0.05 .
Axioms 11 00690 g002
Figure 3. The results of the test for AR error structure between EL and AIC.
Figure 3. The results of the test for AR error structure between EL and AIC.
Axioms 11 00690 g003
Figure 4. Time series graphs of 8 countries.
Figure 4. Time series graphs of 8 countries.
Axioms 11 00690 g004
Table 1. Empirical reject probabilities at τ = 0.05 .
Table 1. Empirical reject probabilities at τ = 0.05 .
ELAIC
Scenarios μ ϕ 300 600 900 1200 300 600 900 1200
S100.50.06010.05180.05410.05540.07150.06930.07620.0768
1 1 n 0.07210.06590.06090.06000.22850.23080.23370.2404
10.07330.06570.05960.05970.23280.23370.23570.2396
0.010.50.06990.05750.05090.05090.07600.07190.07320.0698
1 1 n 0.07100.06110.06080.05880.23010.23230.24070.2301
10.07330.05780.06120.05810.22670.23300.23000.2406
S200.50.05870.04640.04690.05171.00001.00001.00001.0000
1 1 n 0.06520.06200.06290.05700.55310.34580.26350.2243
10.06890.05870.06300.05630.55050.34500.25740.2246
0.010.50.05910.04540.04810.05140.99991.00001.00001.0000
1 1 n 0.06910.05460.05860.06070.54810.33940.26650.2339
10.06940.05890.05870.05780.54540.34630.26610.2254
S300.50.05570.04980.04740.04630.89690.84660.81310.7832
1 1 n 0.07780.06130.06210.05900.59630.32840.22380.1732
10.08620.06350.06040.06010.59780.33180.22520.1789
0.010.50.05500.04960.04830.04610.90050.84450.81950.7848
1 1 n 0.08590.06920.06010.05870.60410.33320.21730.1832
10.08340.06460.05780.05850.59710.33180.22010.1812
Table 2. The power performance of unit root test at τ = 0.05 .
Table 2. The power performance of unit root test at τ = 0.05 .
Right OrderWrong Order
ϕ 30060090012003006009001200
1 1 n 0.04740.04550.04220.04370.04380.04530.04190.0425
1 3 n 0.10290.09410.09230.08620.09980.09180.08970.0851
1 5 n 0.17950.17300.17180.16780.17470.17160.16660.1664
1 10 n 0.47480.45960.44550.43670.44910.44960.43540.4307
1 15 n 0.76200.73430.71590.69830.72760.71680.70280.6864
Table 3. Test results of 8 countries.
Table 3. Test results of 8 countries.
Coutry μ ^ ϕ ^ ELp-ValuesLB-TestAIC
Canada0.03400.9992 A R ( 0 ) 0.2122 7.7250 × 10 3 *** A R ( 3 )
China0.01690.9975 A R ( 0 ) 0.85720.1391 A R ( 0 )
Norway0.04620.9947 A R ( 0 ) 0.4315 1.9450 × 10 6 *** A R ( 2 )
Singapore0.01460.9893 A R ( 0 ) 0.8393 0.0238 ** A R ( 3 )
Thailand0.14530.9953 A R ( 0 ) 0.1666 4.0220 × 10 3 *** A R ( 1 )
Sweden0.02740.9968 A R ( 0 ) 0.2165 5.6510 × 10 3 *** A R ( 3 )
Japan1.93720.9822 A R ( 0 ) 0.23920.8896 A R ( 0 )
Sri Lanka0.17380.9992 A R ( 2 ) 0.2622 7.2230 × 10 6 *** A R ( 2 )
Significance levels: * p 0.1 , ** p 0.05 , *** p 0.01 .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, X.; Liu, X.; Fan, Y.; Tan, L.; Liu, Q. A Unified Test for the AR Error Structure of an Autoregressive Model. Axioms 2022, 11, 690. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11120690

AMA Style

Wei X, Liu X, Fan Y, Tan L, Liu Q. A Unified Test for the AR Error Structure of an Autoregressive Model. Axioms. 2022; 11(12):690. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11120690

Chicago/Turabian Style

Wei, Xinyi, Xiaohui Liu, Yawen Fan, Li Tan, and Qing Liu. 2022. "A Unified Test for the AR Error Structure of an Autoregressive Model" Axioms 11, no. 12: 690. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms11120690

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop