If I have a given set of character stats from rolling with "4d6, drop lowest", how can I calculate the probability of getting that specific result, or better values for each of them?
For example, one of the players in my campaign got 15, 15, 15, 14, 13, 13. In particular, any one of the values is lower than the highest expected value from rolling 4d6, drop lowest (16), but the average is very high. (The order of rolling does not matter, you can sort them highest to lowest or assign them to stats freely afterwards).
How likely is it to get that array or better (in the sense for each die separately, to get a result as high, or higher)?
I'd like to know this die by die, because typically, it is better to have some high dice and some low dice than a higher average of middling dice or just a higher total of points overall. (For example, the player here plays a wizard, and while high Dex and Con etc. are all nice, he really would have liked a 16 so he could get to 18 Int out of the gate).
How would a formula or anydice program look like, where I can just plug in the rolled numbers, and find out how likely it was to get them? (Anydice is my preferred tool, but other tools are OK too).