View source: R/mean2_2014Thulin.R
mean2.2014Thulin | R Documentation |
Given two multivariate data X and Y of same dimension, it tests
H_0 : μ_x = μ_y\quad vs\quad H_1 : μ_x \neq μ_y
using the procedure by Thulin (2014) using random subspace methods. We did not enable parallel computing schemes for this in that it might incur huge computational burden since it entirely depends on random permutation scheme.
mean2.2014Thulin(X, Y, B = 100, nreps = 1000)
X |
an (n_x \times p) data matrix of 1st sample. |
Y |
an (n_y \times p) data matrix of 2nd sample. |
B |
the number of selected subsets for averaging. B≥q 100 is recommended. |
nreps |
the number of permutation iterations to be run. |
a (list) object of S3
class htest
containing:
a test statistic.
p-value under H_0.
alternative hypothesis.
name of the test.
name(s) of provided sample data.
thulin_highdimensional_2014SHT
## CRAN-purpose small example smallX = matrix(rnorm(10*3),ncol=10) smallY = matrix(rnorm(10*3),ncol=10) mean2.2014Thulin(smallX, smallY, B=10, nreps=10) # run the test ## Compare with 'mean2.2011LJW' ## which is based on random projection. n = 33 # number of observations for each sample p = 100 # dimensionality X = matrix(rnorm(n*p), ncol=p) Y = matrix(rnorm(n*p), ncol=p) ## run both methods with 100 permutations mean2.2011LJW(X,Y,nreps=100,method="m") # 2011LJW requires 'm' to be set. mean2.2014Thulin(X,Y,nreps=100)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.