function(x, y, n, m)
{
	ns <- length(x)
	acorr <- c(1:m)
	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
}
