  The mass function of 51 Massive Compact Objects (MACHOs) detected toward the
Galactic bulge is statistically estimated from Einstein ring crossing times
$t_{\rm e}$. For a Gaussian mass function, the best fitting parameters are
$\langle \log (m/M_{\odot}) \rangle = -1.12$ and $\sigma_{\log (m/M_{\odot})} =
0.57$. If the mass function follows a power-law distribution, the best fitting
mass cut-off and power are $m_{\rm cut} = 0.04\ M_{\odot}$ and $p = -2.1$. Both
mass spectra indicate that a significant fraction of events are caused by
MACHOs in substellar mass range. Both best determined mass functions are
compared with that obtained from Hubble Space Telescope (HST) observations. The
power law is marginally favored ($2.3\sigma$) over the Gaussian mass function,
and significantly over HST mass function ($5.5\sigma$), indicating that MACHOs
have a different mass function from stars in the solar neighborhood. In
addition, the fact that all the models have very poor fits to the longest four
events with $t_{e} \geq 70\ {\rm days}$ remains a puzzle.
