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