### Example
dat <- gendata_Fac(n = 300, p = 500)
res <- Factorm(dat$X)
X <- dat$X
# ex1: H01 is false
G1 <- 1:10; # all are nonzero loading vectors
FacRowMaxST(X, G1=G1, alpha=0.05, sub.frac=0.5)
FacRowMaxST(X, q= 6, G1=G1, alpha=0.05, sub.frac=0.5) # specify the true number of factors
# ex2: H01 is true
G1 <- 481:500 # all are zero loading vectors
FacRowMaxST(X, G1=G1, alpha=0.05, sub.frac=0.5)
FacRowMaxST(X, q= 7, G1=G1, alpha=0.05, sub.frac=0.5) # specify a false number of factors
Run the code above in your browser using DataLab