Multivariate autoregressive techniques for constructing confidence regions on the mean vector

We develop a method for constructing confidence regions on the mean vectors of multivariate processes that is based on a vector autoregressive VAR representation of the data-generating process. A confidence-region-construction algorithm for a general autoregressive model is given. We establish the asymptotic validity of the confidence-region estimator that is, the exact achievement of nominal coverage probability as the sample size tends to infinity when the output process is a stationary vector autoregressive process of known, finite order. With respect to confidence-region volume, coverage probability, and execution time, we carry out an experimental performance comparison of VAR versus the methods of Bonferroni Batch Means BBM, Multivariate Batch Means MBM, and Multivariate Spectral Analysis SPA. The experimental results indicate that i VAR delivered confidence regions with the smallest volume; ii BBM delivered confidence regions with the largest volume, the highest coverage and the smallest execution time; iii in small samples, all of the methods might yield confidence-region estimators whose coverage differs significantly from the nominal level; and iv in large samples for which the sample autocorrelation function indicates a vector autoregressive dependence structure, VAR is a viable technique for simulation output analysis.

[1]  James Durbin,et al.  The fitting of time series models , 1960 .

[2]  D. B. Preston Spectral Analysis and Time Series , 1983 .

[3]  Jorma Rissanen,et al.  A Predictive Least-Squares Principle , 1986 .

[4]  N. Levinson The Wiener (Root Mean Square) Error Criterion in Filter Design and Prediction , 1946 .

[5]  Bruce W. Schmeiser,et al.  Batch Size Effects in the Analysis of Simulation Output , 1982, Oper. Res..

[6]  Sheldon M. Ross,et al.  Introduction to Probability Models (4th ed.). , 1990 .

[7]  John M. Charnes,et al.  A comparison of confidence region estimators for multivariate simulation output , 1988, WSC '88.

[8]  John M. Charnes,et al.  Power comparisons for the multivariate batch-means method , 1990, 1990 Winter Simulation Conference Proceedings.

[9]  P. Whittle Prediction and Regulation by Linear Least-Square Methods , 1983 .

[10]  R. W. Andrews,et al.  Arma-Based Confidence Intervals for Simulation Output Analysis , 1984 .

[11]  E. J. Hannan,et al.  Multiple time series , 1970 .

[12]  Andrew F. Seila,et al.  Multivariate inference in stationary simulation using batch means , 1987, WSC '87.

[13]  Paul Kabaila,et al.  On confidence regions for the mean of a multivariate time series , 1985 .

[14]  G. S. Fishman Principles of Discrete Event Simulation , 1978 .

[15]  R. Serfling Approximation Theorems of Mathematical Statistics , 1980 .

[16]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[17]  Stephen D. Roberts,et al.  Implementing a portable FORTRAN Uniform (0,1) generator , 1983 .

[18]  Sheldon M. Ross Introduction to Probability Models. , 1995 .

[19]  B. Nelson,et al.  Multivariate batch means and control variates , 1992 .

[20]  Richard A. Davis,et al.  Time Series: Theory and Methods (2nd ed.). , 1992 .

[21]  H. Lütkepohl COMPARISON OF CRITERIA FOR ESTIMATING THE ORDER OF A VECTOR AUTOREGRESSIVE PROCESS , 1985 .

[22]  Richard A. Davis,et al.  Time Series: Theory and Methods , 2013 .