function(x, alfa, beta, m)
{
	n <- length(x)
	cint <- jay(m, 2)
	for(i in 1:m) {
		xs <- sample(x, n, replace
			 = T)
		cint[i, 1] <- quantile(xs, 
			beta/2)
		cint[i, 2] <- quantile(xs, 1 -
			beta/2)
	}
	lower <- quantile(cint[, 1], alfa/2)
	upper <- quantile(cint[, 2], 1 - 
		alfa/2)
	res <- list(lower, upper)
	names(res) <- c("Lower", "Upper")
	res
}
