Poisson thinning

Poisson thinning splits a Poisson process by independently marking each arrival as retained or rejected. Let the original Homogeneous Poisson process have rate $\lambda >0$, and let $p(s)\in [0,1]$ be a deterministic measurable function of time. Attach independent uniform distributions $U_1,U_2,\dots$ on $(0,1)$, independent of the entire original Poisson process. At arrival time $T_i$, retain the arrival if $U_i\le p(T_i)$. Given all arrival times, these decisions are independent and have the specified retention probabilities. Define the two counting processes using an indicator function:The Poisson thinning theorem gives independent retained and rejected Inhomogeneous Poisson processes. Their rates are $\lambda p(s)$ and $\lambda (1-p(s))$, respectively, in the sense of the intensity function of a point process. In particular, for $0\le a<b$,The independence is of the two whole processes: every finite collection of retained counts is independent of every finite collection of rejected counts. It is not merely zero covariance for one interval. For constant $p$, the streams are Homogeneous Poisson processes of rates $\lambda p$ and $\lambda (1-p)$; the endpoints $p=0,1$ include the identically zero stream.^[ref]

For an exact toy calculation, take $\lambda =2$ per minute and retain with $p(s)=1/4$ on $(0,2]$ and $p(s)=3/4$ on $(2,4]$. In the first two minutes, the retained and rejected expected values are $2\cdot 2\cdot 1/4=1$ and $2\cdot 2\cdot 3/4=3$. In the next two minutes they are $3$ and $1$. Over all four minutes, $R_4,D_4$ are independent Poisson distributions with mean $4$ each, and their sum $N_4$ has mean $8$. ThusThese probabilities and expected values describe the model, not empirical measurements. In the diagram, the original nine synthetic arrivals split into five retained and four rejected arrivals; those realized counts need not equal their expected values.

The mechanism explains why mutually exclusive labels can yield independent counts: randomness of the original total offsets the dependence from allocating a fixed total. For constant $p$ and an interval with original mean $m$, the conditional retained count has a binomial distribution given $N=n$. The joint probability generating function isThis factorization proves the two independent Poisson distributions; disjoint original intervals give independent pairs. For varying $p(s)$, the Independent marking theorem for Poisson point processes gives the same conclusion with the two integrated intensities. The joint red and green Poisson count formula gives the corresponding direct calculation with probability mass functions.

Independence here is unconditional. After conditioning on $N=n$, the identity $D=n-R$ gives covariance $\mathrm{Cov}(R,D\mid N=n)=-np(1-p)$ for constant $p$. Thus the retained and rejected counts are dependent for $n>0$ and $0<p<1$, even though their unconditional covariance is zero. Poisson thinning can model independent routing of arrivals, or independent detection of arrivals with a known efficiency schedule; using only a retention fraction without checking the marking mechanism can give a wrong count law.

Independent marking is the hypothesis that guarantees the Poisson thinning theorem. Equal marginal retention probabilities do not guarantee it. Consider one Bernoulli random variable $B$ with success probability $p\in (0,1)$, independent of a Poisson process, and retain every arrival when $B=1$ and none when $B=0$. Each arrival is retained with probability $p$, but the marks share the same coin. In an interval whose original count has Poisson distribution of mean $m>0$, the retained count is $R=BN$ and the rejected count is $D=(1-B)N$, soThe extra variance shows that $R$ is not Poisson distributed; the negative covariance proves dependence. Its no-arrival probability is $1-p+p{e}^{-m}$, rather than $e^{-pm}$. For the toy values $p=1/2,m=2$, both rules have retained mean $1$, but the common-coin rule has variance $2$ and no-arrival probability $(1+e^{-2})/2\approx 0.567668$, compared with variance $1$ and no-arrival probability $e^{-1}\approx 0.367879$ for independent marking. These formulas follow by conditioning on $B$ and using the law of total variance.

An arrival-dependent or adversarial retention rule can change the law just as strongly. Keeping every second arrival gives $R_t=\lfloor N_t/2\rfloor$. Its successive retained gaps are sums of two independent exponential random variables, hence have the Erlang distribution of shape two and rate $\lambda$, not an exponential distribution. The resulting renewal process and its count variance are analyzed in Every second Poisson arrival has underdispersed counts, and the two dependent streams in alternating split of a Poisson process. As an even more restrictive rule, keeping only the first arrival ever gives $R_t=1_{\{N_t\ge 1\}}$, a Bernoulli distribution at each $t>0$, so it cannot have a positive-mean Poisson distribution. A thinning rule must specify how decisions depend on arrivals and on one another; the independent-marking conclusion cannot be applied automatically to arbitrary deletion.

Let $N$ be a nonnegative integer-valued random variable with finite expected value $\mu$, and independently label each of its items with either of two equally likely labels. If the two resulting counts are independent random variables, then $N$ has a Poisson distribution. Indeed, its probability generating function obeys $G(s)=G((1+s)/2{)}^2$. With $H(t)=G(1-t)$, iteration gives $H(t)=H(t/2^k{)}^{2^k}$. Finite expected value gives $H(h)=1-\mu h+o(h)$, hence $H(t)=e^{-\mu t}$ and $G(s)=e^{\mu (s-1)}$. Conversely, Poisson thinning produces independent Poisson distributions for the two counts.

Assign each of $N\sim \mathrm{Pois}(\lambda )$ claims with independent category labels to one of finitely many categories with probabilities $p_i$. Conditional counts have a multinomial distribution with probability generating function $({\sum }_i{p}_i{z}_i{)}^N$. Averaging over $N$ gives $\exp (\lambda ({\sum }_i{p}_i{z}_i-1))={\prod }_i\exp (\lambda p_i(z_i-1))$. This proves independent category counts and supports discrete-severity compound Poisson distribution representations.

Mark each member of an independent Poisson infection cohort by its symptom-onset interval, using independent delays. Poisson thinning gives independent onset-bin counts, and sums over independent cohorts remain independent Poisson counts. Fixed cohort totals instead induce multinomial dependence. The marking assumptions are essential.

Conditioned on a total of $n$ Poisson points, the number falling in a region of relative intensity $p$ is $\mathrm{Binomial}(n,p)$.

For an M-M-infinity queue with empty initial state, Poisson arrivals of rate $\lambda$ and exponential lifetimes of rate $\mu$ produce a Poisson occupancy of mean $\lambda (1-e^{-\mu t})/\mu$ at time $t$.