pp.plot: von Mises Probability-Probability Plot

In circular: Circular Statistics

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

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.

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.

Value

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

Author(s)

Claudio Agostinelli and Ulric Lund

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)

Example output

Attaching package: 'circular'

The following objects are masked from 'package:stats':

    sd, var

circular documentation built on May 2, 2019, 4:42 p.m.