Automatic bandwidth selection for circular density estimation

Given angular data θ1,…,θn∈[0,2π)θ1,…,θn∈[0,2π) a common objective is to estimate the density. In case that a kernel estimator is used, bandwidth selection is crucial to the performance. A “plug-in rule” for the bandwidth, which is based on the concentration of a reference density, namely, the von Mises distribution is obtained. It is seen that this is equivalent to the usual Euclidean plug-in rule in the case where the concentration becomes large. In case that the concentration parameter is unknown, alternative methods are explored which are intended to be robust to departures from the reference density. Simulations indicate that “wrapped estimators” can perform well in this context. The methods are applied to a real bivariate dataset concerning protein structure.