tauboot <- function(x,data){
signs <- sign(outer(data[x,1], data[x,1], "-") * outer(data[x,2], 
data[x,2], "-"))
tau <- mean(signs[lower.tri(signs)])
tau
}