Paper VIReading edition

Prior-restricted sufficiency and filtering

Reading tools
Reading libraryFind statementsInspect Lean
Switch work

Saved reading

Read side by side

Start reading

Resources & source

Original TeX with line links · Download TeX

The wording and mathematics follow the retained manuscript. This reading format does not add a proof review or extend Lean coverage.

Skip contents and start reading

R3 Prior-restricted sufficiency and filtering

Let \(L\) be a strictly positive finite likelihood and let \({\mathcal B}_L(p)=L\odot p/{\left\langle L,p\right\rangle}\) for a probability vector \(p\). Fix probability vectors \(r^{(a)}\) supported on \(B_a\), and restrict the priors to \(p_i=q_a r_i^{(a)}\) for \(i\in B_a\), with \(q\) a coarse probability vector.

Proposition R3.1: Sufficiency relative to a prior disintegration

Define

\[ \overline L_a = \sum_{i\in B_a}L_ir_i^{(a)}. \]

For every prior in the declared family,

\[ C{\mathcal B}_L(p) = {\mathcal B}_{\overline L}(q). \]

Proof of Proposition R3.1

The evidence in block \(B_a\) is

\[ Z_a =q_a\sum_{i\in B_a}L_ir_i^{(a)} =q_a\overline L_a. \]

Apply Theorem E8.2.

Lemma R3.2: Normalization bound

Let \(a,b\in[0,\infty)^m\) be nonzero, and set \(A={\left\lVert a\right\rVert}_1\) and \(B={\left\lVert b\right\rVert}_1\). Then

\[ \boxed{ {\left\lVert a/A-b/B\right\rVert}_1 \leq \frac{2{\left\lVert a-b\right\rVert}_1}{\max\{A,B\}}. } \]

Proof of Lemma R3.2

Using \(A\) in the denominator gives

\[ {\left\lVert a/A-b/B\right\rVert}_1 \leq \frac{{\left\lVert a-b\right\rVert}_1}{A} + {\left\lVert b\right\rVert}_1{\left\lvert1/A-1/B\right\rvert}. \]

The second term equals \({\left\lvert B-A\right\rvert}/A\), which is at most \({\left\lVert a-b\right\rVert}_1/A\). Interchanging \(a\) and \(b\) gives the analogous bound with \(B\). Taking the smaller of the two bounds gives the result.

For an arbitrary prior \(p\), put \(q=Cp\) and \(Z(p,L)={\left\langle L,p\right\rangle}\). Define \(\ell_a(p)=\sum_{i\in B_a}L_ip_i/q_a\) when \(q_a>0\); choose any positive value when \(q_a=0\). Only the product \(q_a\ell_a(p)\) enters the formulas.

Theorem R3.3: Exact information loss formula

For every prior,

\[ \boxed{ C{\mathcal B}_L(p) = {\mathcal B}_{\ell(p)}(q). } \]

Let \(\overline L\) be any fixed positive coarse likelihood and set

\[ d(p\mid L,\overline L) = \sum_{a=1}^mq_a{\left\lvert\ell_a(p)-\overline L_a\right\rvert}. \]

Then

\[ {\left\lVert C{\mathcal B}_L(p)-{\mathcal B}_{\overline L}(q) \right\rVert}_1 \leq \frac{2d(p\mid L,\overline L)} {\max\{Z(p,L),{\left\langle\overline L,q\right\rangle}\}}. \]

Proof of Theorem R3.3

The packet evidence in block \(a\) is \(q_a\ell_a(p)\), which proves the exact identity. Apply Lemma R3.2 to

\[ a=q\odot\ell(p), \qquad b=q\odot\overline L. \]

Their \(\ell^1\) norms are the two evidences, and \({\left\lVert a-b\right\rVert}_1=d(p\mid L,\overline L)\).

Corollary R3.4: Oscillation bound

Suppose \(\overline L_a\) is any representative likelihood for block \(B_a\), and set

\[ \omega_a = \max_{i\in B_a}{\left\lvert L_i-\overline L_a\right\rvert}. \]

Then

\[ d(p\mid L,\overline L) \leq \sum_aq_a\omega_a. \]

Thus small within block likelihood oscillation gives a small coarse posterior defect whenever the evidence is bounded away from zero.

Proof of Corollary R3.4

The conditional average \(\ell_a(p)\) lies in the convex hull of the likelihood values inside \(B_a\). Hence

\[ {\left\lvert\ell_a(p)-\overline L_a\right\rvert}\leq\omega_a. \]

Multiply by \(q_a\) and sum.

Along the Bayesian flow \(p_i(t)=e^{tg_i}p_i(0)/\sum_j e^{tg_j}p_j(0)\), the following statement holds for every real observable \(h\).

Corollary R3.5: Covariance identity

For any observable \(h\in{\mathbb R}^n\),

\[ \frac{{\,\mathrm d}}{{\,\mathrm d} t}{\mathbb E}_{p(t)}[h] = {\operatorname{Cov}}_{p(t)}(h,g). \]

In particular,

\[ \frac{{\,\mathrm d}}{{\,\mathrm d} t}{\mathbb E}_{p(t)}[g] = {\operatorname{Var}}_{p(t)}(g)\geq0. \]

Proof of Corollary R3.5

Differentiate \(\sum_ih_ip_i(t)\) and insert the replicator equation.

For stochastic matrices \(K,\overline K\) and positive likelihoods \(L_y,\overline L_y\), put \(\Phi_y={\mathcal B}_{L_y}\circ K\) and \(\overline\Phi_y={\mathcal B}_{\overline L_y}\circ\overline K\).

Theorem R3.6: Exact coarse hidden Markov filter

Assume

\[ CK=\overline KC \]

and every likelihood \(L_y\) is constant on each coarse block with block values \(\overline L_y\). Then

\[ C\Phi_y(p)=\overline\Phi_y(Cp) \]

for every prior. For any observation sequence \(y_1,\ldots,y_r\), the complete filtering recursion commutes with aggregation:

\[ C\Phi_{y_r}\circ\cdots\circ\Phi_{y_1} = \overline\Phi_{y_r}\circ\cdots\circ\overline\Phi_{y_1}C. \]

Proof of Theorem R3.6

Exact Markov lumping gives \(CKp=\overline KCp\). The Bayesian sufficiency criterion gives exact commutation of the update. Their composition proves the one step identity. Iterate it over the observation sequence.

Theorem R3.7: One step filter defect

Let \(r=Kp\), \(q=Cr\), and let \(\ell_y(r)\) be the exact coarse likelihood obtained by averaging \(L_y\) inside each block with respect to \(r\). Assume

\[ 0<\overline m\leq(\overline L_y)_a\leq\overline M. \]

Then

\[ \begin{aligned} {\left\lVert C\Phi_y(p)-\overline\Phi_y(Cp)\right\rVert}_1 \leq{}& \frac{2d(r\mid L_y,\overline L_y)} {\max\{{\left\langle L_y,r\right\rangle},{\left\langle\overline L_y,q\right\rangle}\}}\\ &+ \frac{2\overline M}{\overline m} {\left\lVert CKp-\overline KCp\right\rVert}_1. \end{aligned} \]

Proof of Theorem R3.7

Insert the intermediate posterior \({\mathcal B}_{\overline L_y}(CKp)\). The first difference is bounded by Theorem R3.3. The second is bounded by Theorem E8.4 applied to the coarse likelihood.

Source and conversion record

Source SHA-256: 736e615131f57daf341b6cb0600761d58950450947c85d433206148cca228333

Unsupported prose commands stop the build. The mathematical expression sequence is checked against the source, with manuscript macros expanded for display.

Conversion records

Notation in Paper VI

Look up a symbol without leaving this chapter. Select text before opening help to search for it.

Search this work’s notation.

    Symbols can have different meanings in different passages. Read each result in its own context.

    Open the full notation index