von Mises-Fisher kernel density estimation for (hyper-)spherical data | R Documentation |
A von Mises-Fisher kernel is used for the non parametric density estimation.
vmf.kde(x, h, thumb = "none")
x |
A matrix with unit vectors, i.e. the data being expressed in Euclidean cordinates. |
h |
The bandwidth to be used. |
thumb |
If this is "none", the given bandwidth is used. If it is "rot" the rule of thumb suggested by Garcia-Portugues (2013) is used. |
A von Mises-Fisher kernel is used for the non parametric density estimation.
A list including:
h |
The bandwidth used. |
f |
A vector with the kernel density estimate calculated for each of the data points. |
Michail Tsagris.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr and Giorgos Athineou <gioathineou@gmail.com>.
Garcia Portugues, E. (2013). Exact risk improvement of bandwidth selectors for kernel density estimation with directional data. Electronic Journal of Statistics, 7, 1655–1685.
vmfkde.tune, vm.kde, vmf.mle, vmkde.tune
x <- rvmf(100, rnorm(5), 15)
h <- vmfkde.tune(x)[1]
f1 <- vmf.kde(x, h = h, thumb = "none")
f2 <- vmf.kde(x, h = h, thumb = "rot")
f1$h ; f2$h
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.