Calculating the odds of dying when playing the snake-eyes game with a player base of arbitrary size.
For an arbitrary player base of B players, the maximum possible rounds that the game can run is then the whole number part of Log2(B), we will denote the maximum possible number of rounds as m.
We denote the probability of snake-eyes as p. In the case of the Daniel Reeves market p=1/36.
Let Dn= the probability density of the game ending in snake-eyes on round n then Dn=p(1−p)n−1.
The sum of the probability density of the games which end with snake eyes is
=m∑n=1Dn
=m∑n=1p(1−p)n−1
=pm∑n=1(1−p)n−1
=p1−(1−p)m1−(1−p)
=1−(1−p)m
The sum of the probability densities of the games ending in snake-eyes is less than 1 which means that the rounds ending in snake-eyes does not cover the full probability space.
Since there are finitely many dice rolls possible for our population there is always the possibility that the game ends and everyone won. This would require not rolling snake-eyes m times and has a probability density of Sm=(1−p)m which is exactly the term we need to sum to 1 giving us full coverage of the probability space.
Note: That Dm+Sm are the probability densities of the final round ending in a loss or ending in a win and must there for equal (1−p)m−1 the probability of getting to round m by the players winning the prior m−1 rounds.
Here we verify that Dm+Sm=(1−p)m−1
Dm=p(1−p)m−1
Sm=(1−p)m
Dm+Sm=p(1−p)m−1+(1−p)m
=p(1−p)m−1+(1−p)(1−p)m−1
=[p+(1−p)](1−p)m−1
=(1−p)m−1
Denote the number of players in each round n as Pn=2n−1
Denote the total number of players chosen in a game that ends at round n as Tn=2n−1
The probability of losing given that you are chosen to play can be computed by dividing the expected number of players who lost (red-eyed snakes) by the expected total number of players chosen.
First we will develop the numerator of our desired conditional probability. What is the expected number of players who lose in a game, given that cannot go beyond m rounds? This will be the sumproduct of the series of players in each round n and the series of probability densities of the game ending in round n
m∑n=1PnDn=m∑n=12n−1p(1−p)n−1=p×1−2m(1−p)m1−2(1−p)
Now we will develop the numerator of our desired conditional probability. What is the expected total number of people chosen to play? This will be the sumproduct of the series of total players in a game ending in round n and the series of probability densities of the game ending in round n.
And then we thank WolframAlpha for taking care of all the algebraic operations required to simply this fraction, see alternate form at : https://www.wolframalpha.com/input?i2d=true&i=Divide%5Bp*Divide%5B1-Power%5B2%2Cm%5DPower%5B%5C%2840%291-p%5C%2841%29%2Cm%5D%2C1-2%5C%2840%291-p%5C%2841%29%5D%2Cp%5C%2840%292*Divide%5B1-Power%5B2%2Cm-1%5DPower%5B%5C%2840%291-p%5C%2841%29%2Cm-1%5D%2C1-2%5C%2840%291-p%5C%2841%29%5D-Divide%5B1-Power%5B%5C%2840%291-p%5C%2841%29%2Cm-1%5D%2C1-%5C%2840%291-p%5C%2841%29%5D%5C%2841%29%2BPower%5B%5C%2840%291-p%5C%2841%29%2Cm-1%5D%5C%2840%29Power%5B2%2Cm%5D-1%5C%2841%29%5D
We find that the conditional probability of losing given that you are chosen to play from a arbitrary population of players is simply p. The size of the population was arbitrarily selected from the natural numbers so as the population goes to infinity the limit will also be p.
P(Lose|Played)=p∀n∈N therefor the limit of our conditional probability is limn→∞p=p
“The sum of the probability densities of the games ending in snake-eyes is less than
1 which means that the rounds ending in snake-eyes does not cover the full probability space.”
This is contradicted by the problem statement:
“At some point one of those groups will be devoured by snakes”, so there seems to be some error mapping the paradox to the math.
I posit that the supposition that “At some point one of those groups will be devoured by snakes” is erroneous. There exists a non-zero chance that the game goes on forever and infinitely many people win.
The issue is that the quoted supposition collapses the probability field to only those infinitely many universes where the game stops, but there is this one out of infinitely many universes where the game never stops and it has infinitely many winners so we end up with residual term of ∞∞. We cannot assume this term is zero just because of the ∞ in the denominator and disregard this universe.
“At some point one of those groups will be devoured by snakes” is erroneous
I wouldn’t say erroneous but I’ve added this clarification to the original question:
“At some point one of those groups will be devoured by snakes and then I stop” has an implicit “unless I roll snake eyes forever”. I.e., we are not conditioning on the game ending with snake eyes. The probability of an infinite sequences of non-snake-eyes is zero and that’s the sense in which it’s correct to say “at some point snake eyes will happen” but non-snake-eyes forever is possible in the technical sense of “possible”.
It sounds contradictory but “probability zero” and “impossible” are mathematically distinct concepts. For example, consider flipping a coin an infinite number of times. Every infinite sequence like HHTHTTHHHTHT… is a possible outcome but each one has probability zero.
So I think it’s correct to say “if I flip a coin long enough, at some point I’ll get heads” even though we understand that “all tails forever” is one of the infinitely many possible sequences of coin flips.
Going back to each of the finite cases, we can condition the finite case by the population that can support up to m rounds. iteration by the population size presupposes nothing about the game state and we can construct the Bayes Probability table for such games.
For a population M that supports at most m rounds of play the probability that a player will be in any given round n<=m is 2nM and the sum of the probabilities that a player is in round n from n=1→m is m∑n=12n−1M=2m−1M; We can let M=2m−1 because any additional population will not be sufficient to support an m+1 round until the population reaches 2m+1−1, which is just trading m for m+1 where we ultimately will take the limit anyhow.
The horizontal axis of the Bays Probability table now looks like this
[1M2M4M8M...2m−1M]
The vertical axis of the Bays Probability table we can independently look at the odds the game ends at round n for n<=m. This can be due to snake eyes or it can be due to reaching round m with out rolling snake eyes. For the rounds n=1→m where snake eyes were rolled the probability of the game ending on round n is p∗(1−p)n−1 and the probability that a reaches round m with out ever rolling snake eyes is (1−p)m . The sum of all of these possible end states in a game that has at most finite m rounds is (1−p)m+∑mn=1p(1−p)n−1 which equals =1
More over with the full Bayes Probability table we can find other conditional probabilities at infinity by taking the limit as the population grows to allow bigger and bigger maximum rounds of m.
Calculating the odds of dying when playing the snake-eyes game with a player base of arbitrary size.
For an arbitrary player base of B players, the maximum possible rounds that the game can run is then the whole number part of Log2(B), we will denote the maximum possible number of rounds as m.
We denote the probability of snake-eyes as p. In the case of the Daniel Reeves market p=1/36.
Let Dn= the probability density of the game ending in snake-eyes on round n then Dn=p(1−p)n−1.
The sum of the probability density of the games which end with snake eyes is
=m∑n=1Dn
=m∑n=1p(1−p)n−1
=pm∑n=1(1−p)n−1
=p1−(1−p)m1−(1−p)
=1−(1−p)m
The sum of the probability densities of the games ending in snake-eyes is less than 1 which means that the rounds ending in snake-eyes does not cover the full probability space.
Since there are finitely many dice rolls possible for our population there is always the possibility that the game ends and everyone won. This would require not rolling snake-eyes m times and has a probability density of Sm=(1−p)m which is exactly the term we need to sum to 1 giving us full coverage of the probability space.
Note: That Dm+Sm are the probability densities of the final round ending in a loss or ending in a win and must there for equal (1−p)m−1 the probability of getting to round m by the players winning the prior m−1 rounds.
Here we verify that Dm+Sm=(1−p)m−1
Dm=p(1−p)m−1
Sm=(1−p)m
Dm+Sm=p(1−p)m−1+(1−p)m
=p(1−p)m−1+(1−p)(1−p)m−1
=[p+(1−p)](1−p)m−1
=(1−p)m−1
Denote the number of players in each round n as Pn=2n−1
Denote the total number of players chosen in a game that ends at round n as Tn=2n−1
The probability of losing given that you are chosen to play can be computed by dividing the expected number of players who lost (red-eyed snakes) by the expected total number of players chosen.
First we will develop the numerator of our desired conditional probability. What is the expected number of players who lose in a game, given that cannot go beyond m rounds? This will be the sumproduct of the series of players in each round n and the series of probability densities of the game ending in round n
m∑n=1PnDn=m∑n=12n−1p(1−p)n−1=p×1−2m(1−p)m1−2(1−p)
Now we will develop the numerator of our desired conditional probability. What is the expected total number of people chosen to play? This will be the sumproduct of the series of total players in a game ending in round n and the series of probability densities of the game ending in round n.
(m∑n=1DnTn)+SmTm=(m−1∑n=1DnTn)+(DM+Sm)Tm=(m−1∑n=1DnTn)+(1−p)m−1Tm
=(m−1∑n=1p(1−p)n−1(2n−1))+(1−p)m−1(2m−1)
=p(m−1∑n=12n(1−p)n−1−m−1∑n=1(1−p)n−1)+(1−p)m−1(2m−1)
=p(2m−1∑n=12n−1(1−p)n−1−m−1∑n=1(1−p)n−1)+(1−p)m−1(2m−1)
=p(2×1−2m−1(1−p)m−11−2(1−p)−1−(1−p)m−11−(1−p))+(1−p)m−1(2m−1)
Combining the numerator and denominator into the fraction we seek we now have
E(Lose)E(Played)=∑mn=1PnDn(∑mn=1DnTn)+SmTm=p×1−2m(1−p)m1−2(1−p)p(2×1−2m−1(1−p)m−11−2(1−p)−1−(1−p)m−11−(1−p))+(1−p)m−1(2m−1)
Note: that since it is impossible to lose unless you play the conditional probability is equal to the above ratio.
P(Lose|Played)=P(Lose∩Played)P(Played)=P(Lose)P(Played)=E(Lose)E(Played)
And then we thank WolframAlpha for taking care of all the algebraic operations required to simply this fraction, see alternate form at : https://www.wolframalpha.com/input?i2d=true&i=Divide%5Bp*Divide%5B1-Power%5B2%2Cm%5DPower%5B%5C%2840%291-p%5C%2841%29%2Cm%5D%2C1-2%5C%2840%291-p%5C%2841%29%5D%2Cp%5C%2840%292*Divide%5B1-Power%5B2%2Cm-1%5DPower%5B%5C%2840%291-p%5C%2841%29%2Cm-1%5D%2C1-2%5C%2840%291-p%5C%2841%29%5D-Divide%5B1-Power%5B%5C%2840%291-p%5C%2841%29%2Cm-1%5D%2C1-%5C%2840%291-p%5C%2841%29%5D%5C%2841%29%2BPower%5B%5C%2840%291-p%5C%2841%29%2Cm-1%5D%5C%2840%29Power%5B2%2Cm%5D-1%5C%2841%29%5D
We find that the conditional probability of losing given that you are chosen to play from a arbitrary population of players is simply p. The size of the population was arbitrarily selected from the natural numbers so as the population goes to infinity the limit will also be p.
P(Lose|Played)=p ∀ n∈N therefor the limit of our conditional probability is limn→∞p=p
“The sum of the probability densities of the games ending in snake-eyes is less than 1 which means that the rounds ending in snake-eyes does not cover the full probability space.”
This is contradicted by the problem statement: “At some point one of those groups will be devoured by snakes”, so there seems to be some error mapping the paradox to the math.
I posit that the supposition that “At some point one of those groups will be devoured by snakes” is erroneous. There exists a non-zero chance that the game goes on forever and infinitely many people win.
The issue is that the quoted supposition collapses the probability field to only those infinitely many universes where the game stops, but there is this one out of infinitely many universes where the game never stops and it has infinitely many winners so we end up with residual term of ∞∞. We cannot assume this term is zero just because of the ∞ in the denominator and disregard this universe.
I wouldn’t say erroneous but I’ve added this clarification to the original question:
It sounds contradictory but “probability zero” and “impossible” are mathematically distinct concepts. For example, consider flipping a coin an infinite number of times. Every infinite sequence like HHTHTTHHHTHT… is a possible outcome but each one has probability zero.
So I think it’s correct to say “if I flip a coin long enough, at some point I’ll get heads” even though we understand that “all tails forever” is one of the infinitely many possible sequences of coin flips.
I think that’s what makes this a paradox.
Going back to each of the finite cases, we can condition the finite case by the population that can support up to m rounds. iteration by the population size presupposes nothing about the game state and we can construct the Bayes Probability table for such games.
For a population M that supports at most m rounds of play the probability that a player will be in any given round n<=m is 2nM and the sum of the probabilities that a player is in round n from n=1→m is m∑n=12n−1M=2m−1M; We can let M=2m−1 because any additional population will not be sufficient to support an m+1 round until the population reaches 2m+1−1, which is just trading m for m+1 where we ultimately will take the limit anyhow.
The horizontal axis of the Bays Probability table now looks like this
[1M2M4M8M...2m−1M]
The vertical axis of the Bays Probability table we can independently look at the odds the game ends at round n for n<=m. This can be due to snake eyes or it can be due to reaching round m with out rolling snake eyes. For the rounds n=1→m where snake eyes were rolled the probability of the game ending on round n is p∗(1−p)n−1 and the probability that a reaches round m with out ever rolling snake eyes is (1−p)m . The sum of all of these possible end states in a game that has at most finite m rounds is (1−p)m+∑mn=1p(1−p)n−1 which equals =1
So we have m+1 rows for the horizontal axis
⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣pp(1−p)p(1−p)2p(1−p)3⋮p(1−p)m(1−p)m⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦
So the Bayes Probability Table starts to look like this in general.
⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣[1M2M4M8M...2m−1M]⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣pp(1−p)p(1−p)2p(1−p)3⋮p(1−p)m(1−p)m⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦⎡⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢⎣pM2pM4pM8pM...p∗2m−1Mp(1−p)M2p(1−p)M4p(1−p)M8p(1−p)M...p(1−p)∗2m−1Mp(1−p)2M2p(1−p)2M4p(1−p)2M8p(1−p)2M...p(1−p)2∗2m−1Mp(1−p)3M2p(1−p)3M4p(1−p)3M8p(1−p)3M...p(1−p)3∗2m−1M⋮⋮⋮⋮⋱⋮p(1−p)m−1M2p(1−p)m−1M4p(1−p)m−1M8p(1−p)m−1M...p(1−p)m−1∗2m−1M(1−p)mM2(1−p)mM4(1−p)mM8(1−p)mM...(1−p)m∗2m−1M⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦⎤⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥⎦
The total probability of losing is equal to the sum of the diagonal where i=j
pM+2p(1−p)M+4p(1−p)2M+8p(1−p)3M+...+p(1−p)m−1∗2m−1M=pMm∑i=1[(1−p)i−1∗2i−1]
Total probability of being chosen is the sum of the diagonal and all the cells below the diagonal.
m∑i=1[2i−1∗(1−p)mM]+m∑i=1i∑j=1[2j−1∗p(1−p)i−1M]
=(2m−1)(1−p)mM+pMm∑i=1[(2i−1)(1−p)i−1]
So the conditional probability of losing given that you have been selected is
pM∑mi=1[(1−p)i−1∗2i−1](2m−1)(1−p)mM+pM∑mi=1[(2i−1)(1−p)i−1]
=p∑mi=1[(1−p)i−1∗2i−1](2m−1)(1−p)m+p∑mi=1[(2i−1)(1−p)i−1]=p
More over with the full Bayes Probability table we can find other conditional probabilities at infinity by taking the limit as the population grows to allow bigger and bigger maximum rounds of m.
So the conditional probability of losing given that you have been selected from an infinite population
Such as the conditional probability of being chosen in the game that has no snake eyes given you were chosen at all.
Or the conditional probability of dying given that you were chosen in precisely round k
This seems like great work! If we’re allowing to run out of players, the whole paradox collapses.