function(x, y, m)
{
	n <- length(x)
	acorr <- c(1:m)
	ns <- length(x)
	for(i in 1:m) {
		a <- sample(c(1:ns), n, 
			replace = T)
		ax <- x[a,  ]
		ay <- y[a,  ]
		acorr[i] <- cor(ax, ay)
	}
	acorr
}
