LC | R Documentation |
Given two sets of data matrices X and Y, where X is an n1 rows and p cols matrix and Y is an n2 rows and p cols matrix, we conduct hypothesis testing of the covariance matrix between two samples. The null hypothesis is:
H_0 : \Sigma_1 = \Sigma_2
\Sigma_1
and \Sigma_2
are the sample covariance matrices of X and Y respectively. This test method is based on the test method proposed by Li and Chen (2012). When the pval value is less than the significance coefficient (generally 0.05), the null hypothesis is rejected.
LC(X,Y)
X |
A matrix of n1 by p |
Y |
A matrix of n2 by p |
stat |
a test statistic value. |
pval |
a test p_value. |
Li, J. and Chen, S. X. (2012). Two sample tests for high-dimensional covariance matrices. The Annals of Statistics, 40(2):908-940.
## generate X and Y.
p= 500; n1 = 100; n2 = 150
X=matrix(rnorm(n1*p), ncol=p)
Y=matrix(rnorm(n2*p), ncol=p)
## run test
LC(X,Y)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.