  We consider the estimation of the global mode of a density under some decay
rate condition around the global mode. We show that the maximum of a histogram,
with proper choice of bandwidth, achieves the minimax rate that we establish
for the setting that we consider. This is based on knowledge of the decay rate.
Addressing the situation where the decay rate is unknown, we propose a
multiscale variant consisting in the recursive refinement of a histogram, which
is shown to be minimax adaptive. These methods run in linear time, and we prove
in an appendix that this is best possible: There is no estimation procedure
that runs in sublinear time that achieves the minimax rate.
