glhtbf_zz2022 {HDNRA} | R Documentation |
Test proposed by Zhang and Zhu (2022)
Description
Zhang and Zhu (2022)'s test for general linear hypothesis testing (GLHT) problem for high-dimensional data under heteroscedasticity.
Usage
glhtbf_zz2022(Y,G,n,p)
Arguments
Y |
A list of |
G |
A known full-rank coefficient matrix ( |
n |
A vector of |
p |
The dimension of data. |
Details
Suppose we have the following independent high-dimensional samples:
It is of interest to test the following GLHT problem:
where
is a
matrix collecting
mean vectors and
is a known full-rank coefficient matrix with
.
Let be the sample mean vectors and
be the sample covariance matrices.
Zhang and Zhu (2022) proposed the following U-statistic based test statistic:
where ,
, and
is the
th entry of the
matrix
.
Value
A (list) object of S3
class htest
containing the following elements:
- p.value
the
-value of the test proposed by Zhang and Zhu (2022).
- statistic
the test statistic proposed by Zhang and Zhu (2022).
- beta0
the parameter used in Zhang and Zhu (2022)'s test.
- beta1
the parameter used in Zhang and Zhu (2022)'s test.
- df
estimated approximate degrees of freedom of Zhang and Zhu (2022)'s test.
References
Zhang J, Zhu T (2022). “A new normal reference test for linear hypothesis testing in high-dimensional heteroscedastic one-way MANOVA.” Computational Statistics & Data Analysis, 168, 107385. doi:10.1016/j.csda.2021.107385.
Examples
set.seed(1234)
k <- 3
p <- 50
n <- c(25, 30, 40)
rho <- 0.1
M <- matrix(rep(0, k * p), nrow = k, ncol = p)
avec <- seq(1, k)
Y <- list()
for (g in 1:k) {
a <- avec[g]
y <- (-2 * sqrt(a * (1 - rho)) + sqrt(4 * a * (1 - rho) + 4 * p * a * rho)) / (2 * p)
x <- y + sqrt(a * (1 - rho))
Gamma <- matrix(rep(y, p * p), nrow = p)
diag(Gamma) <- rep(x, p)
Z <- matrix(rnorm(n[g] * p, mean = 0, sd = 1), p, n[g])
Y[[g]] <- Gamma %*% Z + t(t(M[g, ])) %*% (rep(1, n[g]))
}
G <- cbind(diag(k - 1), rep(-1, k - 1))
glhtbf_zz2022(Y, G, n, p)