circular (version 0.5-1)

pp.plot: von Mises Probability-Probability Plot

Description

Plots the empirical distribution of a data set against the best fitting von Mises distribution function.

Usage

pp.plot(x, ref.line = TRUE, tol=1e-20,  xlab = "von Mises Distribution", 
  ylab = "Empirical Distribution", control.circular = list(), ...)

Value

a list with the estimated mean and concentration parameter for a von Mises distribution.

Arguments

x

a vector. The object is coerced to class circular.

ref.line

logical, if TRUE a 45 degree reference line is added to the plot. Default is TRUE.

tol

parameter passed to pvonmises.

xlab, ylab

labels of the axis.

control.circular

the attribute of the resulting object.

...

parameters passed to the plot.default function.

Author

Claudio Agostinelli and Ulric Lund

Details

The maximum likelihood estimates of the parameters of the von Mises distribution are computed from the given data set. The empirical distribution function is plotted against a von Mises distribution function with parameters given by the MLEs computed.

References

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

See Also

mle.vonmises

Examples

x <- rvonmises(n=25, mu=circular(0), kappa=3)
pp.plot(x)
x <- c(rvonmises(n=20, mu=circular(0), kappa=7), 
  rvonmises(n=20, mu=circular(pi), kappa=7))
pp.plot(x)