circular (version 0.3-7)

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

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.

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 
 mu=pi, rho=0.7")