Clean Kinematic Samples in Dwarf Spheroidals: An Algorithm for Evaluating Membership and Estimating Distribution Parameters When Contamination is Present
arXiv:0811.1990 · doi:10.1088/0004-6256/137/2/3109
Abstract
(abridged) We develop an algorithm for estimating parameters of a distribution sampled with contamination, employing a statistical technique known as ``expectation maximization'' (EM). Given models for both member and contaminant populations, the EM algorithm iteratively evaluates the membership probability of each discrete data point, then uses those probabilities to update parameter estimates for member and contaminant distributions. The EM approach has wide applicability to the analysis of astronomical data. Here we tailor an EM algorithm to operate on spectroscopic samples obtained with the Michigan-MIKE Fiber System (MMFS) as part of our Magellan survey of stellar radial velocities in nearby dwarf spheroidal (dSph) galaxies. These samples are presented in a companion paper and contain discrete measurements of line-of-sight velocity, projected position, and Mg index for ~1000 - 2500 stars per dSph, including some fraction of contamination by foreground Milky Way stars. The EM algorithm quantifies both dSph and contaminant distributions, returning maximum-likelihood estimates of the means and variances, as well as the probability that each star is a dSph member. Applied to our MMFS data, the EM algorithm identifies more than 5000 probable dSph members. We test the performance of the EM algorithm on simulated data sets that represent a range of sample size, level of contamination, and amount of overlap between dSph and contaminant velocity distributions. The simulations establish that for samples ranging from large (N ~3000) to small (N~30), the EM algorithm distinguishes members from contaminants and returns accurate parameter estimates much more reliably than conventional methods of contaminant removal (e.g., sigma clipping).
Accepted for publication in The Astronomical Journal. Download pdf with full-resolution figures from http://www.ast.cam.ac.uk/~walker/dsph_em.pdf
References in corpus (8)
- The Kinematics of the Ultra-Faint Milky Way Satellites: Solving the Missing Satellite Problem
- Cats and Dogs, Hair and A Hero: A Quintet of New Milky Way Companions
- The accretion origin of the Milky Way's stellar halo
- Stellar Velocities in the Carina, Fornax, Sculptor and Sextans dSph Galaxies: Data from the Magellan/MMFS Survey
- The Velocity Dispersion Profile of the Remote Dwarf Spheroidal Galaxy Leo I: A Tidal Hit and Run?
- Stellar kinematics and metallicities in the Leo I dwarf spheroidal galaxy -- wide field implications for galactic evolution
- Mass modelling of dwarf spheroidal galaxies: the effect of unbound stars from tidal tails and the Milky Way
- The Michigan/MIKE Fiber System Survey of Stellar Radial Velocities in Dwarf Spheroidal Galaxies: Acquisition and Reduction of Data