if (FALSE) {
require(mvtnorm)
K <- 8
gammaVec <- sqrt((1:K)/K)
tVec <- gammaVec^2
alpha = 0.025
c <- 2.072274
cvec <- c/gammaVec
loc <- initLocPeak(alpha,tVec,cvec,type=2,initIntvl=c(1,4))
sbvec <- secondaryBoundary(alpha,tVec,cvec,loc,type=2,
initIntvl=c(1,8))
}
Run the code above in your browser using DataLab