dVm {sdetorus} | R Documentation |
Computes the density of a von Mises in a numerically stable way.
dVm(x, mu, kappa)
x |
evaluation angles, not necessary in |
mu |
circular mean. |
kappa |
non-negative concentration parameter. |
A vector of the same length as x
containing the density.
Jammalamadaka, S. R. and SenGupta, A. (2001) Topics in Circular Statistics. World Scientific, Singapore. doi:10.1142/4031
x <- seq(-pi, pi, l = 200)
plot(x, x, type = "n", ylab = "Density", ylim = c(0, 1))
for (i in 0:20) {
lines(x, dVm(x = x, mu = 0, kappa = 5 * i / 20),
col = rainbow(21)[i + 1])
}