Plots the empirical distribution of a data set against the best fitting von Mises distribution function.
Usage
pp.plot(x, ref.line=TRUE)
Arguments
x
vector of angular measurements in radians.
ref.line
logical flag: if TRUE a 45 degree reference line is added to the plot. Default is TRUE.
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
NULL
Note
A plot is created on the current graphics device.
References
Jammalamadaka, S. Rao and SenGupta, A. (2001). Topics in Circular Statistics, Section 10.2, World Scientific Press, Singapore.
See Also
vm.ml
Examples
data <- rvm(25, 0, 3)
pp.plot(data)
data <- c(rvm(20, 0, 7), rvm(20, pi, 7))
pp.plot(data)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(CircStats)
Loading required package: MASS
Loading required package: boot
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/CircStats/pp.plot.Rd_%03d_medium.png", width=480, height=480)
> ### Name: pp.plot
> ### Title: von Mises Probability-Probability Plot
> ### Aliases: pp.plot
> ### Keywords: hplot
>
> ### ** Examples
>
> data <- rvm(25, 0, 3)
> pp.plot(data)
mu kappa
1 0.1231361 2.441935
> data <- c(rvm(20, 0, 7), rvm(20, pi, 7))
> pp.plot(data)
mu kappa
1 6.140872 0.03234818
>
>
>
>
>
> dev.off()
null device
1
>