It produces a few summary measures for circular data.
Usage
circ.summary(u, rads = FALSE, plot = TRUE)
Arguments
u
The circular data.
rads
If the data are in rads, then this should be TRUE, otherwise FALSE.
plot
If you want to see the data plotted on a cicrle make this TRUE.
Details
It returns the circular mean, mean resultant length, variance, standard deviation and concentration parameter. So, basically it returns the estimated values of the parameters of the von Mises distribution.
Value
A list including:
mesos
The circular mean direction.
confint
The 95% confidence interval for the circular mean direction.
kappa
The concentration parameter.
MRL
The mean resultant length.
circvariance
The circular variance.
circstd
The circular standard deviation.
Author(s)
Michail Tsagris
R implementation and documentation: Michail Tsagris <mtsagris@yahoo.gr> and Giorgos Athineou <athineou@csd.uoc.gr>
References
Mardia, K. V. and Jupp, P. E. (2000). Directional statistics. Chicester: John Wiley & Sons.
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(Directional)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Directional/circ.summary.Rd_%03d_medium.png", width=480, height=480)
> ### Name: circ.summary
> ### Title: Summary statistics for circular data.
> ### Aliases: circ.summary
> ### Keywords: Summary statistics Von Mises distribution
>
> ### ** Examples
>
> x <- rvonmises(50, 2.5, 15, rads = TRUE)
> circ.summary(x, rads = TRUE, plot = TRUE)
$mesos
[1] 2.486222
$confint
[1] 2.483557 2.488887
$kappa
[1] 12.08629
$MRL
[1] 0.9576931
$circvariance
[1] 0.0423069
$circstd
[1] 0.2940337
>
>
>
>
>
> dev.off()
null device
1
>