Computes the combined conditional TV bound using the triangle inequality: $$d_{TV}(K_J | \alpha, 1 + \text{NegBin}) \le B_{\text{Pois}} + B_{\text{lin}}$$
Usage
compute_total_tv_bound(J, alpha, cJ = log(J))Details
The result is capped at 1 since TV distance is bounded by 1.
From Lee (2026, Section 3.3, Theorem 1), the total conditional TV error decomposes as:
Poissonization error: \(S_J | \alpha\) vs \(\text{Poisson}(\lambda_J(\alpha))\)
Linearization error: \(\text{Poisson}(\lambda_J(\alpha))\) vs \(\text{Poisson}(\alpha c_J)\)
Examples
if (FALSE) { # \dontrun{
# Total bound at alpha = E[alpha] under Gamma(2, 1)
compute_total_tv_bound(J = 50, alpha = 2)
# Vectorized
compute_total_tv_bound(J = 50, alpha = c(0.5, 1, 2, 5))
} # }