Inconsistency of resampling algorithms for high-breakdown regression estimators and a new algorithm

• 成果类型: Article
• 署名作者: Hawkins, DM; Olive, DJ
• 署名单位: Southern Illinois University System; Southern Illinois University
• 刊物名称: JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION
• ISSN/ISSBN: 0162-1459
• DOI: 10.1198/016214502753479293
• 发表日期: 2002
• 页码: 136-148
• 关键词: multivariate location; squares; asymptotics; depth

摘要

Because high-breakdown estimators (HBEs) are impractical to compute exactly in large samples, approximate algorithms are used. The algorithm generally produces an estimator with a lower consistency rate and breakdown value than the exact theoretical estimator. This discrepancy grows with the sample size, with the implication that huge computations are needed for good approximations in large high-dimension samples, The workhorse for HBEs has been the ''elemental set,'' or ''basic resampling,'' algorithm. This turns out to be completely ineffective in high dimensions with high levels of contaraination. However, enriching it with a ''concentration'' step turns it into a method that can handle even high levels of contamination, provided that the regression outliers are located on random cases. It remains ineffective if the regression outliers are concentrated on high-leverage cases. We focus on the multiple regression problem, but several of the broad conclusions-notably, those of the inadequacy of fixed numbers of elemental starts-are also relevant to multivariate location and dispersion estimation. We introduce a new algorithm-the ''X-cluster'' method-for large high-dimensional multiple regression datasets that are beyond the reach of standard resampling methods. This algorithm departs sharply from current HBE algorithms in that, even at a constant percentage of contamination, it is more effective the larger the sample, making a compelling case for using it in the large-sample situations that current methods serve poorly. A multipronged analysis using both traditional ordinary least squares and L, methods along with newer resistant techniques will often detect departures from the multiple regression model that cannot be detected by any single estimator.