Skip to contents

Circular mean, standard deviation, variance, quasi-quantiles, 95% confidence angle, standardized skewness and kurtosis

Usage

circular_summary(x, w = NULL, axial = TRUE, na.rm = FALSE)

Arguments

x

numeric vector. Values in degrees.

w

(optional) Weights. A vector of positive numbers and of the same length as x.

axial

logical. Whether the data are axial, i.e. pi-periodical (TRUE, the default) or directional, i.e. \(2 \pi\)-periodical (FALSE).

na.rm

logical value indicating whether NA values in x should be stripped before the computation proceeds.

Value

named vector

Examples

data("nuvel1")
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")
sa.por <- PoR_shmax(san_andreas, PoR, "right")
circular_summary(sa.por$azi.PoR)
#>            n         mean           sd          var          25% quasi-median 
#>  407.0000000  139.3655975   20.0226130    0.2167045  126.5036202  136.1189938 
#>          75%       median        95%CI     skewness     kurtosis            R 
#>  147.8526166  137.5578600    3.5838559   -0.4009978    2.4250500    0.7832955 
circular_summary(sa.por$azi.PoR, w = 1 / san_andreas$unc)
#>            n         mean           sd          var          25% quasi-median 
#>  407.0000000  142.4029884   20.3016715    0.2220562  126.5036202  136.1189938 
#>          75%       median        95%CI     skewness     kurtosis            R 
#>  147.8526166  137.5578600    3.7636758   -0.4526369    1.4620392    0.7832955