wrpcauchy.ml: Wrapped Cauchy Maximum Likelihood Estimates

In CircStats: Circular Statistics, from "Topics in Circular Statistics" (2001)

Description

Computes the maximum likelihood estimates of the location and scale parameters for a wrapped Cauchy distribution.

Usage

wrpcauchy.ml(x, mu0, rho0, acc=1e-015)

Arguments

x	vector of angular data measured in radians.
mu0	initial estimate of the location parameter.
rho0	initial estimate of the scale parameter. rho0 should be between in [0,1).
acc	degree of accuracy in the approximation of the estimates. The default value is 1e-15. See below for details.

Details

An iterative algorithm due to Kent and Tyler (1988) is used. Initial values of the MLE's are required. In the estimation process, estimates of quantities mu1 and mu2 (see Kent and Tyler) are updated iteratively. When both mu1 and mu2 change by less than acc from one iteration to the next, the process ends, and mu and rho are computed from the final estimates of mu1 and mu2.

Value

A dataframe is returned with the variables mu and rho, where mu and rho are the approximations of the MLE of the location and scale parameter, respectively.

References

Kent, J. and Tyler, D. (1988). Maximum likelihood estimation for the wrapped Cauchy distribution. Journal of Applied Statistics, 15, 2, 247-254.

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

See Also

dwrpcauchy, rwrpcauchy

Examples

# Generate data from a wrapped Cauchy distribution.
data <- rwrpcauchy(50, 0, 0.75)
# Compute the sample mean direction and resultant length.
mu0 <- circ.mean(data)
rho0 <- est.rho(data)
# Estimate the MLE's of the Cauchy distribution parameters.
wrpcauchy.ml(data, mu0, rho0)

Example output

Loading required package: MASS
Loading required package: boot
          mu       rho
1 0.04576996 0.7660626

CircStats documentation built on May 2, 2019, 2:24 a.m.