Bootstrap ANOVA for (hyper-)spherical data.
hcfboot(x, ina, B = 999)
hetboot(x, ina, B = 999)This is an "htest"class object. Thus it returns a list including:
The test statistic value.
The degrees of freedom of the test. Since these are bootstrap based tests this is "NA".
The p-value of the test.
A character with the alternative hypothesis.
A character with the test used.
A character vector with two elements.
A matrix with the combined data (from all groups) in Euclidean coordinates, i.e. unit vectors.
The grouping variables. A factor or a numerical vector specifying the groups to which each observation belongs to.
The number of bootstraps to perform.
Michail Tsagris.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr.
The high concentration (hcfboot), or the non equal concentration parameters approach (hetboot) is used.
Mardia K. V. and Jupp P. E. (2000). Directional statistics. Chicester: John Wiley & Sons.
Rumcheva P. and Presnell B. (2017). An improved test of equality of mean directions for the Langevin-von Mises-Fisher distribution. Australian & New Zealand Journal of Statistics, 59(1): 119--135.
Tsagris M. and Alenazi A. (2024). An investigation of hypothesis testing procedures for circular and spherical mean vectors. Communications in Statistics-Simulation and Computation, 53(3): 1387--1408.
hcf.boot, hcf.aov
x <- rvmf(60, rnorm(3), 10)
ina <- rep(1:3, each = 20)
hcfboot(x, ina)
Run the code above in your browser using DataLab