if (FALSE) {
require(mvtnorm)
require(ldbounds)
K <- 6;
tVec <- c(140,328,453,578,659,1080)/1080;
alpha <- 0.025;
cvec.obf <- ldbounds::ldBounds(t=tVec,iuse=c(1),alpha=c(alpha),sides = 1);
cvec <- cvec.obf$upper.bounds;
alphaprime <- ldNominalSig(alpha,tVec,cvec,locPeak=4,type=2,
initIntvl=c(1,4))
}
Run the code above in your browser using DataLab