mixedvonmises: Mixture of von Mises Distributions

Description

Density and random generation for the mixed von Mises circular distribution.

Usage

dmixedvonmises(x, mu1, mu2, kappa1, kappa2, prop)
rmixedvonmises(n, mu1, mu2, kappa1, kappa2, prop, control.circular = list())
pmixedvonmises(q, mu1, mu2, kappa1, kappa2, prop, from=NULL, tol = 1e-020)

Value

dmixedvonmises gives the density, pmixedvonmises gives the distribution function and rmixedvonmises generates random deviates.

Arguments

x, q: a vector. The object is coerced to class circular.
n: number of observations.
mu1: mean direction of one of the two von Mises distributions as a circular object.
mu2: mean direction of the other von Mises distribution as a circular object.
kappa1: concentration parameter of one of the two von Mises distributions.
kappa2: concentration parameter of the other von Mises distribution.
prop: mixing proportion.
from: if NULL is set equal to \(0\) (Notice the difference from the corresponding vonmises distribution). This is the value from which the pmixedvonmises is evaluated. It should be a circular object.
tol: the precision in evaluating the distribution function or the quantile.
control.circular: the attribute of the resulting object.

Author

Claudio Agostinelli and Ulric Lund

Examples

Run this code


x <- rmixedvonmises(n=100, mu1=circular(0), mu2=circular(pi), kappa1=15, 
  kappa2=15, prop=0.5)
plot(x)