random_draws_peterflom.R

makeDraws <- function(numberOfDraws) {
	i <- 1
	myOutput <- rep(NA, numberOfDraws) #preallocate the vector so it all goes faster
	while (i <= numberOfDraws) {
		a <- runif(7, 0, 100) #random draws from a uniform dist 0,100
		b <- runif(7, 0, 100)
		absDiff <- abs(a-b) #calculate the absolute diff
		myOutput[i] <- sum(absDiff) #sum up the diffs and put them in the output vector
		i <- i +1
	}
	myOutput #return the output vector
}

myDraws <- makeDraws(1e6) #start with 1mm draws

min(myDraws)
quantile(myDraws, .0001)