In the original Dungeons & Dragons game (D&D), each character has six ability scores determined by rolling three six-sided dice and taking the sum (i.e., 3d6 taken 6 times). What is the most likely set of scores to be produced with this method?
A few breadcrumbs, I think:
- We can take the probability of any individual 3d6 result from a page like this.
- For any ordered collection of values, we can find the number of k-permutations via the formula for multinomial coefficients.
A few example probabilities:
- A set of all 10's: $0.125^6 \times 6!/6! = 3.8 \times 10^{-6}$
- Five 10's and one 11: $0.125^6 \times 6!/(5!1!) = 2.29 \times 10^{-5}$
- Three 10's and three 11's: $0.125^6 \times 6!/(3!3!) = 7.63 \times 10^{-5}$
- Three 10's, two 11's, and one 12: $0.125^5 \times 0.1157 \times 6!/(3!2!1!) = 0.212 \times 10^{-4}$
Note that each of these examples have sequentially increasing probabilities.
So, what set(s) of scores have the highest likelihood of occurring?
(Note that a few prior questions exist on SE Mathematics about probabilities of rolling particular ability scores for later editions of D&D, using fundamentally different mechanics, and are different queries: e.g., here and here.)