期刊
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-STATISTICAL METHODOLOGY
卷 84, 期 3, 页码 1023-1055出版社
OXFORD UNIV PRESS
DOI: 10.1111/rssb.12540
关键词
changepoint estimation; high-dimensional data; missing data; segmentation; sparsity
资金
- Knox studentship from Trinity College, Cambridge
- EPSRC [EP/P031447/1, EP/N031938/1, EP/T02772X/1]
- ERC [101019498]
- EOX funding from Ecole polytechnique
- European Research Council (ERC) [101019498] Funding Source: European Research Council (ERC)
The proposed method introduces a new way to estimate changepoints in partially observed, high-dimensional time series that undergo a simultaneous change in mean. It demonstrates strong effectiveness and efficiency in both simulated data and real oceanographic dataset analysis.
We propose a new method for changepoint estimation in partially observed, high-dimensional time series that undergo a simultaneous change in mean in a sparse subset of coordinates. Our first methodological contribution is to introduce a 'MissCUSUM' transformation (a generalisation of the popular cumulative sum statistics), that captures the interaction between the signal strength and the level of missingness in each coordinate. In order to borrow strength across the coordinates, we propose to project these MissCUSUM statistics along a direction found as the solution to a penalised optimisation problem tailored to the specific sparsity structure. The changepoint can then be estimated as the location of the peak of the absolute value of the projected univariate series. In a model that allows different missingness probabilities in different component series, we identify that the key interaction between the missingness and the signal is a weighted sum of squares of the signal change in each coordinate, with weights given by the observation probabilities. More specifically, we prove that the angle between the estimated and oracle projection directions, as well as the changepoint location error, are controlled with high probability by the sum of two terms, both involving this weighted sum of squares, and representing the error incurred due to noise and the error due to missingness respectively. A lower bound confirms that our changepoint estimator, which we call MissInspect, is optimal up to a logarithmic factor. The striking effectiveness of the MissInspect methodology is further demonstrated both on simulated data, and on an oceanographic data set covering the Neogene period.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据