wrappedCauchy {circular}

Wrapped Cauchy Density Function

Description

Density, and random generation for the wrapped Cauchy circular distribution.

Usage

dwrappedcauchy(x, mu = circular(0), rho = exp(-1))
rwrappedcauchy(n, mu = circular(0), rho = exp(-1), control.circular=list())

Arguments

x

a vector. The object is coerced to class circular.

n

number of observations.

mu

mean direction of the distribution as a circular object.

rho

concentration parameter of the distribution. rho must be in the interval from 0 to 1.

control.circular

the attribute of the resulting object.

Value

dwrappedcauchy gives the density and rwrappedcauchy generates random deviates.

Author(s)

Claudio Agostinelli and Ulric Lund

References

Jammalamadaka, S. Rao and SenGupta, A. (2001). Topics in Circular Statistics, Section 2.2.7, World Scientific Press, Singapore.

Examples

data1 <- rwrappedcauchy(100, mu=circular(0), rho=0.7, 
  control.circular=list(units="degrees"))
plot(data1)

ff <- function(x) dwrappedcauchy(x, mu=circular(pi), rho=0.7)
curve.circular(ff, join=TRUE, xlim=c(-2, 1), 
  main="Density of a Wrapped Cauchy Distribution \n mu=pi, rho=0.7")