Bootstrap James and Hotelling test for 2 independent sample mean vectors | R Documentation |
Bootstrap James and Hotelling test for 2 independent sample mean vectors.
boot.james(y1, y2, R = 999)
boot.hotel2(y1, y2, R = 999)
y1 |
A numerical matrix with the data of the one sample. |
y2 |
A numerical matrix with the data of the other sample. |
R |
The number of bootstrap samples to use. |
We bootstrap the 2-samples James (does not assume equal covariance matrics) and Hotelling test
(assumes equal covariance matrics). The difference is that the Hotelling test statistic assumes
equaility of the covariance matrices, which if violated leads to inlfated type I errors.
Bootstrap calibration though takes care of this issue. As for the bootstrap calibration, instead
of sampling B times from each sample, we sample sqrt{B}
from each of them and then take all pairs.
Each bootstrap sample is independent of each other, hence there is no violation of the theory
(Chatzipantsiou et al., 2019).
The bootstrap p-value.
Michail Tsagris.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr.
G.S. James (1954). Tests of Linear Hypothese in Univariate and Multivariate Analysis when the Ratios of the Population Variances are Unknown. Biometrika, 41(1/2): 19-43
Efron Bradley and Robert J. Tibshirani (1993). An introduction to the bootstrap. New York: Chapman & Hall/CRC.
Chatzipantsiou C., Dimitriadis M., Papadakis M. and Tsagris M. (2019). Extremely efficient permutation and bootstrap hypothesis tests using R. To appear in the Journal of Modern Applied Statistical Methods.
https://arxiv.org/ftp/arxiv/papers/1806/1806.10947.pdf
welch.tests, trim.mean
boot.james( as.matrix(iris[1:25, 1:4]), as.matrix(iris[26:50, 1:4]) )
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.