hcfcirc.boot: Bootstrap 2-sample mean test for circular data

View source: R/hcfcirc.boot.R

Bootstrap 2-sample mean test for circular dataR Documentation

Bootstrap 2-sample mean test for circular data

Description

Bootstrap 2-sample mean test for circular data.

Usage

hcfcirc.boot(u1, u2, rads = TRUE, B = 999)
lrcirc.boot(u1, u2, rads = TRUE, B = 999)
hclrcirc.boot(u1, u2, rads = TRUE, B = 999)
embedcirc.boot(u1, u2, rads = TRUE, B = 999)
hetcirc.boot(u1, u2, rads = TRUE, B = 999)

Arguments

u1

A numeric vector containing the data of the first sample.

u2

A numeric vector containing the data of the first sample.

rads

If the data are in radians, this should be TRUE and FALSE otherwise.

B

The number of bootstraps to perform.

Details

The high concentration (hcfcirc.boot), the log-likelihood ratio test (lrcirc.boot), high concentration log-likelihood ratio (hclrcirc.boot), embedding approach (embedcirc.boot), or the non equal concentration parameters approach (hetcirc.boot) is used.

Value

This is an "htest"class object. Thus it returns a list including:

statistic

The test statistic value.

parameter

The degrees of freedom of the test. Since these are bootstrap based tests this is "NA".

p.value

The p-value of the test.

alternative

A character with the alternative hypothesis.

method

A character with the test used.

data.name

A character vector with two elements.

Author(s)

Michail Tsagris.

R implementation and documentation: Michail Tsagris mtsagris@uoc.gr.

References

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. (2022). An investigation of hypothesis testing procedures for circular and spherical mean vectors. Communications in Statistics-Simulation and Computation (Accepted for publication).

See Also

hcf.circaov, hcfcircboot, het.aov

Examples

u1 <- rvonmises(20, 2.4, 5)
u2 <- rvonmises(20, 2.4, 10)
hcfcirc.boot(u1, u2)

Directional documentation built on Oct. 12, 2023, 1:07 a.m.