# Is the distribution of golf ball numbers uniform?
golfballs <- c(137, 138, 107, 104) # sample data; n = 486
rgolfballs <- do(1000) * table(rdata(486, 1:4)) # random data assuming uniform dist
# Now compare based on various test statistics
print(statTally(golfballs, rgolfballs, max))
print(statTally(golfballs, rgolfballs, min))
print(statTally(golfballs, rgolfballs, function(x) { diff( range(x)) } ))
print(statTally(golfballs, rgolfballs, sd))
Run the code above in your browser using DataLab