vonMises {circular} | R Documentation |
Density, distribution function and random generation for the von Mises circular distribution.
dvonmises(x, mu, kappa) pvonmises(q, mu, kappa, tol = 1e-020) rvonmises(n, mu, kappa, control.circular=list())
x, q |
a vector. The object is coerced to class
circular . |
n |
number of observations. |
mu |
mean direction of the distribution. The object is coerced to class circular . |
kappa |
non-negative numeric value for the concentration parameter of the distribution. |
control.circular |
the attribute of the resulting object. |
tol |
the precision in evaluating the distributon function. |
dvonmises
gives the density, pvonmises
gives the
distribution function and rvonmises
generates random deviates.
Since version 0.3-5 the random deviates are generated using a C code.
Claudio Agostinelli, Ulric Lund and Harry Southworth
Jammalamadaka, S. Rao and SenGupta, A. (2001). Topics in Circular Statistics, Section 2.2.4, World Scientific Press, Singapore.
data1 <- rvonmises(100, circular(0), 10, control.circular=list(units="degrees")) plot(data1) ff <- function(x) dvonmises(x, mu=circular(pi), kappa=10) curve.circular(ff, join=TRUE, xlim=c(-2.3, 1), main="Density of a VonMises Distribution \n mu=pi, kappa=10") ff <- function(x) pvonmises(x, mu=circular(pi), kappa=10) curve.circular(ff, join=FALSE, xlim=c(-2, 2), ylim=c(-2, 1), to=(2*pi-3*.Machine$double.eps), modulo="asis", nosort=TRUE, main="Probability of a VonMises Distribution \n mu=pi, kappa=10")