  The recent positron excess in cosmic rays (CR) observed by the PAMELA
satellite may be a signal for dark matter (DM) annihilation. When these
measurements are combined with those from FERMI on the total ($e^++e^-$) flux
and from PAMELA itself on the $\bar p/p$ ratio, these and other results are
difficult to reconcile with traditional models of DM, including the
conventional mSUGRA version of Supersymmetry even if boosts as large as
$10^{3-4}$ are allowed. In this paper, we combine the results of a previously
obtained scan over a more general 19-parameter subspace of the MSSM with a
corresponding scan over astrophysical parameters that describe the propagation
of CR. We then ascertain whether or not a good fit to this CR data can be
obtained with relatively small boost factors while simultaneously satisfying
the additional constraints arising from gamma ray data. We find that a specific
subclass of MSSM models where the LSP is mostly pure bino and annihilates
almost exclusively into $\tau$ pairs comes very close to satisfying these
requirements. The lightest $\tilde \tau$ in this set of models is found to be
relatively close in mass to the LSP and is in some cases the nLSP. These models
lead to a significant improvement in the overall fit to the data by an amount
$\Delta \chi^2 \sim 1/$dof in comparison to the best fit without Supersymmetry
while employing boosts $\sim 100$. The implications of these models for future
experiments are discussed.
