The Living Thing / Notebooks :

Nearly sufficient statistics

:completeness: 0 Status: draft Modified: 2018-03-19 17:15:35 +1100 Rubric: How about “Sufficient sufficiency,” is that taken?


I’m working through a small realisation, for my own interest, which has been helpful in my understanding of variational Bayes, - indeed in fact, relating it to non-Bayes variational inference. By starting from the idea of sufficient statistics, we come to the idea of variational inference in a natural way, via some other interesting points.

I doubt this insight is novel, but I will work through it as if it is, for the sake of my own education.

See also mixture models, probabilistic deep learning, directed graphical models, other probability metrics

To mention: Bayesian likelihood principle, Pitman–Koopman–Darmois theorem, connection to degrees of freedom.

Sufficient statistics in exponential families

Let’s start with sufficient statistics in exponential families, which, for reasons of historical pedagogy, are the Garden of Eden of Inference, the Garden of Edenference for short. I suspect that deep in their hearts, all statisticians regard themselves as prodigal exiles from of the exponential family, and long for the innocence of that Garden of Edenference.

Anyway, informally speaking, here’s what’s going on with the inference problems involving sufficient statistics. We are interested in estimating some parameter of inference, \(\theta\) using realisations \(x\) of some random process \(X\sim \mathbb{P}(x|\theta).\)

Then \(T(x)\) is a sufficient statistic for \(\theta\) iff

$$\mathbb{P}(x|T(x),\theta)= \mathbb{P}(x|T(x)).$$

That is, our inference about \(\theta\) depends on the data only through the sufficient statistic.

(mention size of sufficient statistic)

Fisher–Neyman factorization theorem

$$ \mathbb{P}(x;\theta)=h(x)g(T(x);\theta) $$

Famously, Maximum Likelihood estimators of exponential family models are highly compressible, in that these have sufficient statistics - these are low-dimensional functions of the data which characterise all the information in the complete data, with respect to the parameter estimates. Many models and data sets and estimation methods do not have this feature, even parametric models with very few parameters.

This can be a PITA when your data is very big and you wish to get benefit from that, and yet you can’t fit the data in memory; The question then arises - when can I do better? Can I find a “nearly sufficient” statistic, which is smaller than my data and yet does not worsen my error substantially? Can I quantify this nearness to the original?