Response

The aim of the paper is to consider the implicit restrictions imposed when adopting an AR(1) error term in the context of the linear regression model. It is shown that these restrictions amount to assuming a largely identical temporal structure for all the variables involved in the specification. Implicit in this is the assumption that these variables are mutually Granger non-causal. The main implication of this result is that in most cases when residual autocorrelation is detected boththe OLS and GLS estimators are biased and inconsistent.     

Slow-explosive AR(1) processes converging to random walk

Abstract

This article investigates slow-explosive AR(1) processes, which converge to a random walk (RW) process with logarithm rates, to fill the gap between nearly non-stationary AR(1) and moderately deviated AR(1) processes, and derives the asymptotics of the least squares estimator using central limit theorems for (reduced) U-statistic. We successfully establish the smooth link between the nearly non-stationary AR(1) and the moderately deviated AR(1) processes. Some novel results are reported, which include the convergence of the least squares estimator to a biased fractional Brownian motion.     

Goodness-of-fit tests for binomial AR(1) processes

The binomial AR(1) model describes a nonlinear process with a first-order autoregressive (AR(1)) structure and a binomial marginal distribution. To develop goodness-of-fit tests for the binomial AR(1) model, we investigate the observed marginal distribution of the binomial AR(1) process, and we tackle its autocorrelation structure. Motivated by the family of power-divergence statistics for handling discrete multivariate data, we derive the asymptotic distribution of certain categorized power-divergence statistics for the case of a binomial AR(1) process. Then we consider Bartlett's formula, which is widely used in time series analysis to provide estimates of the asymptotic covariance between sample autocorrelations, but which is not applicable when the underlying process is nonlinear. Hence, we derive a novel Bartlett-type formula for the asymptotic distribution of the sample autocorrelations of a binomial AR(1) process, which is then applied to develop tests concerning the autocorrelation structure. Simulation studies are carried out to evaluate the size and power of the proposed tests under diverse alternative process models. Several real examples are used to illustrate our methods and findings.     

Binomial AR(1) processes: moments,cumulants, and estimation

The modelling and analysis of count-data time series are areas of emerging interest with various applications in practice. We consider the particular case of the binomial AR(1) model, which is well suited for describing binomial counts with a first-order autoregressive serial dependence structure. We derive explicit expressions for the joint (central) moments and cumulants up to order 4. Then, we apply these results for expressing moments and asymptotic distribution of the squared difference estimator as an alternative to the sample autocovariance. We also analyse the asymptotic distribution of the conditional least-squares estimators of the parameters of the binomial AR(1) model. The finite-sample performance of these estimators is investigated in a simulation study, and we apply them to real data about computerized workstations.     

Extended binomial AR(1) processes with generalized binomial thinning operator

Abstract

In this article, we introduce an extended binomial AR(1) model based on the generalized binomial thinning operator. This operator relaxes the independence assumption of the binomial thinning operator and contains dependent Bernoulli counting series. The new model contains the binomial AR(1) model as a particular case. Some probabilistic and statistical properties are explored. Estimators of the model parameters are derived by conditional maximum likelihood (CML), conditional least squares (CLS) and weighted conditional least squares (WCLS) methods. Some asymptotic properties and numerical results of the estimators are studied. The good performance of the new model is illustrated, among other competitive models in the literature, by an application to the monthly drunken driving counts.     

中国财政收入与GDP之间关系的协整分析与误差修正模型研究

一、引言财政收入是政府部门的公共收入,是国民收入分配中用于保证政府行使其公共职能、实施公共政策以及提供公共服务的资金需求。国内生产总值是反映一个国家(地区)在一定时期内国民经济活动最终成果的总量指标。从生产的角度看,它是国民经济各部门新创造的增加值的总和;从使     

Inference for the tail index of a GARCH(1,1) model and an AR(1) model with ARCH(1) errors

For a GARCH(1,1) sequence or an AR(1) model with ARCH(1) errors, one can estimate the tail index by solving an estimating equation with unknown parameters replaced by the quasi maximum likelihood estimation, and a profile empirical likelihood method can be employed to effectively construct a confidence interval for the tail index. However, this requires that the errors of such a model have at least a finite fourth moment. In this article, we show that the finite fourth moment can be relaxed by employing a least absolute deviations estimate for the unknown parameters by noting that the estimating equation for determining the tail index is invariant to a scale transformation of the underlying model.     

On Small Sample Properties of the Wald,LR and LM Tests in a Linear Model with AR(1) Errors

When the error terms are autocorrelated, the conventional t-tests for individual regression coefficients mislead us to over-rejection of the null hypothesis. We examine, by Monte Carlo experiments, the small sample properties of the unrestricted estimator of ρ and of the estimator of ρ restricted by the null hypothesis. We compare the small sample properties of the Wald, likelihood ratio and Lagrange multiplier test statistics for individual regression coefficients. It is shown that when the null hypothesis is true, the unrestricted estimator of ρ is biased. It is also shown that the Lagrange multiplier test using the maximum likelihood estimator of ρ performs better than the Wald and likelihood ratio tests.     

On least-squares bias in the AR(p) models: Bias correction using the bootstrap methods

In the case where the lagged dependent variables are included in the regression model, it is known that the ordinary least squares estimates (OLSE) are biased in small sample and that the bias increases as the number of the irrelevant variables increases. In this paper, based on the bootstrap methods, an attempt is made to obtain the unbiased estimates in autoregressive and non-Gaussian cases. We propose the residual-based bootstrap method in this paper. Some simulation studies are performed to examine whether the proposed estimation procedure works well or not. We obtain the results that it is possible to recover the true parameter values and that the proposed procedure gives us the less biased estimators than OLSE. This paper is a substantial revision of Tanizaki (2000). The normality assumption is adopted in Tanizaki (2000), but it is not required in this paper. The authors are grateful to an anonymous referee for valuable suggestions and comments. This research was partially supported by Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research (C)(2) #14530033, 2002–2005, for H. Tanizaki and Grants-in-Aid for the 21st Century COE Program.     

Efficiency and Validity Analyses of Two-Stage Estimation Procedures and Derived Testing Procedures in Quantitative Linear Models with AR(1) Errors

Abstract

In a quantitative linear model with errors following a stationary Gaussian, first-order autoregressive or AR(1) process, Generalized Least Squares (GLS) on raw data and Ordinary Least Squares (OLS) on prewhitened data are efficient methods of estimation of the slope parameters when the autocorrelation parameter of the error AR(1) process, ρ, is known. In practice, ρ is generally unknown. In the so-called two-stage estimation procedures, ρ is then estimated first before using the estimate of ρ to transform the data and estimate the slope parameters by OLS on the transformed data. Different estimators of ρ have been considered in previous studies. In this article, we study nine two-stage estimation procedures for their efficiency in estimating the slope parameters. Six of them (i.e., three noniterative, three iterative) are based on three estimators of ρ that have been considered previously. Two more (i.e., one noniterative, one iterative) are based on a new estimator of ρ that we propose: it is provided by the sample autocorrelation coefficient of the OLS residuals at lag 1, denoted r(1). Lastly, REstricted Maximum Likelihood (REML) represents a different type of two-stage estimation procedure whose efficiency has not been compared to the others yet. We also study the validity of the testing procedures derived from GLS and the nine two-stage estimation procedures. Efficiency and validity are analyzed in a Monte Carlo study. Three types of explanatory variable x in a simple quantitative linear model with AR(1) errors are considered in the time domain: Case 1, x is fixed; Case 2, x is purely random; and Case 3, x follows an AR(1) process with the same autocorrelation parameter value as the error AR(1) process. In a preliminary step, the number of inadmissible estimates and the efficiency of the different estimators of ρ are compared empirically, whereas their approximate expected value in finite samples and their asymptotic variance are derived theoretically. Thereafter, the efficiency of the estimation procedures and the validity of the derived testing procedures are discussed in terms of the sample size and the magnitude and sign of ρ. The noniterative two-stage estimation procedure based on the new estimator of ρ is shown to be more efficient for moderate values of ρ at small sample sizes. With the exception of small sample sizes, REML and its derived F-test perform the best overall. The asymptotic equivalence of two-stage estimation procedures, besides REML, is observed empirically. Differences related to the nature, fixed or random (uncorrelated or autocorrelated), of the explanatory variable are also discussed.     

Modeling Technology as a Dynamic Error Components Process: The Case of the Inter-country Agricultural Production Function

In this paper, we propose a dynamic error-components model to represent the unobserved level of technology. This specification implies a well-defined common factor dynamic model for per capita output that can be tested explicitly. The model is applied to data on aggregates of agricultural inputs and outputs for groups of countries from the OECD, Africa (AF), Latin America (LA) as well as centrally planned countries, over a period of 31 years. We find that the proposed model fits the data better than alternative static specifications and satisfies the implied common factor restrictions in two of the samples. The results suggest that although technological change seems to have been a faster process for less developed countries relative to the OECD countries, it has not been fast enough to reduce appreciably the enormous differences in average technological levels that still persist between them.     

On the statistical inference of a machine-generated autoregressive AR(1) model

We have obtained the asymptotic bias and the limiting distribution for the Yule–Walker estimator of the autoregressive parameter under a considerably weaker assumption than that of independence in the noise sequence. Among other things, these suggest robustness of the classical results and throw some light on the use of simulations based on pseudorandom numbers in verifying these results.     

An Approximation Model of the Collective Risk Model with INAR(1) Claim Process

Cossette et al. (2010 Cossette, H., Marceau, E., Maume-Deschamps, V. (2010). Discerte-time risk models based on time series for count random variables. ASTIN Bull. 40:123–150.[Crossref], [Web of Science ®] , [Google Scholar], 2011 Cossette, H., Marceau, E., Toureille, F. (2011). risk models based on time series for count random variables. Insur. Math. Econ. 48:19–28.[Crossref], [Web of Science ®] , [Google Scholar]) gave a novel collective risk model where the total numbers of claims satisfy the first-order integer-valued autoregressive process. For a risk model, it is interesting to investigate the upper bound of ruin probability. However, the loss increments of the above model are dependent; it is difficult to derive the upper bound of ruin probability. In this article, we propose an approximation model with stationary independent increments. The upper bound of ruin probability and the adjustment coefficient are derived. The approximation model is illustrated via four simulated examples. Results show that the gap of the approximation model and dependent model can be ignored by adjusting values of parameters.     

Canonical correlation analysis for the vector AR(1) model with ARCH innovations

This paper extends the results of canonical correlation analysis of Anderson [2002. Canonical correlation analysis and reduced-rank regression in autoregressive models. Ann. Statist. 30, 1134–1154] to a vector AR(1) process with a vector ARCH(1) innovations. We obtain the limiting distributions of the sample matrices, the canonical correlations and the canonical vectors of the process. The extension is important because many time series in economics and finance exhibit conditional heteroscedasticity. We also use simulation to demonstrate the effects of ARCH innovations on the canonical correlation analysis in finite sample. Both the limiting distributions and simulation results show that overlooking the ARCH effects in canonical correlation analysis can easily lead to erroneous inference.     

On the estimation of missing values in AR(1) model with exponential innovations

We consider estimation of a missing value for a stationary autoregressive process of order one with exponential innovations and compare two methods of estimation of the missing value, with respect to Pitman's measure of closeness (PMC).     

An empirical study on technical analysis: GARCH (1, 1) model

One of the deficits of the common Bollinger band is that it fails to consider the fat tails/leptokurtosis often exists in financial time series. An adjusted Bollinger band generated by rolling GARCH regression method is proposed in this study. The performance of the adjusted Bollinger band strategy on EUR, GBP, JPY, and AUD vs. USD foreign exchange trading is evaluated. Results show that in general, the adjusted Bollinger band performs better than the traditional one in terms of success ratios, net successes, and profit. In addition, no matter there is transaction cost or not, only adjusted Bollinger strategies are recommended for investors. Adjusted Bollinger band strategies with MA 5 or 10 are recommended for EUR, GBP, and JPY. Adjusted Bollinger strategy with MA 20 is the recommended strategies for AUD.     

A First-Order Spatial Integer-Valued Autoregressive SINAR(1, 1) Model

Binomial thinning operator has a major role in modeling one-dimensional integer-valued autoregressive time series models. The purpose of this article is to extend the use of such operator to define a new stationary first-order spatial non negative, integer-valued autoregressive SINAR(1, 1) model. We study some properties of this model like the mean, variance and autocorrelation function. Yule-Walker estimator of the model parameters is also obtained. Some numerical results of the model are presented and, moreover, this model is applied to a real data set.     

A variable-neighbourhood search algorithm for finding optimal run orders in the presence of serial correlation

The responses obtained from response surface designs that are run sequentially often exhibit serial correlation or time trends. The order in which the runs of the design are performed then has an impact on the precision of the parameter estimators. This article proposes the use of a variable-neighbourhood search algorithm to compute run orders that guarantee a precise estimation of the effects of the experimental factors. The importance of using good run orders is demonstrated by seeking D-optimal run orders for a central composite design in the presence of an AR(1) autocorrelation pattern.     

Zero-Inflated NGINAR(1) process

In this paper, we develop a zero-inflated NGINAR(1) process as an alternative to the NGINAR(1) process (Risti?, Nasti?, and Bakouch 2009 Risti?, M. M., A. S. Nasti?, and H. S. Bakouch. 2009. A new geometric first-order integer-valued autoregressive (NGINAR(1)) process. Journal of Statistical Planning and Inference 139:2218–26.[Crossref], [Web of Science ®] , [Google Scholar]) when the number of zeros in the data is larger than the expected number of zeros by the geometric process. The proposed process has zero-inflated geometric marginals and contains the NGINAR(1) process as a particular case. In addition, various properties of the new process are derived such as conditional distribution and autocorrelation structure. Yule-Walker, probability based Yule-Walker, conditional least squares and conditional maximum likelihood estimators of the model parameters are derived. An extensive Monte Carlo experiment is conducted to evaluate the performances of these estimators in finite samples. Forecasting performances of the model are discussed. Application to a real data set shows the flexibility and potentiality of the new model.