I am very curious about simulation of an event where an estimating/sampling technique is used.
In this example, the goal is to simulate a subset of all the roughly $2.6$ million $5$ card poker hands from a standard $52$ card deck and determine how many are $2$ pair hands.
We already know that the correct count is about $123,552 / 2,598,960$ which is about $4.75$% but I am using this example because it has a reasonable # of outcomes to be simulated in any proportion from $0$ to $100$% and because we can use it to check the accuracy of the partial simulation easily.
I would like if someone could simulate a subset of these hands using some simulation software and perhaps build a small table with 5 columns, namely:
- Number of hands evaluated
- % of total possible hands
- Number of winners found
- Extrapolated winners expected
- Actual error
So for example, a good starting point would be to simulate $1$% of the roughly $2.6$ million possible outcomes so that would be about $26,000$. So for our table we would fill in the following:
- ? (we would expect around $1236$)
- $100$ * whatever we actually get in # $3$.
- how far off are we from the expected $1236$ in %?
I'd also be curious of what type of sampling method you used such as random or something else. For example, would you just choose $5$ random cards $26,000$ times and check for a 2 pair? Would that be the best method in this case or would something other than random be more appropriate?
I would like to see the results of some other sample sizes too like $0.01$% ($260$ hands), $0.1$% ($2600$ hands) and $2$% ($52,000$ hands).
Can Mathematica and/or other simulation software handle this problem easily? If so, does it automatically check that the same random hand is not generated multiple times for the same simulation or does that require separate programming or is it not a concern since it would be rare if the number of samples is small (such as $1$% or less of the total # of possible $5$ card hands)?