(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).

We present stellar velocity dispersion profiles for seven Milky Way dwarf spheroidal (dSph) satellite galaxies. We have measured 8394 line-of-sight velocities (+/- 2.5 km/s) for 6804 stars from high-resolution spectra obtained at the Magellan and MMT telescopes. We combine these new data with previously published velocities to obtain the largest available kinematic samples, which include more than 5500 dSph members. All the measured dSphs have stellar velocity dispersion of order 10 km/s that remains approximately constant with distance from the dSph center, out to and in some cases beyond the radius at which the mean surface brightness falls to the background level. Assuming dSphs reside within dark matter halos characterized by the NFW density profile, we obtain reasonable fits to the empirical velocity dispersion profiles. These fits imply that, among the seven dSphs, M_vir ~ 10^[8-9] M_sun. The mass enclosed at a radius of 600 pc, the region common to all data sets, ranges from (2-7) x 10^7 M_sun .

We introduce a stellar velocity survey of dwarf spheroidal galaxies, undertaken using the Michigan/MIKE Fiber System (MMFS) at the Magellan/Clay 6.5 m telescope at Las Campanas Observatory. As of 2006 November we have used MMFS to collect 6415 high-resolution (R= 20000-25000) spectra from 5180 stars in four dwarf spheroidal galaxies: Carina, Fornax, Sculptor and Sextans. Spectra sample the range 5140-5180 Angstroms, which includes the prominent magnesium triplet absorption feature. We measure radial velocity (RV) to a median precision of 2.0 km/s for stars as faint as V ~ 20.5. From the spectra we also are able to measure the strength of iron and magnesium absorption features using spectral indices that correlate with effective temperature, surface gravity and chemical abundance. Measurement of line strength allows us to identify interloping foreground stars independently of velocity, and to examine the metallicity distribution among dSph members. Here we present detailed descriptions of MMFS, our target selection and spectroscopic observations, the data reduction procedure, and error analysis. We compare our RV results to previously published measurements for individual stars. In some cases we find evidence for a mild, velocity-dependent offset between the RVs we measure using the magnesium triplet and previously published RV measurements derived from the infrared calcium triplet. In companion papers we will present the complete data sets and kinematic analyses of these new observations.

We present multifiber echelle radial velocity results for 551 stars in the Sextans dwarf spheroidal galaxy and identify 294 stars as probable Sextans members. The projected velocity dispersion profile of the binned data remains flat to a maximum angular radius of $30^{\prime}$. We introduce a nonparametric technique for estimating the projected velocity dispersion surface, and use this to search for kinematic substructure. Our data do not confirm previous reports of a kinematically distinct stellar population at the Sextans center. Instead we detect a region near the Sextans core radius that is kinematically colder than the overall Sextans sample with 95% confidence.

We present new radial velocity results for 176 stars in the Fornax dwarf spheroidal galaxy, of which at least 156 are probable Fornax members. We combine with previously published data to obtain a radial velocity sample with 206 stars, of which at least 176 are probable Fornax members. We detect the hint of rotation about an axis near Fornax's morphological minor axis, though the significance of the rotation signal in the galactic rest frame is sensitive to the adopted value of Fornax's proper motion. Regardless, the observed stellar kinematics are dominated by random motions, and we do not find kinematic evidence of tidal disruption. The projected velocity dispersion profile of the binned dataset remains flat over the sampled region, which reaches a maximum angular radius of 65 arcmin. Single-component King models in which mass follows light fail to reproduce the observed flatness of the velocity dispersion profile. Two-component (luminous plus dark matter) models can reproduce the data, provided the dark component extends sufficiently beyond the luminous component, and the central dark matter density is of the same order as the central luminous density. These requirements suggest a more massive, darker Fornax than standard core-fitting analyses have previously concluded, with M/L_V over the sampled region reaching 10 to 40 times the M/L_V of the luminous component. We also apply a non-parametric mass estimation technique, introduced in a companion paper. Though it is designed to operate on datasets containing velocities for $>$1000 stars, the estimation yields preliminary results suggesting M/L_V ~ 15 inside r < 1.5 kpc.

Thanks to instrumental advances, new, very large kinematic datasets for nearby dwarf spheroidal (dSph) galaxies are on the horizon. A key aim of these datasets is to help determine the distribution of dark matter in these galaxies. Past analyses have generally relied on specific dynamical models or highly restrictive dynamical assumptions. We describe a new, non-parametric analysis of the kinematics of nearby dSph galaxies designed to take full advantage of the future large datasets. The method takes as input the projected positions and radial velocities of stars known to be members of the galaxies, but does not use any parametric dynamical model, nor the assumption that the mass distribution follows that of the visible matter. The problem of estimating the radial mass distribution, M(r) (the mass interior to true radius r), is converted into a problem of estimating a regression function non-parametrically. From the Jeans Equation we show that the unknown regression function is subject to fundamental shape restrictions which we exploit in our analysis using statistical techniques borrowed from isotonic estimation and spline smoothing. Simulations indicate that M(r) can be estimated to within a factor of two or better with samples as small as 1000 stars over almost the entire radial range sampled by the kinematic data. The technique is applied to a sample of 181 stars in the Fornax dSph galaxy. We show that the galaxy contains a significant, extended dark halo some ten times more massive than its baryonic component. Though applied here to dSph kinematics, this approach can be used in the analysis of any kinematically hot stellar system in which the radial velocity field is discretely sampled.