4

My approach so far is as follows:

If $X_{n-2}$ is the number of heads in $n-2$ tosses then $E[X_{n-2}] = \frac{n-2}{2}$. So up to the $n-2$ toss it is all the same regardless if we end our experiment in $n$ or $n-1$ tosses. Then if we want to ensure that we get more heads in $n-1$ than in $n$ tosses:

$P(X_{n-1} > X_{n}) = \frac{1}{2} \frac{1}{2.2} = \frac{1}{8}$

where $\frac{1}{2}$ is the probability of obtaining heads in the last toss out of $n-1$ and $\frac{1}{2.2}$ is the probability of obtaining two tails in the last two of $n$ tosses ?

Does this make any sense ? And what is the general approach for problems like this ?

Thanks in advance!

2 Answers2

5

Consider $X$ and $Y$ i.i.d. binomial $(n-1,\frac12)$ and $Z$ Bernoulli $\frac12$ independent of $(X,Y)$, then the desired probability is

$$p=P(X>Y+Z)$$

To compute $p$, note that the event $[X\leqslant Y+Z]$ occurs when either $X\leqslant Y$, or $X=Y+1$ and $Z=1$, thus, $$P(X\leqslant Y+Z)=P(X\leqslant Y)+\frac12P(X=Y+1)$$ Using the fact that $P(X<Y)+P(X=Y)+P(X>Y)=1$ by the law of total probability and that $P(X<Y)=P(X>Y)$ by symmetry, one gets $$2p=1-P(X=Y)-P(X=Y+1)$$ Now, $$P(X=Y)=\sum_kP(X=k)P(Y=k)=\frac1{4^n}\sum_k{n\choose k}^2$$ and a classical trick yields $$\sum_k{n\choose k}^2=\sum_k{n\choose k}{n\choose n-k}=[t^n](1+t)^n(1+t)^n={2n\choose n}$$ Likewise, $$P(X=Y+1)=\sum_kP(X=k+1)P(Y=k)=\frac1{4^n}\sum_k{n\choose k+1}{n\choose k}$$ and the same trick yields $$\sum_k{n\choose k+1}{n\choose k}=\sum_k{n\choose k}{n\choose n-k-1}=[t^{n-1}](1+t)^n(1+t)^n={2n\choose n-1}$$ Thus, $$2p=1-\frac1{4^n}{2n\choose n}-\frac1{4^n}{2n\choose n-1}$$ or, equivalently,

$$p=\frac12-\frac1{2^{2n+1}}\frac{2n+1}{n+1}{2n\choose n}$$

In terms of Catalan numbers, this reads $$p=\frac12-\frac{2n+1}{2^{2n+1}}C_n$$

Did
  • 279,727
3

I do not think that your approach is correct. A general approach to a problem like this would be the following.

We know that the binomial distribution counts the number of heads in $n$ coin tosses. Let $X_{n} \sim Bin(n,\frac{1}{2})$ (so the number of heads in $n$ coin tosses). You want the probability that $X_{n-1} > X_{n}$. This can be written as \begin{align*} & \mathbb{P}(X_{n-1} > X_{n}) \\ = & \sum_{k=1}^{n-1} \mathbb{P}(X_{n-1} > X_{n} \mid X_{n-1} = k) \mathbb{P}(X_{n-1} = k) \\ = & \sum_{k=1}^{n-1} \mathbb{P}(X_{n} < k) \mathbb{P}(X_{n-1} = k) \\ = & \sum_{k=1}^{n-1} \left[ \sum_{m = 0}^{k-1} \binom{n}{m} \left(\frac{1}{2}\right)^{n} \right] \binom{n-1}{k} \left(\frac{1}{2}\right)^{n-1}\\ = & \left(\frac{1}{2}\right)^{2n-1} \sum_{k=1}^{n-1} \sum_{m = 0}^{k-1} \binom{n}{m} \binom{n-1}{k} \end{align*} I do not know how this sum can simplified or if this even is possible.

Filling in some small $n$ we get that $\mathbb{P}(X_{1} > X_{2}) = \frac{1}{8}$ and $\mathbb{P}(X_{2} > X_{3}) = \frac{3}{16}$

  • $\frac1{2} \Sigma C(k)/4^k$ – karakfa Sep 20 '18 at 20:02
  • 1
    This is the obvious thing to do +1. The rather cryptic comment made by @karakfa indicates that an equivalent formulation is $$\sum_{k=1}^{n-1} \frac{C_k}{2^{2k+1}},$$ where $C_k$ is the $k^{\text{th}}$ Catalan number. That's easy to verify experimentally and I assume that some binomial identities could transform your answer into karakfa's. Perhaps karakfa has a more fundamental justification? – Mark McClure Sep 20 '18 at 20:37
  • I don't think that the result written by @karakfa is correct. If you try it for $n = 2017$ for example you will get probability of $0.987$ which seems fishy to me ... – Vasil Yordanov Sep 20 '18 at 22:04
  • for large $n$ it should converge to $\frac{1}2$ – karakfa Sep 20 '18 at 23:52
  • @VasilYordanov If I flip 2017 times and you flip 2016 times, then there should scarcely be a difference. So, clearly, the result should converge to 1/2, which is what I get. I suspect you started your sum at $k=0$, rather than $k=1$? – Mark McClure Sep 21 '18 at 01:17
  • 1
    I believe that the answer given from @user585522 can be simplified a bit: $$\sum_{k=1}^{n-1} \left[ \sum_{m = 0}^{k-1} \binom{n}{m} \left(\frac{1}{2}\right)^{n} \right] \binom{n-1}{k} \left(\frac{1}{2}\right)^{n-1} = \$$

    $$\left(\frac{1}{2}\right)^{2n-1} (2^n-1) \sum_{m=0}^{k-1} \binom{n}{m}$$

    – Vasil Yordanov Sep 21 '18 at 06:19