require(mvtnorm)
K <- 4
alpha <- 0.025
tVec <- (1:K)/K
boundaryVector <- primaryBoundaryVec(alpha,tVec,initIntvl=c(1,4),
OBF=TRUE,LanDeMets=FALSE,digits=3,printOut=TRUE)
boundaryVector <- primaryBoundaryVec(alpha,tVec,initIntvl=c(1,4),
OBF=FALSE,LanDeMets=FALSE,digits=3,printOut=TRUE)
boundaryVector <- primaryBoundaryVec(alpha,tVec,initIntvl=c(1,8),
OBF=TRUE,LanDeMets=TRUE,digits=3,printOut=TRUE)
boundaryVector <- primaryBoundaryVec(alpha,tVec,initIntvl=c(1,4),
OBF=FALSE,LanDeMets=TRUE,digits=3,printOut=TRUE)
Run the code above in your browser using DataLab