NS <- c(3,5,10,25)
B <- 10000
mean(sample(hf,n)-mean(sample(chow,n)))
})
})
B <- 10000
res <- sapply(NS, function(n){
sapply(1:B, function(j){mean(sample(hf,n)-mean(sample(chow,n)))
})
})
res <- sapply(NS, function(n){
sapply(1:B, function(j){