total <- sum(ants) dif <- NULL for (i in 1:1000) { #sample sample1 <- sample(ants,6) #dif dif[i] <- abs(mean(sample1)-((total-sum(sample1))/4)) }