{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,31]],"date-time":"2026-10-31T14:58:49Z","timestamp":1793458729988,"version":"4.5.0"},"reference-count":48,"publisher":"Oxford University Press (OUP)","issue":"3","license":[{"start":{"date-parts":[[2021,10,13]],"date-time":"2021-10-13T00:00:00Z","timestamp":1634083200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2022,8,24]]},"abstract":"<jats:title>Summary<\/jats:title>\n               <jats:p>Lag windows are commonly used in time series analysis, econometrics, steady-state simulation and Markov chain Monte Carlo to estimate time-average covariance matrices. In the presence of positive correlation in the underlying process, estimators of this matrix almost always exhibit significant negative bias, leading to undesirable finite-sample properties. We propose a new family of lag windows specifically designed to improve finite-sample performance by offsetting this negative bias. Any existing lag window can be adapted into a lugsail equivalent with no additional assumptions. We use these lag windows in spectral variance estimators and demonstrate their advantages in a linear regression model with autocorrelated and heteroskedastic residuals. We further employ the lugsail lag windows in weighted batch means estimators because of their computational efficiency on large simulation output. We obtain bias and variance results for these multivariate estimators and significantly weaken the mixing condition on the process. Superior finite-sample properties are demonstrated in a vector autoregressive process and a Bayesian logistic regression model.<\/jats:p>","DOI":"10.1093\/biomet\/asab049","type":"journal-article","created":{"date-parts":[[2021,10,10]],"date-time":"2021-10-10T11:06:38Z","timestamp":1633863998000},"page":"735-750","source":"Crossref","is-referenced-by-count":16,"title":["Lugsail lag windows for estimating time-average covariance matrices"],"prefix":"10.1093","volume":"109","author":[{"given":"D","family":"Vats","sequence":"first","affiliation":[{"name":"Department of Mathematics and Statistics, Indian Institute of Technology Kanpur , Kanpur 208016, India"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"J M","family":"Flegal","sequence":"additional","affiliation":[{"name":"Department of Statistics, University of California Riverside , 1428 Olmsted Hall, Riverside, California 92521, U.S.A"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2021,10,13]]},"reference":[{"key":"2022082420281065300_B1","doi-asserted-by":"crossref","first-page":"397","DOI":"10.1002\/nav.20215","article-title":"Exact expected values of variance estimators for simulation","volume":"54","author":"Aktaran-Kalayc\u0131,","year":"2007","journal-title":"Naval Res. Logist."},{"key":"2022082420281065300_B2","volume-title":"The Statistical Analysis of Time Series.","author":"Anderson,","year":"1971"},{"key":"2022082420281065300_B3","doi-asserted-by":"crossref","first-page":"817","DOI":"10.2307\/2938229","article-title":"Heteroskedasticity and autocorrelation consistent covariance matrix estimation","volume":"59","author":"Andrews,","year":"1991","journal-title":"Econometrica"},{"key":"2022082420281065300_B4","doi-asserted-by":"crossref","first-page":"477","DOI":"10.1007\/s10463-007-0154-0","article-title":"Higher-order accurate polyspectral estimation with flat-top lag-windows","volume":"61","author":"Berg,","year":"2009","journal-title":"Ann. Inst. Statist. Math."},{"key":"2022082420281065300_B5","doi-asserted-by":"crossref","first-page":"540","DOI":"10.1109\/5.135378","article-title":"Estimating and interpreting the instantaneous frequency of a signal. II. Algorithms and applications","volume":"80","author":"Boashash,","year":"1992","journal-title":"Proc. IEEE"},{"key":"2022082420281065300_B6","doi-asserted-by":"crossref","first-page":"1076","DOI":"10.1080\/01621459.2016.1189337","article-title":"Automatic optimal batch size selection for recursive estimators of time-average covariance matrix","volume":"112","author":"Chan,","year":"2017","journal-title":"J. Am. Statist. Assoc."},{"key":"2022082420281065300_B7","doi-asserted-by":"crossref","first-page":"866","DOI":"10.1111\/sjos.12279","article-title":"High-order corrected estimator of asymptotic variance with optimal bandwidth","volume":"44","author":"Chan,","year":"2017","journal-title":"Scand. J. Statist."},{"key":"2022082420281065300_B8","first-page":"302","article-title":"Multivariate inference in stationary simulation using batch means","volume-title":"Proc. 19th Conf. Winter Simulation","author":"Chen,","year":"1987"},{"key":"2022082420281065300_B9","doi-asserted-by":"crossref","first-page":"1288","DOI":"10.1287\/mnsc.43.9.1288","article-title":"Large-sample results for batch means","volume":"43","author":"Chien,","year":"1997","journal-title":"Manag. Sci."},{"key":"2022082420281065300_B10","doi-asserted-by":"crossref","first-page":"184","DOI":"10.1016\/j.jmva.2017.05.009","article-title":"Multivariate initial sequence estimators in Markov chain Monte Carlo","volume":"159","author":"Dai,","year":"2017","journal-title":"J. Mult. Anal."},{"key":"2022082420281065300_B11","doi-asserted-by":"crossref","first-page":"282","DOI":"10.1287\/opre.43.2.282","article-title":"Mean-square consistency of the variance estimator in steady-state simulation output analysis","volume":"43","author":"Damerdji,","year":"1995","journal-title":"Oper. Res."},{"key":"2022082420281065300_B12","doi-asserted-by":"crossref","first-page":"262","DOI":"10.1017\/S0266466600162061","article-title":"A strong consistency proof for heteroskedasticity and autocorrelation consistent covariance matrix estimators","volume":"16","author":"De Jong,","year":"2000","journal-title":"Economet. Theory"},{"key":"2022082420281065300_B13","doi-asserted-by":"publisher","DOI":"10.3386\/t0195","article-title":"Inferences from parametric and non-parametric covariance matrix estimation procedures","volume-title":"Technical Working Paper 0195, National Bureau of Economic Research","author":"den Haan,","year":"1996"},{"key":"2022082420281065300_B14","first-page":"299","article-title":"A practitioner\u2019s guide to robust covariance matrix estimation","volume-title":"Handbook of Statistics","author":"den Haan,","year":"1997"},{"key":"2022082420281065300_B15","doi-asserted-by":"publisher","DOI":"10.3386\/t0255","article-title":"Robust covariance matrix estimation with data-dependent VAR prewhitening order","volume-title":"Technical Working Paper 0255, National Bureau of Economic Research","author":"den Haan,","year":"2000"},{"key":"2022082420281065300_B16","doi-asserted-by":"crossref","first-page":"459","DOI":"10.1109\/WSC.2015.7408187","article-title":"Jackknifed variance estimators for simulation output analysis","volume-title":"2015 Winter Simulation Conference (WSC)","author":"Dinge\u00e7,","year":"2015"},{"key":"2022082420281065300_B17","doi-asserted-by":"crossref","first-page":"1034","DOI":"10.1214\/09-AOS735","article-title":"Batch means and spectral variance estimators in Markov chain Monte Carlo","volume":"38","author":"Flegal,","year":"2010","journal-title":"Ann. Statist."},{"key":"2022082420281065300_B18","first-page":"473","article-title":"Practical Markov chain Monte Carlo (with Discussion)","volume":"7","author":"Geyer,","year":"1992","journal-title":"Statist. Sci."},{"key":"2022082420281065300_B19","doi-asserted-by":"crossref","first-page":"180","DOI":"10.1214\/aoap\/1177005777","article-title":"The asymptotic validity of sequential stopping rules for stochastic simulations","volume":"2","author":"Glynn,","year":"1992","journal-title":"Ann. Appl. Prob."},{"key":"2022082420281065300_B20","first-page":"252","article-title":"On the optimal window shape for genomic signal processing","volume-title":"2008 Int. Conf. Computer and Communication Engineering","author":"Gunawan,","year":"2008"},{"key":"2022082420281065300_B21","doi-asserted-by":"crossref","DOI":"10.1002\/9780470316429","volume-title":"Multiple Time Series","author":"Hannan,","year":"1970"},{"key":"2022082420281065300_B22","doi-asserted-by":"crossref","first-page":"55","DOI":"10.1016\/j.csda.2017.09.007","article-title":"Alternative HAC covariance matrix estimators with improved finite sample properties","volume":"119","author":"Hartigan,","year":"2018","journal-title":"Comp. Statist. Data Anal."},{"key":"2022082420281065300_B23","doi-asserted-by":"crossref","first-page":"1117","DOI":"10.1214\/14-AOS1301","article-title":"Covariance matrix estimation and linear process bootstrap for multivariate time series of possibly increasing dimension","volume":"43","author":"Jentsch,","year":"2015","journal-title":"Ann. Statist."},{"key":"2022082420281065300_B24","doi-asserted-by":"crossref","first-page":"299","DOI":"10.1214\/154957804100000051","article-title":"On the Markov chain central limit theorem","volume":"1","author":"Jones,","year":"2004","journal-title":"Prob. Surv."},{"key":"2022082420281065300_B25","doi-asserted-by":"crossref","first-page":"2093","DOI":"10.1111\/1468-0262.00366","article-title":"Heteroskedasticity-autocorrelation robust standard errors using the Bartlett kernel without truncation","volume":"70","author":"Kiefer,","year":"2002","journal-title":"Econometrica"},{"key":"2022082420281065300_B26","doi-asserted-by":"crossref","first-page":"1350","DOI":"10.1017\/S026646660218604X","article-title":"Heteroskedasticity-autocorrelation robust testing using bandwidth equal to sample size","volume":"18","author":"Kiefer,","year":"2002","journal-title":"Economet. Theory"},{"key":"2022082420281065300_B27","doi-asserted-by":"crossref","first-page":"1130","DOI":"10.1017\/S0266466605050565","article-title":"A new asymptotic theory for heteroskedasticity-autocorrelation robust tests","volume":"21","author":"Kiefer,","year":"2005","journal-title":"Economet. Theory"},{"key":"2022082420281065300_B28","doi-asserted-by":"crossref","first-page":"85","DOI":"10.1016\/S0167-7152(99)00090-5","article-title":"Monte Carlo error estimation for multivariate Markov chains","volume":"46","author":"Kosorok,","year":"2000","journal-title":"Statist. Prob. Lett."},{"key":"2022082420281065300_B29","doi-asserted-by":"crossref","first-page":"1003","DOI":"10.1214\/aop\/1176994565","article-title":"Almost sure invariance principles for partial sums of mixing B-valued random variables","volume":"8","author":"Kuelbs,","year":"1980","journal-title":"Ann. Prob."},{"key":"2022082420281065300_B30","doi-asserted-by":"crossref","first-page":"541","DOI":"10.1080\/07350015.2018.1506926","article-title":"HAR inference: Recommendations for practice","volume":"36","author":"Lazarus,","year":"2018","journal-title":"J. Bus. Econ. Statist."},{"key":"2022082420281065300_B31","doi-asserted-by":"crossref","first-page":"3397","DOI":"10.1214\/18-EJS1483","article-title":"Weighted batch means estimators in Markov chain Monte Carlo","volume":"12","author":"Liu,","year":"2018","journal-title":"Electron. J. Statist."},{"key":"2022082420281065300_B32","doi-asserted-by":"crossref","first-page":"1331","DOI":"10.1016\/j.jeconom.2007.01.019","article-title":"A theory of robust long-run variance estimation","volume":"141","author":"M\u00fcller,","year":"2007","journal-title":"J. Economet."},{"key":"2022082420281065300_B33","doi-asserted-by":"crossref","first-page":"311","DOI":"10.1080\/07350015.2014.931238","article-title":"HAC corrections for strongly autocorrelated time series","volume":"32","author":"M\u00fcller,","year":"2014","journal-title":"J. Bus. Econ. Statist."},{"key":"2022082420281065300_B34","doi-asserted-by":"crossref","first-page":"703","DOI":"10.2307\/1913610","article-title":"A simple, positive semi-definite, heteroskedasticity and autocorrelation consistent covariance matrix","volume":"55","author":"Newey,","year":"1987","journal-title":"Econometrica"},{"key":"2022082420281065300_B35","doi-asserted-by":"crossref","first-page":"329","DOI":"10.1214\/aoms\/1177706962","article-title":"On consistent estimates of the spectrum of a stationary time series","volume":"28","author":"Parzen,","year":"1957","journal-title":"Ann. Math. Statist."},{"key":"2022082420281065300_B36","doi-asserted-by":"crossref","first-page":"837","DOI":"10.1111\/j.1468-2354.2006.00398.x","article-title":"Spectral density estimation and robust hypothesis testing using steep origin kernels without truncation","volume":"47","author":"Phillips,","year":"2006","journal-title":"Int. Econ. Rev."},{"key":"2022082420281065300_B37","doi-asserted-by":"crossref","first-page":"67","DOI":"10.1111\/j.1467-9892.1995.tb00223.x","article-title":"Bias-corrected nonparametric spectral estimation","volume":"16","author":"Politis,","year":"1995","journal-title":"J. Time Ser. Anal."},{"key":"2022082420281065300_B38","doi-asserted-by":"crossref","first-page":"41","DOI":"10.1016\/0378-3758(95)00069-0","article-title":"On flat-top kernel spectral density estimators for homogeneous random fields","volume":"51","author":"Politis,","year":"1996","journal-title":"J. Statist. Plan. Infer."},{"key":"2022082420281065300_B39","volume-title":"Spectral Analysis and Time Series","author":"Priestley,","year":"1981"},{"key":"2022082420281065300_B40","year":"2022","journal-title":"R: A Language and Environment for Statistical Computing"},{"key":"2022082420281065300_B41","first-page":"3","article-title":"The relative importance of bias and variability in the estimation of the variance of a statistic","volume":"42","author":"Simonoff,","year":"1993","journal-title":"J. R. Statist. Soc. D"},{"key":"2022082420281065300_B42","doi-asserted-by":"crossref","first-page":"110","DOI":"10.1287\/mnsc.41.1.110","article-title":"Optimal mean-squared-error batch sizes","volume":"41","author":"Song,","year":"1995","journal-title":"Manag. Sci."},{"key":"2022082420281065300_B43","doi-asserted-by":"crossref","first-page":"175","DOI":"10.1111\/j.0012-9682.2008.00822.x","article-title":"Optimal bandwidth selection in heteroskedasticity\u2013autocorrelation robust testing","volume":"76","author":"Sun,","year":"2008","journal-title":"Econometrica"},{"key":"2022082420281065300_B44","doi-asserted-by":"crossref","first-page":"587","DOI":"10.2307\/1427459","article-title":"Non-linear time series and Markov chains","volume":"22","author":"Tj\u00f8stheim,","year":"1990","journal-title":"Adv. Appl. Prob."},{"key":"2022082420281065300_B45","doi-asserted-by":"crossref","first-page":"1860","DOI":"10.3150\/16-BEJ914","article-title":"Strong consistency of multivariate spectral variance estimators in Markov chain Monte Carlo","volume":"24","author":"Vats,","year":"2018","journal-title":"Bernoulli"},{"key":"2022082420281065300_B46","doi-asserted-by":"crossref","first-page":"321","DOI":"10.1093\/biomet\/asz002","article-title":"Multivariate output analysis for Markov chain Monte Carlo","volume":"106","author":"Vats,","year":"2019","journal-title":"Biometrika"},{"key":"2022082420281065300_B47","doi-asserted-by":"crossref","first-page":"817","DOI":"10.2307\/1912934","article-title":"A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity","volume":"48","author":"White,","year":"1980","journal-title":"Econometrica"},{"key":"2022082420281065300_B48","doi-asserted-by":"crossref","first-page":"1054","DOI":"10.1017\/S026646661400067X","article-title":"Optimal bandwidth selection for robust generalized method of moments estimation","volume":"31","author":"Wilhelm,","year":"2015","journal-title":"Economet. Theory"}],"container-title":["Biometrika"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/biomet\/advance-article-pdf\/doi\/10.1093\/biomet\/asab049\/41987583\/asab049.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biomet\/article-pdf\/109\/3\/735\/45512212\/asab049.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/biomet\/article-pdf\/109\/3\/735\/45512212\/asab049.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,8,24]],"date-time":"2022-08-24T20:29:28Z","timestamp":1661372968000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/biomet\/article\/109\/3\/735\/6395353"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,10,13]]},"references-count":48,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2021,10,13]]},"published-print":{"date-parts":[[2022,8,24]]}},"URL":"https:\/\/doi.org\/10.1093\/biomet\/asab049","relation":{},"ISSN":["0006-3444","1464-3510"],"issn-type":[{"value":"0006-3444","type":"print"},{"value":"1464-3510","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2022,9,1]]},"published":{"date-parts":[[2021,10,13]]}}}