Last data update: 2014.03.03

R: Circular Summary Statistics
circ.summaryR Documentation

Circular Summary Statistics

Description

Computes circular summary statistics including the sample size, mean direction and mean resultant length.

Usage

circ.summary(x)

Arguments

x

vector of data points measured in radians.

Details

Each observation is treated as a unit vector or a point on the unit circle. The resultant vector of the observations is found, and the direction of the resultant vector is returned as well as its length divided by the sample size.

Value

Returns a data frame with variables n, the sample size; mean.dir, the sample mean direction; and rho, the sample mean resultant length.

References

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

See Also

circ.mean, circ.disp, est.kappa, est.rho.

Examples

# Compute summary statistics of a random sample of observations. 
data <- runif(50, 0, pi)
circ.summary(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/circ.summary.Rd_%03d_medium.png", width=480, height=480)
> ### Name: circ.summary
> ### Title: Circular Summary Statistics
> ### Aliases: circ.summary
> ### Keywords: misc
> 
> ### ** Examples
> 
> # Compute summary statistics of a random sample of observations. 
> data <- runif(50, 0, pi)
> circ.summary(data)
   n mean.dir       rho
1 50 1.783193 0.6752342
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>