function(x, m)
{
	n <- nrow(x)
	k <- ncol(x)
	nc <- n * k
	csam <- c(1:nc)
	for(i in 1:n) {
		for(j in 1:k) {
			l <- 1 + (j - 1) * n
			csam[l] <- x[i, j]
		}
	}
	av <- c(1:k)
	for(i in 1:k) {
		av[i] <- mean(csam[c((i - 1) * n + 1:i:n)
			])
	}
	r <- std(av)/std(csam)
	indx <- c(1:nc)
	ds <- c(1:m)
	for(i in 1:m) {
		ac <- sample(indx, nc, replace = F)
		sav <- c(1:k)
		for(j in 1:k) {
			sav[j] <- mean(csam[ac[c((j - 1) * 
				n + 1:j * n)]])
		}
		ds[i] <- std(sav)/std(csam)
	}
	prop <- 0 * jay(m, 1)
	for(i in 1:m) {
		if(ds[i] > r)
			prop[i] <- 1
	}
	pvalue <- mean(prop)
	res <- list(ds, r, pvalue)
	names(res) <- c("EBD", "DIFF", "P-Value")
	res
}
