Being Monte Carlo
methods, particle filters have over Laplace approximations the advantage of much greater
flexibility. This comes however at the expense of a higher representational and computational
complexity. To combine the flexibility of particle filtering with the computational
advantage of parametric methods, we propose a variational approximation