Variational message passing (VMP) is an approximate inference technique for continuous- or discrete-valued Bayesian networks, with conjugate-exponential parents, developed by John Winn. VMP was developed as a means of generalizing the approximate variational methods used by such techniques as Latent Dirichlet allocation and works by updating an approximate distribution at each node through messages in the node's Markov blanket.
Contents
Likelihood Lower Bound
Given some set of hidden variables
So, if we define our lower bound to be
then the likelihood is simply this bound plus the relative entropy between
where
Determining the Update Rule
The likelihood estimate needs to be as large as possible; because it's a lower bound, getting closer
where
Messages in Variational Message Passing
Parents send their children the expectation of their sufficient statistic while children send their parents their natural parameter, which also requires messages to be sent from the co-parents of the node.
Relationship to Exponential Families
Because all nodes in VMP come from exponential families and all parents of nodes are conjugate to their children nodes, the expectation of the sufficient statistic can be computed from the normalization factor.
VMP Algorithm
The algorithm begins by computing the expected value of the sufficient statistics for that vector. Then, until the likelihood converges to a stable value (this is usually accomplished by setting a small threshold value and running the algorithm until it increases by less than that threshold value), do the following at each node:
- Get all messages from parents
- Get all messages from children (this might require the children to get messages from the co-parents)
- Compute the expected value of the nodes sufficient statistics
Constraints
Because every child must be conjugate to its parent, this limits the types of distributions that can be used in the model. For example, the parents of a Gaussian distribution must be a Gaussian distribution (corresponding to the Mean) and a gamma distribution (corresponding to the precision, or one over