R/aeS.R

In meca7653/HDtest-2018-8: High Dimensional Hypothesis Testing for Mean Vectors, Covariance Matrices, and White Noise of Vector Time Series

aeS = function(ft, Sn, W, M, alpha){
  p = dim(ft)[1]
  n = dim(ft)[2]
  xi = t(mvrnorm(n = M, mu = rep(0, n), Sn))/sqrt(n)
  G1 = ft %*% xi
  G = W*G1
  rrr <- apply(abs(G),2, max)
  cv = quantile(rrr, 1-alpha)
  return(list(cv = cv, stat = rrr))
}