# Shape Restricted Estimation- In the Search for Dark Matter

Dr. Michael Woodroofe

L. J. Savage Professor of Mathematics and Statistics

Department of Statistics

University of Michigan

In addition to visible matter, like stars, there is matter that
astronomers cannot see. They know that it is there from its
gravitational effects on the positions and velocity of stars. It's
nature and even location are not well known at present, but there is a
wide belief that the vast majority of matter in the Universe is dark.
The talk centers on statistical issues that arise in the non-parametric
estimation of dark matter in the dwarf spheroidal galaxies that are
satellites of the Milky Way. Statistically, this is an inverse problem
with missing data. The distribution of (total) mass together with the
law of gravity determines the positions and velocities of stars, and the
statistical problem is to recover the distribution of mass from
observed positions and velocities. Complete positions and velocities
cannot be observed here--only the velocities in the line of sight (from
Earth) and the projections of positions on the plane orthogonal to the
line of sight. The distribution of projected position can be estimated
quite precisely from large ($\sim 150,000$ stars) photometric studies.
Velocities require some time to measure, however, and previous studies
have been too small ($\sim 30$ stars) to allow estimation in any
generality. Instead, astronomers have been forced to adopt highly
structured parametric models. New instrumentation making much larger
samples ($\sim 1500$ stars) available and non-parametric methods
feasible. Assuming spherical symmetry and isotropy, let $M(r)$ denote
the total mass within $r$ of the center of a dwarf spheroidal galaxy.
Then Jeans' Equation (from classical potential theory) shows that $M(r) =
c\Psi''(r)/f(r)$, where $c$ is a known constant, $f$ is a (nearly)
known non-negative function, and $\Psi$ is a function for which an
unbiased, $\sqrt{n}$-consistent estimator exists. From Jeans' Equation,
$\Psi$ is a convex function, and $\Psi''(r)/f(r)$ is non-decreasing in
$r$. The unbiased estimator, $\Psi^{\#}$ say, does not satisfy these
shape restrictions and can be improved by imposing them. The talk
centers on methods for accomplishing this, and properties of the
resulting estimators. Some preliminary work that relaxes the isotropy
assumption will also be discussed. This is joint work with Jayanta Pal
and Xiao Wang (UM Statistics) and Mario Mateo and Matt Walker (UM
Astronomy).