Modeling and forecasting of stock index volatility with APARCH models under ordered restriction

This article examines volatility models for modeling and forecasting the Standard & Poor 500 (S&P 500) daily stock index returns, including the autoregressive moving average, the Taylor and Schwert generalized autoregressive conditional heteroscedasticity (GARCH), the Glosten, Jagannathan and Runkle GARCH and asymmetric power ARCH (APARCH) with the following conditional distributions: normal, Student's t and skewed Student's t‐distributions. In addition, we undertake unit root (augmented Dickey–Fuller and Phillip–Perron) tests, co‐integration test and error correction model. We study the stationary APARCH (p) model with parameters, and the uniform convergence, strong consistency and asymptotic normality are prove under simple ordered restriction. In fitting these models to S&P 500 daily stock index return data over the period 1 January 2002 to 31 December 2012, we found that the APARCH model using a skewed Student's t‐distribution is the most effective and successful for modeling and forecasting the daily stock index returns series. The results of this study would be of great value to policy makers and investors in managing risk in stock markets trading.     

Spurious Fixed Effects Regression*

This article shows that spurious regression results can occur for a fixed effects model with weak time series variation in the regressor and/or strong time series variation in the regression errors when the first‐differenced and Within‐OLS estimators are used. Asymptotic properties of these estimators and the related t‐tests and model selection criteria are studied by sending the number of cross‐sectional observations to infinity. This article shows that the first‐differenced and Within‐OLS estimators diverge in probability, that the related t‐tests are inconsistent, that R²s converge to zero in probability and that AIC and BIC diverge to ?∞ in probability. The results of the article warn that one should not jump to the use of fixed effects regressions without considering the degree of time series variations in the data.     

Hypothesis testing based on goodness-of-fit in the moving average time series model

Small sample properties of t-tests are compared with those of tests based on relative goodness- of-fit in the context of the first order moving average time series model. Monte Carlo experiments reported in the paper suggest that the actual size of these t-tests greatly exceeds theoretical large sample significance levels, while conformity of goodness-of-fit statistics to the appropriate chi-square or F-distributions is much closer. The evidence presented suggests that practitioners are well advised to employ goodness-of-fit tests as a check on results of t-tests particularly when the latter indicate ‘significance’.     

Further Evidence on the Statistical Properties of Real GNP*

The well‐known lack of power of unit‐root tests has often been attributed to the short length of macroeconomic variables and also to data‐generating processes (DGPs) departing from the I(1)–I(0) models. This paper shows that by using long spans of annual real gross national product (GNP) and GNP per capita (133 years), high power can be achieved, leading to the rejection of both the unit‐root and the trend‐stationary hypothesis. More flexible representations are then considered, namely, processes containing structural breaks (SB) and fractional orders of integration (FI). Economic justification for the presence of these features in GNP is provided. It is shown that both FI and SB formulations are in general preferred to the autoregressive integrated moving average (ARIMA) [I(1) or I(0)] formulations. As a novelty in this literature, new techniques are applied to discriminate between FI and SB. It turns out that the FI specification is preferred, implying that GNP and GNP per capita are non‐stationary, highly persistent but mean‐reverting series. Finally, it is shown that the results are robust when breaks in the deterministic component are allowed for in the FI model. Some macroeconomic implications of these findings are also discussed.     

Exploiting tail shape biases to discriminate between stable and student t alternatives

The nonnormal stable laws and Student t distributions are used to model the unconditional distribution of financial asset returns, as both models display heavy tails. The relevance of the two models is subject to debate because empirical estimates of the tail shape conditional on either model give conflicting signals. This stems from opposing bias terms. We exploit the biases to discriminate between the two distributions. A sign estimator for the second‐order scale parameter strengthens our results. Tail estimates based on asset return data match the bias induced by finite‐variance unconditional Student t data and the generalized autoregressive conditional heteroscedasticity process.     

Prewhitening Bias in HAC Estimation*

Heteroskedasticity and autocorrelation consistent (HAC) estimation commonly involves the use of prewhitening filters based on simple autoregressive models. In such applications, small sample bias in the estimation of autoregressive coefficients is transmitted to the recolouring filter, leading to HAC variance estimates that can be badly biased. The present paper provides an analysis of these issues using asymptotic expansions and simulations. The approach we recommend involves the use of recursive demeaning procedures that mitigate the effects of small‐sample autoregressive bias. Moreover, a commonly used restriction rule on the prewhitening estimates (that first‐order autoregressive coefficient estimates, or largest eigenvalues, >0.97 be replaced by 0.97) adversely interferes with the power of unit‐root and [ Kwiatkowski, Phillips, Schmidt and Shin (1992) Journal of Econometrics, Vol. 54, pp. 159–178] (KPSS) tests. We provide a new boundary condition rule that improves the size and power properties of these tests. Some illustrations of the effects of these adjustments on the size and power of KPSS testing are given. Using prewhitened HAC estimates and the new boundary condition rule, the KPSS test is consistent, in contrast to KPSS testing that uses conventional prewhitened HAC estimates [ Lee, J. S. (1996) Economics Letters, Vol. 51, pp. 131–137].     

Distribution of residual autocorrelations for multiplicative seasonal ARMA models with uncorrelated but nonindependent error terms

In this paper, we consider portmanteau tests for testing the adequacy of multiplicative seasonal autoregressive moving‐average models under the assumption that the errors are uncorrelated but not necessarily independent. We relax the standard independence assumption on the error terms in order to extend the range of applications of the seasonal autoregressive moving‐average models. We study the asymptotic distributions of residual and normalized residual empirical autocovariances and autocorrelations under weak assumptions on noise. We establish the asymptotic behavior of the proposed statistics. A set of Monte Carlo experiments and an application to monthly mean total sunspot number are presented.     

A model of fractional cointegration, and tests for cointegration using the bootstrap

The paper proposes a framework for modelling cointegration in fractionally integrated processes, and considers methods for testing the existence of cointegrating relationships using the parametric bootstrap. In these procedures, ARFIMA models are fitted to the data, and the estimates used to simulate the null hypothesis of non-cointegration in a vector autoregressive modelling framework. The simulations are used to estimate p-values for alternative regression-based test statistics, including the F goodness-of-fit statistic, the Durbin–Watson statistic and estimates of the residual d. The bootstrap distributions are economical to compute, being conditioned on the actual sample values of all but the dependent variable in the regression. The procedures are easily adapted to test stronger null hypotheses, such as statistical independence. The tests are not in general asymptotically pivotal, but implemented by the bootstrap, are shown to be consistent against alternatives with both stationary and nonstationary cointegrating residuals. As an example, the tests are applied to the series for UK consumption and disposable income. The power properties of the tests are studied by simulations of artificial cointegrating relationships based on the sample data. The F test performs better in these experiments than the residual-based tests, although the Durbin–Watson in turn dominates the test based on the residual d.     

On Tests for a Normal Mean with Known Coefficient of Variation

Hinkley (1977) derived two tests for testing the mean of a normal distribution with known coefficient of variation (c.v.) for right alternatives. They are the locally most powerful (LMP) and the conditional tests based on the ancillary statistic for μ. In this paper, the likelihood ratio (LR) and Wald tests are derived for the one‐ and two‐sided alternatives, as well as the two‐sided version of the LMP test. The performances of these tests are compared with those of the classical t, sign and Wilcoxon signed rank tests. The latter three tests do not use the information on c.v. Normal approximation is used to approximate the null distribution of the test statistics except for the t test. Simulation results indicate that all the tests maintain the type‐I error rates, that is, the attained level is close to the nominal level of significance of the tests. The power functions of the tests are estimated through simulation. The power comparison indicates that for one‐sided alternatives the LMP test is the best test whereas for the two‐sided alternatives the LR or the Wald test is the best test. The t, sign and Wilcoxon signed rank tests have lower power than the LMP, LR and Wald tests at various alternative values of μ. The power difference is quite large in several simulation configurations. Further, it is observed that the t, sign and Wilcoxon signed rank tests have considerably lower power even for the alternatives which are far away from the null hypothesis when the c.v. is large. To study the sensitivity of the tests for the violation of the normality assumption, the type I error rates are estimated on the observations of lognormal, gamma and uniform distributions. The newly derived tests maintain the type I error rates for moderate values of c.v.     

Inference with dependent data using cluster covariance estimators

This paper presents an inference approach for dependent data in time series, spatial, and panel data applications. The method involves constructing t and Wald statistics using a cluster covariance matrix estimator (CCE). We use an approximation that takes the number of clusters/groups as fixed and the number of observations per group to be large. The resulting limiting distributions of the t and Wald statistics are standard t and F distributions where the number of groups plays the role of sample size. Using a small number of groups is analogous to ‘fixed-b’ asymptotics of [Kiefer and Vogelsang, 2002] and [Kiefer and Vogelsang, 2005] (KV) for heteroskedasticity and autocorrelation consistent inference. We provide simulation evidence that demonstrates that the procedure substantially outperforms conventional inference procedures.     

A proof of a relationship between the generalized variances for associated autoregressive and moving average processes

Summary In this paper we give a simple proof of the result that, for any integer,r, given two processes of orderr, one autoregressive and the other moving average but both with the same parameters, then the generalized variance of all ordersk2r, for the autoregressive process, is exactly equal to the infinite order generalized variance for the moving average process.     

Asymptotic distributions of some scale estimators in nonlinear models

Often in the robust analysis of regression and time series models there is a need for having a robust scale estimator of a scale parameter of the errors. One often used scale estimator is the median of the absolute residuals s ₁. It is of interest to know its limiting distribution and the consistency rate. Its limiting distribution generally depends on the estimator of the regression and/or autoregressive parameter vector unless the errors are symmetrically distributed around zero. To overcome this difficulty it is then natural to use the median of the absolute differences of pairwise residuals, s ₂, as a scale estimator. This paper derives the asymptotic distributions of these two estimators for a large class of nonlinear regression and autoregressive models when the errors are independent and identically distributed. It is found that the asymptotic distribution of a suitably standardizes s ₂ is free of the initial estimator of the regression/autoregressive parameters. A similar conclusion also holds for s ₁ in linear regression models through the origin and with centered designs, and in linear autoregressive models with zero mean errors.  This paper also investigates the limiting distributions of these estimators in nonlinear regression models with long memory moving average errors. An interesting finding is that if the errors are symmetric around zero, then not only is the limiting distribution of a suitably standardized s ₁ free of the regression estimator, but it is degenerate at zero. On the other hand a similarly standardized s ₂ converges in distribution to a normal distribution, regardless of the errors being symmetric or not. One clear conclusion is that under the symmetry of the long memory moving average errors, the rate of consistency for s ₁ is faster than that of s ₂.     

A Comparison of Johansen's,Bierens’ and the Subspace Algorithm Method for Cointegration Analysis*

The methods listed in the title are compared by means of a simulation study and a real world application. The aspects compared via simulations are the performance of the tests for the cointegrating rank and the quality of the estimated cointegrating space. The subspace algorithm method, formulated in the state space framework and thus applicable for vector autoregressive moving average (VARMA) processes, performs at least comparably to the Johansen method. Both the Johansen procedure and the subspace algorithm cointegration analysis perform significantly better than Bierens’ method. The real‐world application is an investigation of the long‐run properties of the one‐sector neoclassical growth model for Austria. The results do not fully support the implications of the model with respect to cointegration. Furthermore, the results differ greatly between the different methods. The amount of variability depends strongly upon the number of variables considered and huge differences occur for the full system with six variables. Therefore we conclude that the results of such applications with about five or six variables and 100 observations, which are typical in the applied literature, should possibly be interpreted with more caution than is commonly done.     

Diagnostic analysis for a vector autoregressive model under Student′s t‐distributions

In this paper, we use the local influence method to study a vector autoregressive model under Student^′s t‐distributions. We present the maximum likelihood estimators and the information matrix. We establish the normal curvature diagnostics for the vector autoregressive model under three usual perturbation schemes for identifying possible influential observations. The effectiveness of the proposed diagnostics is examined by a simulation study, followed by our data analysis using the model to fit the weekly log returns of Chevron stock and the Standard & Poor's 500 Index as an application.     

Power of Tests for Unit Roots in the Presence of a Linear Trend*

Dickey and Fuller [Econometrica (1981) Vol. 49, pp. 1057–1072] suggested unit‐root tests for an autoregressive model with a linear trend conditional on an initial observation. TPower of tests for unit roots in the presence of a linear trendightly different model with a random initial value in which nuisance parameters can easily be eliminated by an invariant reduction of the model. We show that invariance arguments can also be used when comparing power within a conditional model. In the context of the conditional model, the Dickey–Fuller test is shown to be more stringent than a number of unit‐root tests motivated by models with random initial value. The power of the Dickey–Fuller test can be improved by making assumptions to the initial value. The practitioner therefore has to trade‐off robustness and power, as assumptions about initial values are hard to test, but can give more power.     

A flexible parametric GARCH model with an application to exchange rates

Many asset prices, including exchange rates, exhibit periods of stability punctuated by infrequent, substantial, often one‐sided adjustments. Statistically, this generates empirical distributions of exchange rate changes that exhibit high peaks, long tails, and skewness. This paper introduces a GARCH model, with a flexible parametric error distribution based on the exponential generalized beta (EGB) family of distributions. Applied to daily US dollar exchange rate data for six major currencies, evidence based on a comparison of actual and predicted higher‐order moments and goodness‐of‐fit tests favours the GARCH‐EGB2 model over more conventional GARCH‐t and EGARCH‐t model alternatives, particularly for exchange rate data characterized by skewness. Copyright © 2001 John Wiley & Sons, Ltd.     

Fluctuation Tests for a Change in Persistence*

In this paper, we develop a set of new persistence change tests which are similar in spirit to those of Kim [Journal of Econometrics (2000) Vol. 95, pp. 97–116], Kim et al. [Journal of Econometrics (2002) Vol. 109, pp. 389–392] and Busetti and Taylor [Journal of Econometrics (2004) Vol. 123, pp. 33–66]. While the exisiting tests are based on ratios of sub‐sample Kwiatkowski et al. [Journal of Econometrics (1992) Vol. 54, pp. 158–179]‐type statistics, our proposed tests are based on the corresponding functions of sub‐sample implementations of the well‐known maximal recursive‐estimates and re‐scaled range fluctuation statistics. Our statistics are used to test the null hypothesis that a time series displays constant trend stationarity [I(0)] behaviour against the alternative of a change in persistence either from trend stationarity to difference stationarity [I(1)], or vice versa. Representations for the limiting null distributions of the new statistics are derived and both finite‐sample and asymptotic critical values are provided. The consistency of the tests against persistence change processes is also demonstrated. Numerical evidence suggests that our proposed tests provide a useful complement to the extant persistence change tests. An application of the tests to US inflation rate data is provided.     

Multivariate t and beta distributions associated with the multivariate F distribution

Relationships between F, skew t and beta distributions in the univariate case are in this paper extended in a natural way to the multivariate case. The result is two new distributions: a multivariate t/skew t distribution (on ℜ^m) and a multivariate beta distribution (on (0,1)^m). A special case of the former distribution is a new multivariate symmetric t distribution. The new distributions have a natural relationship to the standard multivariate F distribution (on (ℜ⁺)^m) and many of their properties run in parallel. We look at: joint distributions, mathematically and graphically; marginal and conditional distributions; moments; correlations; local dependence; and some limiting cases. Received: March 2001     

Modelling dependence using skew t copulas: Bayesian inference and applications

We construct a copula from the skew t distribution of Sahu et al. ( 2003 ). This copula can capture asymmetric and extreme dependence between variables, and is one of the few copulas that can do so and still be used in high dimensions effectively. However, it is difficult to estimate the copula model by maximum likelihood when the multivariate dimension is high, or when some or all of the marginal distributions are discrete‐valued, or when the parameters in the marginal distributions and copula are estimated jointly. We therefore propose a Bayesian approach that overcomes all these problems. The computations are undertaken using a Markov chain Monte Carlo simulation method which exploits the conditionally Gaussian representation of the skew t distribution. We employ the approach in two contemporary econometric studies. The first is the modelling of regional spot prices in the Australian electricity market. Here, we observe complex non‐Gaussian margins and nonlinear inter‐regional dependence. Accurate characterization of this dependence is important for the study of market integration and risk management purposes. The second is the modelling of ordinal exposure measures for 15 major websites. Dependence between websites is important when measuring the impact of multi‐site advertising campaigns. In both cases the skew t copula substantially outperforms symmetric elliptical copula alternatives, demonstrating that the skew t copula is a powerful modelling tool when coupled with Bayesian inference. Copyright © 2010 John Wiley & Sons, Ltd.     

Bootstrap unit root tests in panels with cross-sectional dependency

We apply bootstrap methodology to unit root tests for dependent panels with N cross-sectional units and T time series observations. More specifically, we let each panel be driven by a general linear process which may be different across cross-sectional units, and approximate it by a finite order autoregressive integrated process of order increasing with T. As we allow the dependency among the innovations generating the individual series, we construct our unit root tests from the estimation of the system of the entire N cross-sectional units. The limit distributions of the tests are derived by passing T to infinity, with N fixed. We then apply bootstrap method to the approximated autoregressions to obtain critical values for the panel unit root tests, and establish the asymptotic validity of such bootstrap panel unit root tests under general conditions. The proposed bootstrap tests are indeed quite general covering a wide class of panel models. They in particular allow for very general dynamic structures which may vary across individual units, and more importantly for the presence of arbitrary cross-sectional dependency. The finite sample performance of the bootstrap tests is examined via simulations, and compared to that of commonly used panel unit root tests. We find that our bootstrap tests perform relatively well, especially when N is small.