On the Computation and Finite Sample Behavior of the Constrained M-Estimates for Multivariate Location and Scatter

Abstract

Constrained M (CM) estimates of multivariate location and scatter [Kent, J. T., Tyler, D. E. (1996). Constrained M-estimation for multivariate location and scatter. Ann. Statist. 24:1346–1370] are defined as the global minimum of an objective function subject to a constraint. These estimates combine the good global robustness properties of the S estimates and the good local robustness properties of the redescending M estimates. The CM estimates are not explicitly defined. Numerical methods have to be used to compute the CM estimates. In this paper, we give an algorithm to compute the CM estimates. Using the algorithm, we give a small simulation study to demonstrate the capability of the algorithm finding the CM estimates, and also to explore the finite sample behavior of the CM estimates. We also use the CM estimators to estimate the location and scatter parameters of some multivariate data sets to see the performance of the CM estimates dealing with the real data sets that may contain outliers.     

The S-estimator in the change-point random model with long memory

This paper considers two-phase random design linear regression models. Errors and regressors are stationary long-range-dependent Gaussian processes. The regression parameters, the scale parameter and the change-point are estimated using a method introduced by Rousseeuw and Yohai [Robust regression by means of S-estimators, in Robust and Nonlinear Time Series Analysis, J. Franke, W. Hrdle, and R.D. Martin, eds., Lecture Notes in Statistics, Vol. 26, Springer, New York, 1984, pp. 256–272], which is called the S-estimator and has the property be more robust than the classical estimators in the sense that the outliers do not bias the estimation results. Some asymptotic results, including the strong consistency and the convergence rate of the S-estimator are proved. Simulations and an application to the Nile River data are also presented. It is shown via Monte Carlo simulations that the S-estimator is better than two other estimators that are proposed in the literature.     

Robust estimation of the mean vector for high-dimensional data set using robust clustering

The first step in statistical analysis is the parameter estimation. In multivariate analysis, one of the parameters of interest to be estimated is the mean vector. In multivariate statistical analysis, it is usually assumed that the data come from a multivariate normal distribution. In this situation, the maximum likelihood estimator (MLE), that is, the sample mean vector, is the best estimator. However, when outliers exist in the data, the use of sample mean vector will result in poor estimation. So, other estimators which are robust to the existence of outliers should be used. The most popular robust multivariate estimator for estimating the mean vector is S-estimator with desirable properties. However, computing this estimator requires the use of a robust estimate of mean vector as a starting point. Usually minimum volume ellipsoid (MVE) is used as a starting point in computing S-estimator. For high-dimensional data computing, the MVE takes too much time. In some cases, this time is so large that the existing computers cannot perform the computation. In addition to the computation time, for high-dimensional data set the MVE method is not precise. In this paper, a robust starting point for S-estimator based on robust clustering is proposed which could be used for estimating the mean vector of the high-dimensional data. The performance of the proposed estimator in the presence of outliers is studied and the results indicate that the proposed estimator performs precisely and much better than some of the existing robust estimators for high-dimensional data.     

Almost sure representation for multivariate m-estimators

We generalize Carroll (1978)'s almost sure expansion for one-dimensional M-estimators of location and scale to multivariate M-estimators with the same error order. Bahadur type representations for M-estimators and S-estimators of location and scatter are given as an application. In this case, the error rate is a slight improvement over that of He and Wang (1993).     

A small-sample correction factor for S-estimators

S-estimators are frequently used as robust estimators of regression and of location and dispersion. Under certain differentiability conditions, S-estimators of multivariate location and dispersion parameters are consistent [Davies PL. Asymtotic behaviour of S-estimators of multivariate location parameters and dispersion matrices. Ann Stat. 1987;15(3):1269–1292]. However, it has been observed that the S-estimators of dispersion parameters give biased results in the case of small-sample data sets. In this work, we constructed formulas based on simulation studies, which allow us to compute small-sample correction factors for all sample sizes and dimensions for S-estimators of dispersion parameters without having to carry out any new simulations. We considered real data to illustrate the effects of the small-sample correction factor.     

Robust estimation of complicated profiles using wavelets

Some quality characteristics are well defined when treated as the response variables and their relationships are identified to some independent variables. This relationship is called a profile. The parametric models, such as linear models, may be used to model the profiles. However, due to the complexity of many processes in practical applications, it is inappropriate to model the process using parametric models. In these cases non parametric methods are used to model the processes. One of the most applicable non parametric methods used to model complicated profiles is the wavelet. Many authors considered the use of the wavelet transformation only for monitoring the processes in phase II. The problem of estimating the in-control profile in phase I using wavelet transformation is not deeply addressed. Usually classical estimators are used in phase I to estimate the in-control profiles, even when the wavelet transformation is used. These estimators are suitable if the data do not contain outliers. However, when the outliers exist, these estimators cannot estimate the in-control profile properly. In this research, a robust method of estimating the in-control profiles is proposed, which is insensitive to the presence of outliers and could be applied when the wavelet transformation is used. The proposed estimator is the combination of the robust clustering and the S-estimator. This estimator is compared with the classical estimator of the in-control profile in the presence of outliers. The results from a large simulation study show that using the proposed method, one can estimate the in-control profile precisely when the data are contaminated either locally or globally.     

S-estimator in partially linear regression models

In this paper, a robust estimator is proposed for partially linear regression models. We first estimate the nonparametric component using the penalized regression spline, then we construct an estimator of parametric component by using robust S-estimator. We propose an iterative algorithm to solve the proposed optimization problem, and introduce a robust generalized cross-validation to select the penalized parameter. Simulation studies and a real data analysis illustrate that the our proposed method is robust against outliers in the dataset or errors with heavy tails.