Confidence Bands for the Gradient Lines of a Density Function

We consider the problem of estimating the gradient ascent line of a density originating at a given point. Going beyond mere consistency, we establish a weak convergence result for a plugin estimator based on a kernel density estimator of the density. We then leverage that result to construct a confidence region for the gradient ascent line, including by bootstrap.