function(x, m)
{
	ft <- max(diag(var(x)))/min(diag(var(
		x)))
	n <- nrow(x)
	k <- ncol(x)
	fts <- c(1:m)
	for(i in 1:m) {
		xs <- jay(n, k)
		for(j in 1:k) {
			xs[, j] <- sample(x[,
				j], n, 
				replace = T)
		}
		fts[i] <- max(diag(var(xs)))/
			min(diag(var(xs)))
	}
	prop <- 0 * jay(m, 1)
	for(i in 1:m) {
		if(fts[i] >= ft)
			prop[i] <- 1
	}
	pvalue <- sum(prop)/m
	res <- list(ft, pvalue)
	names(res) <- c("FMAX", "P-Value")
	res
}
