ON THE RELATIONSHIP BETWEEN GENERALIZED LEAST SQUARES AND GAUSSIAN ESTIMATION OF VECTOR ARMA MODELS

Abstract. This paper investigates theoretical aspects of the relationship between the generalized least squares and Gaussian estimation schemes for vector autoregressive moving-average models. The asymptotic convergence of the generalized least squares estimator to the Gaussian estimator is established and an alternative numerical method for implementing the generalized least squares scheme is proposed. Finally, some simulation results are presented to illustrate the theory.     

Bootstrapping Time Series Regressions with Integrated Processes

This paper studies the bootstrap procedures for time series regressions with integrated processes. Both estimation and hypothesis testing are studied. It is shown that the suggested bootstrap approximations to the distribution of the least squares estimator and the regression test statistic are asymptotically valid. A Monte Carlo experiment is conducted to evaluate the finite sample performance of these bootstrap procedures. The simulation results indicate that the bootstrap method provides reasonably good approximation to the distribution of the least squares estimator, and gives proper size and satisfactory power.     

ALTERNATIVE ESTIMATORS AND UNIT ROOT TESTS FOR THE AUTOREGRESSIVE PROCESS

Abstract. We compare several estimators for the second-order autoregressive process and compare the associated tests for a unit root. Monte Carlo results are reported for the ordinary least squares estimator, the simple symmetric least squares estimator and the weighted symmetric least squares estimator. The weighted symmetric least squares estimator of the autoregressive parameters generally has smaller mean square error than that of the ordinary least squares estimator, particularly when one root is close to one in absolute value. For the second-order model with known zero intercept, the one-sided ordinary least squares test for a unit root is more powerful than the symmetric tests. For the model with an estimated intercept, the one-sided weighted symmetric least squares test is the most powerful test.     

Least squares estimation of ARCH models with missing observations

A least squares estimator for ARCH models in the presence of missing data is proposed. Strong consistency and asymptotic normality are derived. Monte Carlo simulation results are analysed and an application to real data of a Chilean stock index is reported.     

ESTIMATION OF COEFFICIENTS OF TIME SERIES REGRESSION WITH A NONSTATIONARY ERROR PROCESS

Abstract. We treat a problem of estimating unknown coefficients of a time series regression when the variance of the error changes with time, i.e. when a process which the error term obeys is nonstationary. First, we show the weak consistency of the ordinary least squares estimator for the coefficients of a polynomial regression under some assumptions on the covariance structure of the error process. Next, we propose a nonparametric method for estimating the variance of the error process and a weighted least squares estimator of the regression coefficients, which is constructed by using the estimator of the variance. We investigate statistical properties of our proposed estimator in the following way. We consider the prediction of a future value of a linear trend by using our proposed estimator and evaluate its prediction error. By simulation studies, we compare the prediction error of the predictor constructed by using our proposed estimator with the prediction errors obtained for other estimators including the ordinary least squares estimator when the variance of the error process increases with time and the sample sizes are small. As a result, our proposed estimator seems to be reasonable.     

Using least squares to generate forecasts in regressions with serial correlation

Abstract. The topic of serial correlation in regression models has attracted a great deal of research in the last 50 years. Most of these studies have assumed that the structure of the error covariance matrix Ω was known or could be consistently estimated from the data. In this article, we describe a new procedure for generating forecasts for regression models with serial correlation based on ordinary least squares and on an approximate representation of the form of the autocorrelation. We prove that the predictors from this specification are asymtotically efficient under some regularity conditions. In addition, we show that there is not much to be gained in trying to identify the correct form of the serial correlation since efficient forecasts can be generated using autoregressive approximations of the autocorrelation. A large simulation study is also used to compare the finite sample predictive efficiencies of this new estimator vis‐à‐vis estimators based on ordinary least squares and generalized least squares.     

ESTIMATION OF LINEAR REGRESSION MODEL WITH AUTOCORRELATED DISTURBANCES

Abstract. In this paper we have derived the large sample asymptotic approximation for the variance-covariance matrix of the two stage Prais-Winston estimator of the regression coefficients. The efficiency properties of this estimator with respect to ordinary least squares, and generalized least squares with a known autocorrelation coefficient are then analysed numerically. The results are useful for the practitioners dealing with moderate size sample data.     

Explosive Random‐Coefficient AR(1) Processes and Related Asymptotics for Least‐Squares Estimation

Abstract. Large sample properties of the least‐squares and weighted least‐squares estimates of the autoregressive parameter of the explosive random‐coefficient AR(1) process are discussed. It is shown that, contrary to the standard AR(1) case, the least‐squares estimator is inconsistent whereas the weighted least‐squares estimator is consistent and asymptotically normal even when the error process is not necessarily Gaussian. Conditional asymptotics on the event that a certain limiting random variable is non‐zero is also discussed.     

TESTING FOR A UNIT ROOT IN AN AR(1) TIME SERIES USING IRREGULARLY OBSERVED DATA

Abstract. For an AR(1) model having a unit root with nonconsecutively observed or missing data we consider the ordinary least squares estimator, the one-step Newton-Raphson estimator and an ordinary least squares type estimator which is a simple approximation of the Newton-Raphson estimator. It is shown that the limiting distributions of these estimators of the unit root are the same as those of the regression estimators as tabulated by Dickey and Fuller (Distribution of the estimators for autoregressive time series with a unit root. J. Am. Statist. Assoc. 74 (1979), 427–31) for the complete data situation. Simulation results show that our proposed unit root tests perform very well for small samples.     

THE ESTIMATION AND APPLICATION OF LONG MEMORY TIME SERIES MODELS

Abstract. The definitions of fractional Gaussian noise and integrated (or fractionally differenced) series are generalized, and it is shown that the two concepts are equivalent. A new estimator of the long memory parameter in these models is proposed, based on the simple linear regression of the log periodogram on a deterministic regressor. The estimator is the ordinary least squares estimator of the slope parameter in this regression, formed using only the lowest frequency ordinates of the log periodogram. Its asymptotic distribution is derived, from which it is evident that the conventional interpretation of these least squares statistics is justified in large samples. Using synthetic data the asymptotic theory proves to be reliable in samples of 50 observations or more. For three postwar monthly economic time series, the estimated integrated series model provides more reliable out-of-sample forecasts than do more conventional procedures.     

基于岭估计的化工过程测量置信水平研究(英文)

Ordinary least squares （OLS） algorithm is widely applied in process measurement, because the sensor model used to estimate unknown parameters can be approximated through multivariate linear model. However, with few or noisy data or multi-collinearity, unbiased OLS leads to large variance. Biased estimators, especially ridge es-timator, have been introduced to improve OLS by trading bias for variance. Ridge estimator is feasible as an esti-mator with smaller variance. At the same confidence level, with additive noise as the normal random variable, the less variance one estimator has, the shorter the two-sided symmetric confidence interval is. However, this finding is limited to the unbiased estimator and few studies analyze and compare the confidence levels between ridge estima-tor and OLS. This paper derives the matrix of ridge parameters under necessary and sufficient conditions based on which ridge estimator is superior to OLS in terms of mean squares error matrix, rather than mean squares error. Then the confidence levels between ridge estimator and OLS are compared under the condition of OLS fixed sym-metric confidence interval, rather than the criteria for evaluating the validity of different unbiased estimators. We conclude that the confidence level of ridge estimator can not be directly compared with that of OLS based on the criteria available for unbiased estimators, which is verified by a simulation and a laboratory scale experiment on a single parameter measurement.     

Practical improvements to autocovariance least‐squares

Identifying disturbance covariances from data is a critical step in estimator design and controller performance monitoring. Here, the autocovariance least‐squares (ALS) method for this identification is examined. For large industrial models with poorly observable states, the process noise covariance is high dimensional and the optimization problem is poorly conditioned. Also, weighting the least‐squares problem with the identity matrix does not provide minimum variance estimates. Here, ALS method to resolve these two challenges is modified. Poorly observable states using the singular value decomposition (SVD) of the observability matrix is identified and removed, thus decreasing the computational time. Using a new feasible‐generalized least‐squares estimator that approximates the optimal weighting from data, the variance of the estimates is significantly reduced. The new approach on industrial data sets provided by Praxair is successfully demonstrated. The disturbance model identified by the ALS method produces an estimator that performs optimally over a year‐long period. © 2015 American Institute of Chemical Engineers AIChE J, 61: 1840–1855, 2015     

Finite-sample comparison of robust estimators for nonlinear regression using Monte Carlo simulation: Part I. Univariate response models

Classical least squares can be strongly affected due to the inevitable occurrence of departures from its model assumptions, most notably those from the distributional assumptions. Robust estimators, on the other hand, will resist them. Unfortunately, the multiplicity of alternative robust regression estimators that have been suggested in the literature over the years is a source of confusion for practitioners of regression analysis. Moreover, little is known about their small-sample performance in the nonlinear regression setting, in particular on the chemical engineering field. A simulation study comparing six such estimators (namely LMS, LTS, LTD, MM-, τ-, and L_p-norm) together with the usual least squares estimator is presented. The results obtained provide guidance as to the choice of an appropriate estimator.     

ESTIMATING THE ARCH PARAMETERS BY SOLVING LINEAR EQUATIONS

Abstract. This paper discusses the asymptotics of two-stage least squares estimator of the parameters of ARCH models. The estimator is easy to obtain since it involves solving two sets of linear equations. At the same time, the estimator has the same asymptotic efficiency as that of the widely used quasi-maximum likelihood estimator. Simulation results show that, even for small sample size, the performance of our estimator compared to the quasi-maximum likelihood estimator is better.     

Computer Algebra Derivation of the Bias of Linear Estimators of Autoregressive Models

Abstract. A symbolic method which can be used to obtain the asymptotic bias and variance coefficients to order O(1/n) for estimators in stationary time series is discussed. Using this method, the large‐sample bias of the Burg estimator in the AR(p) for p = 1, 2, 3 is shown to be equal to that of the least squares estimators in both the known and unknown mean cases. Previous researchers have only been able to obtain simulation results for the Burg estimator's bias because this problem is too intractable without using computer algebra. The asymptotic bias coefficient to O(1/n) of Yule–Walker as well as least squares estimates is also derived in AR(3) models. Our asymptotic results show that for the AR(3), just as in the AR(2), the Yule–Walker estimates have a large bias when the parameters are near the nonstationary boundary. The least squares and Burg estimates are much better in this situation. Simulation results confirm our findings.     

APPROXIMATE DISTRIBUTION OF PARAMETER ESTIMATORS FOR FIRST-ORDER AUTOREGRESSIVE MODELS

Abstract. Formulae for the exact bias and mean square error for the least squares for forward-backward least squares estimators are obtained based on the explicit expressions for the moment-generating and characteristic functions of quadratic form in the first-order autoregressive process. Asymptotic expressions for their cumulants and the maximum likelihood estimator are given. Approximations of the distributions of the above estimators are proposed based on the Ornstein-Ulenbeck process. A simple computational procedure for the exact distribution is developed, and some numerical comparisons are given which support the overall good accuracy of the approximation and confirm that the maximum likelihood estimator performs better than the others.     

A Note on Mean-squared Prediction Errors of the Least Squares Predictors in Random Walk Models

An asymptotic expression for the mean-squared prediction error (MSPE) of the least squares predictor is obtained in the random walk model. It is shown that the term of order 1/ n in this error, where n is the sample size, is twice as large as the one obtained from the first-order autoregressive (AR(1)) model satisfying the stationary assumption. Moreover, while the correlation between the squares of the (normalized) regressor variable and normalized least squares estimator is asymptotically negligible in the stationary AR(1) model, we have found that the correlation has significantly negative value in the random walk model. To obtain these results, a new methodology, which is found to be useful in dealing with the moment properties of a strongly dependent process, is introduced.     

A refined efficiency rate for ordinary least squares and generalized least squares estimators for a linear trend with autoregressive errors

When a straight line is fitted to time series data, generalized least squares (GLS) estimators of the trend slope and intercept are attractive as they are unbiased and of minimum variance. However, computing GLS estimators is laborious as their form depends on the autocovariances of the regression errors. On the other hand, ordinary least squares (OLS) estimators are easy to compute and do not involve the error autocovariance structure. It has been known for 50 years that OLS and GLS estimators have the same asymptotic variance when the errors are second‐order stationary. Hence, little precision is gained by using GLS estimators in stationary error settings. This article revisits this classical issue, deriving explicit expressions for the GLS estimators and their variances when the regression errors are drawn from an autoregressive process. These expressions are used to show that OLS methods are even more efficient than previously thought. Specifically, we show that the convergence rate of variance differences is one polynomial degree higher than that of least squares estimator variances. We also refine Grenander's (1954) variance ratio. An example is presented where our new rates cannot be improved upon. Simulations show that the results change little when the autoregressive parameters are estimated.     

SOME PROPERTIES OF THE MAXIMUM LIKELIHOOD ESTIMATOR IN THE SIMULTANEOUS SWITCHING AUTOREGRESSIVE MODEL

Abstract. The simultaneous switching autoregressive (SSAR) model proposed by Kunitomo and Sato (A non-linearity in economic time series and disequilibrium econometric models. In Theory and Application of Mathematical Statistics (ed. A. Takemura). Tokyo:University of Tokyo Press (in Japanese), 1994; Asymmetry in economic time series and simultaneous switching autoregressive model. Struct. Change Econ. Dyn. , forthcoming (1994).) is a Markovian non-linear time series model. We investigate the finite sample as well as the asymptotic properties of the least squares estimator and the maximum likelihood (ML) estimator. Due to a specific simultaneity involved in the SSAR model, the least squares estimator is badly biased. However, the ML estimator under the assumption of Gaussian disturbances gives reasonable estimates.     

SIGNIFICANT VARIABLE SELECTION AND AUTOREGRESSIVE ORDER DETERMINATION FOR TIME‐SERIES PARTIALLY LINEAR MODELS

This paper is concerned with the regression coefficient and autoregressive order shrinkage and selection via the smoothly clipped absolute deviation (SCAD) penalty for a partially linear model with time‐series errors. By combining the profile semi‐parametric least squares method and SCAD penalty technique, a new penalized estimation for the regression and autoregressive parameters in the model is proposed. We show that the asymptotic property of the resultant estimator is the same as if the order of autoregressive error structure and non‐zero regression coefficients are known in advance, thus achieving the oracle property in the sense of Fan and Li (2001). In addition, based on a prewhitening technique, we construct a two‐stage local linear estimator (TSLLE) for the non‐parametric component. It is shown that the TSLLE is more asymtotically efficient than the one that ignores the autoregressive time‐series error structure. Some simulation studies are conducted to illustrate the finite sample performance of the proposed procedure. An example of application on electricity usage data is also illustrated. Copyright © 2014 Wiley Publishing Ltd