Benjamini-Hochberg leave-one-out identity

Let $R$ be the number of rejections by the Benjamini-Hochberg procedure at level $\alpha$, and let $R_i$ be the rejection count after replacing $p_i$ by zero. Then $R_i\ge 1$ andIf $i$ was rejected, lowering its p-value does not change any ordered p-value beyond rank $R$, so $R_i=R$. Conversely, if $p_i\le \alpha R_i/m$, the other $R_i-1$ rejected values together with $p_i$ satisfy the original threshold; monotonicity gives $R=R_i$. For a valid true-null p-value independent of the others, conditioning on $R_i$ bounds the expectation of the right side by $\alpha /m$. Summing over true nulls proves false discovery rate at most $m_0\alpha /m$.

For the Benjamini-Hochberg procedure, replace a true-null p-value $P_i$ by zero and call the new rejection count $R^{(i)}$. The Benjamini-Hochberg leave-one-out identity identifies rejection with $R=r$ with the event $P_i\le \alpha r/m$ and $R^{(i)}=r$. If $P_i$ is uniform and independent of all the other p-values, each possible rejection count contributes $\alpha \mathbb{P}(R^{(i)}=r)/m$ to its expected false-discovery fraction. Summation gives $\alpha /m$ per true null and the displayed exact false discovery rate. Super-uniform independent nulls yield the corresponding upper bound.