"BcInt"<- function(x,nboot,theta,...,alpha = c(.025,.05,.1,.9,.95,.975))
{

# adapted from bcanon( ) inside the bootstrap library by Mai Zhou 
# only do bias correction, no accelaration

call <- match.call()
n <- length(x)
thetahat <- theta(x,...)
bootsam <- t(matrix(sample(x,size=n*nboot,replace=T),ncol=nboot))

thetastar <- apply(bootsam,1,theta,...)
z0 <- qnorm(sum(thetastar<thetahat)/nboot)
zalpha <- qnorm(alpha)
tt <- pnorm(z0+z0+zalpha)
ooo <- trunc(tt*nboot)                   # Can we improve trunc
confpoints <- sort(thetastar)[ooo]       # so that take linear interpolation
confpoints <- cbind(alpha,confpoints)
dimnames(confpoints)[[2]] <- c("alpha","BcInt point")
return(confpoints, z0, call=call)
 
}