Last chance! 50% off unlimited learning
Sale ends in
It produces a few summary measures for grouped circular data.
group.vm(group, fi, rads = FALSE)
A matrix denoting the classes. Each row consists of two numbers, the lower and upper points of each class.
The frequency of each class of data.
If the data are in rads, then this should be TRUE, otherwise FALSE.
A list including:
The circular mean direction.
The 95% confidence interval for the circular mean direction.
The concentration parameter.
The mean resultant length.
The circular variance.
The circular standard deviation.
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. The mena resultant length though is group corrected.
Pewsey Arthur, Markus Neuhauser and Graeme D. Ruxton (2013). Circular statistics in R. Oxford University Press.
Mardia K. V. and Jupp P. E. (2000). Directional statistics. Chicester: John Wiley & Sons.
# NOT RUN {
x <- rvonmises(200, 3, 10)
a <- circ.summary(x, rads = TRUE, plot = FALSE)
group <- seq(min(x) - 0.1, max(x) + 0.1, length = 6)
y <- cut(x, breaks = group, length = 6)
group <- matrix( c( group[1], rep(group[2:5], each = 2), group[6]), ncol = 2, byrow = TRUE)
fi <- as.vector( table(y) )
b <- group.vm(group, fi, rads = TRUE)
a
b
# }
Run the code above in your browser using DataLab