🤖 AI Summary
This work addresses the long-standing challenge of efficient and stable bandwidth selection in Beta kernel density estimation, where existing iterative optimization methods are computationally expensive and prone to instability. The authors propose the first closed-form bandwidth selector, derived by approximating the unweighted asymptotic mean integrated squared error (AMISE) via the method of moments, reducing bandwidth computation to constant time complexity. To handle integrability issues near boundaries—particularly for U-shaped and J-shaped distributions—they introduce tailored boundary-aware heuristics. In Monte Carlo simulations, the proposed method achieves accuracy comparable to numerical optimization while accelerating computation by over 35,000-fold. Applied to real-world socioeconomic data, it effectively mitigates the boundary bias commonly observed with Gaussian kernels, substantially improving both estimation efficiency and stability.
📝 Abstract
The Beta kernel estimator offers a theoretically superior alternative to the Gaussian kernel for unit interval data, eliminating boundary bias without requiring reflection or transformation. However, its adoption remains limited by the lack of a reliable bandwidth selector; practitioners currently rely on iterative optimization methods that are computationally expensive and prone to instability. We derive the ``\rot,''a fast, closed-form bandwidth selector based on the unweighted Asymptotic Mean Integrated Squared Error (AMISE) of a beta reference distribution. To address boundary integrability issues, we introduce a principled heuristic for U-shaped and J-shaped distributions. By employing a method-of-moments approximation, we reduce the bandwidth selection complexity from iterative optimization to $\mathcal{O}(1)$. Extensive Monte Carlo simulations demonstrate that our rule matches the accuracy of numerical optimization while delivering a speedup of over 35,000 times. Real-world validation on socioeconomic data shows that it avoids the ``vanishing boundary''and ``shoulder''artifacts common to Gaussian-based methods. We provide a comprehensive, open-source Python package to facilitate the immediate adoption of the Beta kernel as a drop-in replacement for standard density estimation tools.