Total TV Error Bound (Conditional)
compute_total_tv_bound.RdComputes the combined conditional TV bound using the triangle inequality: $$d_{TV}(\mathcal{L}(S_J | \alpha), \text{Poisson}(\alpha c_J)) \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.
In the submitted manuscript this is Appendix D, Corollary
cor:combined-conditional (Equations D14–D16). 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))
} # }