Robust Multivariate Mean Estimation: the Optimality of Trimmed Mean
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Efficient Inflation Estimation
1 Efficient Inflation Estimclticn by Michael F. Bryan, Stephen G. Cecchetti, and Rodney L. Wiggins II Working Paper 9707 EFFICIENT JIVFLATION ESTIMATION by Michael F. Bryan, Stephen G. Cecchetti, and Rodney L. Wiggins I1 Michael F. Bryan is assistant vice president and economist at the Federal Reserve Bank of Cleveland. Stephen G. Cecchetti is executive vice president and director of research at the Federal Reserve Bank of New York and a research associate at the National Bureau of Econonlic Research. Rodney L. Wiggins I1 is a member of the Research Department of the Federal Reserve Bank of Cleveland. The authors gratefully acknowledge the comments and assistance of Todd Clark, Ben Craig, and Scott Roger. Working papers of the Federal Reserve Bank of Cleveland are preliminary materials circulated to stimulate discussion and critical comment. The views stated herein are those of the authors and are not necessarily those of the Federal Reserve Bank of Cleveland or of the Board of Governors of the Federal Reserve System. Federal Reserve Bank of Cleveland working papers are distributed for the purpose of promoting discussion of research in progress. These papers may not have been subject to the formal editorial review accorded official Federal Reserve Bank of Cleveland publications. Working papers are now available electronically through the Cleveland Fed's home page on the World Wide Web: http:Ilwww.clev.frb.org. August 1997 Abstract This paper investigates the use of trimmed means as high-frequency estimators of inflation. The known characteristics of price change distributions, specifically the . observation that they generally exhibit high levels of kurtosis, imply that simple averages of price data are unlikely to produce efficient estimates of inflation. -
A Globally Optimal Solution to 3D ICP Point-Set Registration
1 Go-ICP: A Globally Optimal Solution to 3D ICP Point-Set Registration Jiaolong Yang, Hongdong Li, Dylan Campbell, and Yunde Jia Abstract—The Iterative Closest Point (ICP) algorithm is one of the most widely used methods for point-set registration. However, being based on local iterative optimization, ICP is known to be susceptible to local minima. Its performance critically relies on the quality of the initialization and only local optimality is guaranteed. This paper presents the first globally optimal algorithm, named Go-ICP, for Euclidean (rigid) registration of two 3D point-sets under the L2 error metric defined in ICP. The Go-ICP method is based on a branch-and-bound (BnB) scheme that searches the entire 3D motion space SE(3). By exploiting the special structure of SE(3) geometry, we derive novel upper and lower bounds for the registration error function. Local ICP is integrated into the BnB scheme, which speeds up the new method while guaranteeing global optimality. We also discuss extensions, addressing the issue of outlier robustness. The evaluation demonstrates that the proposed method is able to produce reliable registration results regardless of the initialization. Go-ICP can be applied in scenarios where an optimal solution is desirable or where a good initialization is not always available. Index Terms—3D point-set registration, global optimization, branch-and-bound, SE(3) space search, iterative closest point F 1 INTRODUCTION given an initial transformation (rotation and transla- Point-set registration is a fundamental problem in tion), it alternates between building closest-point cor- computer and robot vision. -
Trimmed Mean Group Estimation
Trimmed Mean Group Estimation Yoonseok Lee∗ Donggyu Sul† Syracuse University University of Texas at Dallas May 2020 Abstract This paper develops robust panel estimation in the form of trimmed mean group estimation. It trims outlying individuals of which the sample variances of regressors are either extremely small or large. As long as the regressors are statistically independent of the regression coefficients, the limiting distribution of the trimmed estimator can be obtained in a similar way to the standard mean group estimator. We consider two trimming methods. The first one is based on the order statistic of the sample variance of each regressor. The second one is based on the Mahalanobis depth of the sample variances of regressors. We apply them to the mean group estimation of the two-way fixed effect model with potentially heterogeneous slope parameters and to the commonly correlated regression, and we derive limiting distribution of each estimator. As an empirical illustration, we consider the effect of police on property crime rates using the U.S. state-level panel data. Keywords: Trimmed mean group estimator, Robust estimator, Heterogeneous panel, Two- way fixed effect, Commonly correlated estimator JEL Classifications: C23, C33 ∗Address: Department of Economics and Center for Policy Research, Syracuse University, 426 Eggers Hall, Syracuse, NY 13244. E-mail: [email protected] †Address: Department of Economics, University of Texas at Dallas, 800 W. Campbell Road, Richardson, TX 75080. E-mail: [email protected] 1Introduction Though it is popular in practice to impose homogeneity in panel data regression, the homogeneity restriction is often rejected (e.g., Baltagi, Bresson, and Pirotte, 2008). -
Programs Tramo (Time Series Regression with Arima Noise
PROGRAMS TRAMO (TIME SERIES REGRESSION WITH ARIMA NOISE, MISSING OBSERVATIONS AND OUTLIERS) AND SEATS (SIGNALEXTRACTIONS AREMA TIME SERIES) INSTRUCTIONS FOR THE USER (BETA VERSION: JUNIO 1997) Víctor Gómez * Agustín Maravall ** SGAPE-97001 Julio 1997 * Ministerio de Economía y Hacienda ** Banco de España This software is available at the following Internet address: http://www.bde.es The Working Papers of the Dirección General de Análisis y Programación Presupuestaria are not official statements of the Ministerio de Economía y Hacienda Abstract Brief summaries and user Instructions are presented for the programs-* TRAMO ("Time 'Series Regression with ARTMA Noise, Missing Observations and Outliers") and SEATS ("Signal Extraction in ARIMA Time Series"). TRAMO is a program for estimation and forecasting of regression models with possibly nonstationary (ARIMA) errors and any sequence of missing val- ues. The program interpolates these values, identifies and corrects for several types of outliers, and estimates special effects such as Trading Day and Easter and, in general, intervention variable type of effects. Fully automatic model identification and outlier correction procedures are available. SEATS is a program for estimation of unobserved components in time se- ries following the so-called ARlMA-model-based method. The trend, seasonal, irregular, and cyclical components are estimated and forecasted with signal extraction techniques applied to ARIMA models. The standard errors of the es- timates and forecasts are obtained and the model-based structure is exploited to answer questions of interest in short-term analysis of the data. The two programs are structured so as to be used together both for in- depth analysis of a few series (as presently done at the Bank of Spain) or for automatic routine applications to a large number of series (as presently done at Eurostat). -
Learning Entangled Single-Sample Distributions Via Iterative Trimming
Learning Entangled Single-Sample Distributions via Iterative Trimming Hui Yuan Yingyu Liang Department of Statistics and Finance Department of Computer Sciences University of Science and Technology of China University of Wisconsin-Madison Abstract have n data points from n different distributions with a common mean but different variances (the mean and all the variances are unknown), and our goal is to es- In the setting of entangled single-sample dis- timate the mean. tributions, the goal is to estimate some com- mon parameter shared by a family of distri- This setting is motivated for both theoretical and prac- butions, given one single sample from each tical reasons. From the theoretical perspective, it goes distribution. We study mean estimation and beyond the typical i.i.d. setting and raises many inter- linear regression under general conditions, esting open questions, even on basic topics like mean and analyze a simple and computationally ef- estimation for Gaussians. It can also be viewed as ficient method based on iteratively trimming a generalization of the traditional mixture modeling, samples and re-estimating the parameter on since the number of distinct mixture components can the trimmed sample set. We show that the grow with the number of samples. From the practical method in logarithmic iterations outputs an perspective, many modern applications have various estimation whose error only depends on the forms of heterogeneity, for which the i.i.d. assumption noise level of the dαne-th noisiest data point can lead to bad modeling of their data. The entan- where α is a constant and n is the sample gled single-sample setting provides potentially better size. -
1.2 Examples of Descriptive Statistics in Statistics, a Summary Statistic Is a Single Numerical Measure of an At- Tribute of a Sample
Arkansas Tech University MATH 3513: Applied Statistics I Dr. Marcel B. Finan 1.2 Examples of Descriptive Statistics In Statistics, a summary statistic is a single numerical measure of an at- tribute of a sample. In this section we discuss commonly used summary statistics. Sample Mean The sample mean (also known as average or arithmetic mean) is a mea- sure of the \center" of the data. Suppose that a sample consists of the data values x1; x2; ··· ; xn: Then the sample mean is given by the formula n x1 + x2 + ··· + xn 1 X x = = x : n n i i=1 Example 1.2.1 A random sample of 10 ATU students reported missing school 7, 6, 8, 4, 2, 7, 6, 7, 6, 5 days. Find the sample mean. Solution. The sample mean is 7+6+8+4+2+7+6+7+6+5 x = = 5:8 days 10 The Sample Variance and Standard Deviation The Standard deviation is a measure of the spread of the data around the sample mean. When the spread of data is large we expect most of the sample values to be far from the mean and when the spread is small we expect most of the sample values to be close to the mean. Suppose that a sample consists of the data values x1; x2; ··· ; xn: One way of measuring how these values are spread around the mean is to compute the deviations of these values from the mean, i.e., x1 − x; x2 − x; ··· ; xn − x and then take their average, i.e., find how far, on average, is each data value from the mean. -
Chapter 4 Truncated Distributions
Chapter 4 Truncated Distributions This chapter presents a simulation study of several of the confidence intervals first presented in Chapter 2. Theorem 2.2 on p. 50 shows that the (α, β) trimmed mean Tn is estimating a parameter μT with an asymptotic variance 2 2 equal to σW /(β−α) . The first five sections of this chapter provide the theory needed to compare the different confidence intervals. Many of these results will also be useful for comparing multiple linear regression estimators. Mixture distributions are often used as outlier models. The following two definitions and proposition are useful for finding the mean and variance of a mixture distribution. Parts a) and b) of Proposition 4.1 below show that the definition of expectation given in Definition 4.2 is the same as the usual definition for expectation if Y is a discrete or continuous random variable. Definition 4.1. The distribution of a random variable Y is a mixture distribution if the cdf of Y has the form k FY (y)= αiFWi (y) (4.1) i=1 <α < , k α ,k≥ , F y where 0 i 1 i=1 i =1 2 and Wi ( ) is the cdf of a continuous or discrete random variable Wi, i =1, ..., k. Definition 4.2. Let Y be a random variable with cdf F (y). Let h be a function such that the expected value Eh(Y )=E[h(Y )] exists. Then ∞ E[h(Y )] = h(y)dF (y). (4.2) −∞ 104 Proposition 4.1. a) If Y is a discrete random variable that has a pmf f(y) with support Y,then ∞ Eh(Y )= h(y)dF (y)= h(y)f(y). -
Statistical Inference of Truncated Normal Distribution Based on the Generalized Progressive Hybrid Censoring
entropy Article Statistical Inference of Truncated Normal Distribution Based on the Generalized Progressive Hybrid Censoring Xinyi Zeng and Wenhao Gui * Department of Mathematics, Beijing Jiaotong University, Beijing 100044, China; [email protected] * Correspondence: [email protected] Abstract: In this paper, the parameter estimation problem of a truncated normal distribution is dis- cussed based on the generalized progressive hybrid censored data. The desired maximum likelihood estimates of unknown quantities are firstly derived through the Newton–Raphson algorithm and the expectation maximization algorithm. Based on the asymptotic normality of the maximum likelihood estimators, we develop the asymptotic confidence intervals. The percentile bootstrap method is also employed in the case of the small sample size. Further, the Bayes estimates are evaluated under various loss functions like squared error, general entropy, and linex loss functions. Tierney and Kadane approximation, as well as the importance sampling approach, is applied to obtain the Bayesian estimates under proper prior distributions. The associated Bayesian credible intervals are constructed in the meantime. Extensive numerical simulations are implemented to compare the performance of different estimation methods. Finally, an authentic example is analyzed to illustrate the inference approaches. Keywords: truncated normal distribution; generalized progressive hybrid censoring scheme; ex- pectation maximization algorithm; Bayesian estimate; Tierney and Kadane approximation; impor- tance sampling Citation: Zeng, X.; Gui, W. Statistical Inference of Truncated Normal 1. Introduction Distribution Based on the 1.1. Truncated Normal Distribution Generalized Progressive Hybrid Normal distribution has played a crucial role in a diversity of research fields like Entropy 2021 23 Censoring. , , 186. reliability analysis and economics, as well as many other scientific developments. -
Power-Shrinkage and Trimming: Two Ways to Mitigate Excessive Weights
ASA Section on Survey Research Methods Power-Shrinkage and Trimming: Two Ways to Mitigate Excessive Weights Chihnan Chen1, Nanhua Duan2, Xiao-Li Meng3, Margarita Alegria4 Boston University1 UCLA2 Harvard University3 Cambridge Health Alliance and Harvard Medical School4 Abstract 2 Power-Shrinkage and Trimming Large-scale surveys often produce raw weights with very The survey weights are constructed as the reciprocals of large variations. A standard approach is to perform the selection probabilities. Let w be the survey weight, n some form of trimming, as a way to reduce potentially be the sample size, and y be the variable of interest. Let large variances in various survey estimators. The amount p 2 [0; 1] be the power-shrinkage parameter and T 2 [0; 1] of trimming is usually determined by considerations of be the trimming threshold de¯ned in terms of percentile. bias-variance trade-o®. While bias-variance trade-o® is The original weighted estimator, power-shrinkage estima- a sound principle, the trimming method itself is popular tor and the trimmed estimator for the expectation, E(y), only because of its simplicity, not because it has statisti- are given by cally desirable properties. In this paper we investigate a Xn more principled method by shrinking the variability of the wiyi log of the weights. We use the same mean squared error ² Original weighted estimator :y ¹ = i=1 (MSE) criterion to determine the amount of shrinking. w Xn The shrinking on the log scale implies shrinking weights wi by a power parameter p between [0,1]. Our investiga- i=1 tion suggests an empirical way to predict the optimal Xn choice. -
26 Sep 2019 the Trimmed Mean in Non-Parametric Regression
The Trimmed Mean in Non-parametric Regression Function Estimation Subhra Sankar Dhar1, Prashant Jha2, and Prabrisha Rakhshit3 1Department of Mathematics and Statistics, IIT Kanpur, India 2Department of Mathematics and Statistics, IIT Kanpur, India 3Department of Statistics, Rutgers University, USA September 27, 2019 Abstract This article studies a trimmed version of the Nadaraya-Watson estimator to estimate the unknown non-parametric regression function. The characterization of the estimator through minimization problem is established, and its pointwise asymptotic distribution is derived. The robustness property of the proposed estimator is also studied through breakdown point. Moreover, as the trimmed mean in the location model, here also for a wide range of trim- ming proportion, the proposed estimator poses good efficiency and high breakdown point for various cases, which is out of the ordinary property for any estimator. Furthermore, the usefulness of the proposed estimator is shown for three benchmark real data and various simulated data. Keywords: Heavy-tailed distribution; Kernel density estimator; L-estimator; Nadaraya-Watson estimator; Robust estimator. arXiv:1909.10734v2 [math.ST] 26 Sep 2019 1 Introduction We have a random sample (X1,Y1),..., (Xn,Yn), which are i.i.d. copies of (X,Y ), and the regres- sion of Y on X is defined as Yi = g(Xi)+ ei, i =1, . , n, (1) where g( ) is unknown, and e are independent copies of error random variable e with E(e X)=0 · i | and var(e X = x)= σ2(x) < for all x. Note that the condition E(e X) = 0 is essentially the | ∞ | 2000 AMS Subject Classification: 62G08 1 identifiable condition for mean regression, and it varies over the different procedures of regression. -
Appendix Methods for Calculating the Mean of a Lognormal Distribution
Appendix Methods for Calculating the Mean of a Lognormal Distribution All of the following methods are based upon a dis- NOTE: The statistical mean is theoretically the pre- tribution in which n = 9, and the constituent values ferred expression. However, I recommend against its are: general use, for the following reasons: NATURAL LOGS 1. For field-size distributions, we do not ordinarily P900/0 = 0.6MM KEY VALUES (In) anticipate that a given trend or basin will have an P800/0 = 1.25MM P99% = O.lMM infinite number of fields; instead we generally P700/0 = 2.1MM P900/0 = 0.6MM find tens to perhaps a few hundreds of fields. P600/0 = 3.3MM P84% = 0.95MM Accordingly, the continuous distribution, when n P500/0 = 5.0MM P500/0 = 5.0MM 1.6094 = 00, seems inappropriate, and leads to an inflated P400/0 = 7.6MM PI6% = 26.5MM mean field size. P300/0 = 12.1MM P100/0 = 43.0MM 2. For prospect parameters such as Area, Average P200/0 = 20.0MM PI% = 250.0MM Net Pay, Gross Rock Volume, HC-recovery in P100/0 = 43.0MM barrels per acre-foot (bbl/af) or thousand cubic Total = 95.95MM feet per acre-foot (mcf/ af), Prospect Reserves, Initial Production Rates, etc., the recommended 1. Arithmetic Mean (x) probabilistic estimating connotations of P99%, x = 95.95 + 9 = 10.66 MM P900/0, P500/0, P100/0, and PI% treat those very 2. Swanson's Mean (Msw): 0 0 0 large outcomes greater than PI% as practically 0.3 (P90 /0) + 0.4 (P50 /0) + 0.3 (P10 /0) and geologically impossible. -
Measuring Core Inflation in India: an Application of Asymmetric-Trimmed Mean Approach
Measuring Core Inflation in India: An Application of Asymmetric-Trimmed Mean Approach Motilal Bicchal* Naresh Kumar Sharma ** Abstract The paper seek to obtain an optimal asymmetric trimmed means-core inflation measure in the class of trimmed means measures when the distribution of price changes is leptokurtic and skewed to the right for any given period. Several estimators based on asymmetric trimmed mean approach are constructed, and evaluated by the conditions set out in Marques et al. (2000). The data used in the study is the monthly 69 individual price indices which are constituent components of Wholesale Price Index (WPI) and covers the period, April 1994 to April 2009, with 1993-94 as the base year. Results of the study indicate that an optimally trimmed estimator is found when we trim 29.5 per cent from the left-hand tail and 20.5 per cent from the right-hand tail of the distribution of price changes. Key words: Core inflation, underlying inflation, asymmetric trimmed mean JEL Classification: C43, E31, E52 * Ph.D Research scholar and ** Professor at Department of Economics, University of Hyderabad, P.O Central University, Hyderabad-500046, India. Corresponding author: [email protected] We are very grateful to Carlos Robalo Marques of Banco de Portugal for providing the program for computing the asymmetric trimmed-mean as well as for his valuable comments and suggestions. We also wish to thank Bandi Kamaiah for going through previous draft and for offering comments, which improved the presentation of this paper. 1. Introduction Various approaches to measuring core inflation have been discussed in the literature.