glhtbf_zzg2022 {HDNRA} | R Documentation |
Test proposed by Zhang et al. (2022)
Description
Zhang et al. (2022)'s test for general linear hypothesis testing (GLHT) problem for high-dimensional data under heteroscedasticity.
Usage
glhtbf_zzg2022(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
.
Zhang et al. (2022) proposed the following test statistic:
where with
, and
with
being the sample mean vectors.
They showed that under the null hypothesis, and a chi-squared-type mixture have the same normal or non-normal limiting distribution.
Value
A (list) object of S3
class htest
containing the following elements:
- p.value
the
-value of the test proposed by Zhang et al. (2022)
- statistic
the test statistic proposed by Zhang et al. (2022).
- beta
the parameters used in Zhang et al. (2022)'s test.
- df
estimated approximate degrees of freedom of Zhang et al. (2022)'s test.
References
Zhang J, Zhou B, Guo J (2022).
“Linear hypothesis testing in high-dimensional heteroscedastic one-way MANOVA: A normal reference -norm based test.”
Journal of Multivariate Analysis, 187, 104816.
doi:10.1016/j.jmva.2021.104816.
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_zzg2022(Y, G, n, p)