Try R: > s <- sample(c(3,3,3,3,3,4), size=100, replace=T) [...] output elided > table(s) 3 4 86 14 In order: this samples size times from a

Try R: > s <- sample(c(3,3,3,3,3,4), size=100, replace=T) [...] output elided > table(s) 3 4 86 14 In order: this samples size times from a vector created by c with 5 3s and one 4 (you could also say c(rep(3,5), 4) with replacement and assigns the sample to a vector called s. table(s) produces counts of unique values.