Tutorials › Stochastic Processes › Retained Counts Conditional on the Original Count

Properties of Poisson Processes · Tutorial 422 of 994

Retained Counts Conditional on the Original Count

Learn why fixing the original count gives a binomial retained count, while allowing that count to vary produces a Poisson law.

Intermediate 10 min read

What You'll Learn

  • Define the retained count using independent Bernoulli marks on Poisson arrivals.
  • Derive its conditional binomial distribution, including endpoint cases.
  • Compute conditional means and variances from the Bernoulli indicators.
  • Derive the unconditional Poisson distribution by mixing over the original count.
  • Distinguish conditioning on the original count from conditioning on the retained count.
  • Identify and correct a common error that treats a conditional count as Poisson.

Fixing the Original Count Changes the Distribution

Suppose events arrive according to a Poisson process, and each event is independently retained with a fixed probability. The retained count has two different descriptions, depending on what is known. If the original number of arrivals is fixed at \(n\), the retained count is binomial: it counts successes in \(n\) independent retention trials. If the original count is allowed to vary, averaging over its Poisson distribution gives a Poisson retained count. These statements concern different probability laws, not competing answers to the same conditioning question.

Work on a probability space \((\Omega,\mathcal{F},\mathbb{P})\). Let \((N_t)_{t\ge0}\) be a Poisson process with rate \(\lambda\ge0\). Write \(S_k\) for its \(k\)th arrival time, and let \(I_k\) be the indicator that its \(k\)th arrival is retained. Assume that the \(I_k\) are independent Bernoulli random variables with \(\mathbb{P}(I_k=1)=p\), where \(0\le p\le1\), and that the entire sequence of marks is independent of the Poisson process. The retained count by time \(t\) is \(R(t)=\sum_{k=1}^{N_t} I_k\), with an empty sum equal to zero. For definiteness, the information available through time \(t\) is the filtration \(\mathcal{F}_t=\sigma(N_s,I_k\mathbf{1}_{\{S_k\le s\}}:0\le s\le t,\ k\ge1)\). Set \(\theta=\lambda t\). When \(\theta>0\), \(N_t\) has the Poisson distribution with mean \(\theta\), and every event \(\{N_t=n\}\), \(n\ge0\), has positive probability.

Theorem (Conditional retained-count law): For \(\theta>0\), \(n\ge0\), and \(0\le p\le1\), conditional on \(N_t=n\), \(R(t)\) has the binomial distribution with parameters \(n\) and \(p\). Thus, for \(r=0,\ldots,n\), \(\mathbb{P}(R(t)=r\mid N_t=n)=\binom{n}{r}p^r(1-p)^{n-r}\). For \(r\) outside this range the probability is zero. At \(p=0\) the retained count is zero, and at \(p=1\) it equals \(n\).

Proof. On the event \(N_t=n\), the retained count is the sum \(I_1+\cdots+I_n\). Independence of the marks from \(N_t\) means that conditioning on \(N_t=n\) leaves these \(n\) marks independent, each with probability \(p\) of being one. For any particular set of \(r\) indices to be retained and the other \(n-r\) indices to be rejected, the conditional probability is \(p^r(1-p)^{n-r}\). There are \(\binom nr\) such sets, and their events are disjoint. Summing their probabilities gives the displayed formula. If \(r<0\) or \(r>n\), no such set exists. When \(p=0\), all \(n\) indicators equal zero with probability one; when \(p=1\), all equal one. These direct cases also cover the endpoint interpretations of the formula. \(\square\)

This is a conditional distribution: it describes uncertainty in the marks after the original count has been fixed. In particular, the number of possible retained counts is limited to \(0,\ldots,n\). Its conditional mean and variance follow from the same sum of indicators.

Conditional moments: Whenever \(\theta>0\) and \(n\ge0\), \(\mathbb{E}[R(t)\mid N_t=n]=np\) and \(\operatorname{Var}(R(t)\mid N_t=n)=np(1-p)\).

Proof. On \(\{N_t=n\}\), write \(R(t)=\sum_{k=1}^n I_k\). Each indicator has conditional mean \(p\) and conditional variance \(p(1-p)\), because the marks are independent of \(N_t\). Linearity of conditional expectation gives \(\mathbb{E}[R(t)\mid N_t=n]=\sum_{k=1}^n p=np\). Conditional independence makes the covariance between distinct indicators zero, so the variance of the sum is \(\sum_{k=1}^n p(1-p)=np(1-p)\). If \(n=0\), the sum is empty and both moments are zero, as the formulas state. \(\square\)

Worked Example: A Retained Count Given Four Arrivals

Suppose \(\lambda=\tfrac32\), \(t=2\), and \(p=\tfrac13\). Find the probability of exactly two retained events given that four original events occurred, and compute the conditional mean and variance.

First, \(\theta=\lambda t=(\tfrac32)(2)=3\), so \(\mathbb{P}(N_t=4)=e^{-3}3^4/4!>0\), and conditioning on \(N_t=4\) is meaningful. The conditional law is binomial with \(n=4\) and \(p=\tfrac13\). Therefore

$$ \mathbb{P}(R(t)=2\mid N_t=4) =\binom42\left(\frac13\right)^2\left(\frac23\right)^2 =6\cdot\frac19\cdot\frac49 =\frac{24}{81} =\frac{8}{27}. $$

The conditional mean is \(np=4\cdot\tfrac13=\tfrac43\). The conditional variance is \(np(1-p)=4\cdot\tfrac13\cdot\tfrac23=\tfrac89\). As a check on the probability calculation, the four conditional probabilities, for \(r=0,1,2,3,4\), are respectively \(\tfrac{16}{81},\tfrac{32}{81},\tfrac{24}{81},\tfrac{8}{81},\tfrac{1}{81}\). Their sum is \(\tfrac{16+32+24+8+1}{81}=1\), and the probability at \(r=2\) agrees with the calculation above.

Allowing the Original Count to Vary

Without conditioning on \(N_t\), the retained count is a mixture over all possible original counts. Each term in that mixture uses the conditional binomial law, weighted by the probability of the corresponding Poisson count. The resulting distribution is not binomial with a fixed number of trials: the number of trials itself is random.

Theorem (Unconditional retained-count law): If \(\theta=\lambda t>0\), then \(R(t)\) has a Poisson distribution with mean \(\theta p\). That is, for \(r\ge0\), \(\mathbb{P}(R(t)=r)=e^{-\theta p}(\theta p)^r/r!\). If \(\theta=0\) or \(p=0\), then \(R(t)=0\) with probability one. If \(p=1\), then \(R(t)=N_t\).

Proof. First take \(\theta>0\) and \(0<p<1\). For \(r\ge0\), the event \(R(t)=r\) can occur only when \(N_t=n\) for some \(n\ge r\). Conditioning on \(N_t\), and using the Poisson probability for \(N_t=n\) and the conditional binomial probability, gives \(\mathbb{P}(R(t)=r)=\sum_{n=r}^{\infty}e^{-\theta}\frac{\theta^n}{n!}\binom nr p^r(1-p)^{n-r}\). Set \(m=n-r\), so \(n!=(m+r)!\) and \(\binom{m+r}{r}=(m+r)!/(r!m!)\). Substitution and cancellation yield \[ \mathbb{P}(R(t)=r) =e^{-\theta}\frac{(\theta p)^r}{r!} \sum_{m=0}^{\infty}\frac{(\theta(1-p))^m}{m!} =e^{-\theta}\frac{(\theta p)^r}{r!}e^{\theta(1-p)} =e^{-\theta p}\frac{(\theta p)^r}{r!}. \] The series identity used here is the power-series expansion of the exponential function. If \(\theta=0\), there are no arrivals by time \(t\), so \(R(t)=0\). If \(p=0\), no arrival is retained. If \(p=1\), every arrival is retained, so \(R(t)=N_t\), which has mean \(\theta=\theta p\). These cases establish the result at the boundaries as well. \(\square\)

The conditional and unconditional laws have different supports and parameters. Given \(N_t=n\), the count cannot exceed \(n\), and its mean is \(np\). Unconditionally, there is no fixed upper bound on the retained count when \(\theta p>0\); its Poisson mean is \(\theta p\). The unconditional result follows by averaging over the original count, not by replacing \(n\) in the conditional formula with its mean.

Worked Example: Conditional and Unconditional Probabilities

Let \(\lambda t=2\) and \(p=\tfrac12\). Compare the probability of two retained events given \(N_t=4\) with the unconditional probability of two retained events. Also compute the conditional and unconditional means and variances.

Conditionally on \(N_t=4\), \(R(t)\) is binomial with parameters \(4\) and \(\tfrac12\), so \(\mathbb{P}(R(t)=2\mid N_t=4)=\binom42(\tfrac12)^2(\tfrac12)^2=6\cdot\tfrac1{16}=\tfrac38\). The conditional mean is \(4\cdot\tfrac12=2\), and the conditional variance is \(4\cdot\tfrac12\cdot\tfrac12=1\).

Unconditionally, the retained count is Poisson with mean \((\lambda t)p=2\cdot\tfrac12=1\). Hence \(\mathbb{P}(R(t)=2)=e^{-1}1^2/2!=e^{-1}/2\), and both its mean and variance are \(1\). In particular, the unconditional probability is not the conditional probability: the former uses the full Poisson mixture, whereas the latter fixes the original count at four.

As a further check, the conditional moments average to the unconditional mean and variance through the laws of total expectation and total variance. Since \(N_t\) is Poisson with mean and variance \(2\), \(\mathbb{E}[R(t)]=\mathbb{E}[\mathbb{E}(R(t)\mid N_t)]=\mathbb{E}[N_t/2]=1\). For the variance, \(\operatorname{Var}(R(t))=\mathbb{E}[\operatorname{Var}(R(t)\mid N_t)]+\operatorname{Var}(\mathbb{E}[R(t)\mid N_t])=\mathbb{E}[N_t/4]+\operatorname{Var}(N_t/2)=2/4+2/4=1\). This matches the unconditional Poisson calculation.

Reversing the Conditioning

Conditioning on the original count and conditioning on the retained count are not interchangeable. The first leaves \(n\) Bernoulli retention decisions. The reverse question asks how many original arrivals could have produced the observed retained count, and it allows any number of rejected arrivals. The next calculation makes that distinction explicit.

Worked Example: Original Count Given One Retained Event

Take \(\lambda t=2\) and \(p=\tfrac12\). Find the conditional distribution of \(N_t\) given \(R(t)=1\), and compare it with the binomial law in the forward direction.

The event \(R(t)=1\) has positive probability: the unconditional retained count is Poisson with mean \(1\), so \(\mathbb{P}(R(t)=1)=e^{-1}\). For \(n\ge1\), Bayes' formula gives \(\mathbb{P}(N_t=n\mid R(t)=1)=\mathbb{P}(N_t=n)\mathbb{P}(R(t)=1\mid N_t=n)/\mathbb{P}(R(t)=1)\). Substitute the Poisson probability with mean \(2\) and the binomial probability with parameters \(n,\tfrac12\): \[ \mathbb{P}(N_t=n\mid R(t)=1) =\frac{e^{-2}2^n/n!\ \cdot\ n(1/2)^n}{e^{-1}} =\frac{e^{-1}}{(n-1)!}. \] For \(n=0\), this conditional probability is zero because one retained event requires at least one original event.

Thus \(N_t-1\), given \(R(t)=1\), has a Poisson distribution with mean \(1\): writing \(m=n-1\ge0\), the probability is \(e^{-1}/m!\). The probabilities sum to \(e^{-1}\sum_{m=0}^{\infty}1/m!=1\). For example, \(\mathbb{P}(N_t=3\mid R(t)=1)=e^{-1}/2\), obtained by setting \(n=3\) in the formula. This reverse conditional law is not the binomial distribution of \(R(t)\) given \(N_t=n\).

A Common Error and Its Correction

A frequent mistake is to treat the conditional retained count as Poisson with mean \(np\). For instance, with \(N_t=2\) fixed and \(p=\tfrac12\), the correct conditional probability of no retained events is \(\mathbb{P}(R(t)=0\mid N_t=2)=(1-\tfrac12)^2=\tfrac14\). If one incorrectly substitutes the conditional mean \(np=1\) into a Poisson formula, one obtains \(e^{-1}\) instead. The error is that a Poisson law with mean \(np\) describes neither the fixed-\(n\) retention experiment nor its conditional distribution. The correct conditional law is binomial; the Poisson law arises only after averaging over a Poisson-distributed original count.

When solving a retention problem, first identify what is being conditioned on. If the original count is specified as \(n\), use the binomial probability and restrict the retained count to \(0\le r\le n\). If no original count is specified, average the binomial probabilities over the Poisson law, or use the resulting Poisson distribution with mean \(\lambda t p\). Keeping those two steps separate prevents confusing a conditional law with its mixture.

Check Your Understanding

For each question, identify whether the original count is fixed or random before choosing a distribution.

  1. Given \(N_t=5\) and retention probability \(p=\tfrac25\), write the probability that exactly two arrivals are retained.
  2. If \(\lambda t=3\) and \(p=\tfrac13\), state the unconditional distribution, its mean, and the probability of zero retained arrivals.
  3. For fixed \(N_t=n\), explain why \(R(t)\) cannot exceed \(n\), even when \(\lambda t\) is large.
  4. When \(\lambda t=2\) and \(p=\tfrac12\), compute \(\mathbb{E}[R(t)\mid N_t=4]\) and compare it with \(\mathbb{E}[R(t)]\).
  5. Explain why the conditional law of \(R(t)\) given \(N_t=n\) does not directly give the law of \(N_t\) given \(R(t)=r\).