circular (version 0.5-1)

wrappedCauchy: 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())

Value

dwrappedcauchy gives the density and rwrappedcauchy generates random deviates.

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.

Author

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")