Thursday, July 29, 2010

Kolia Sadeghi : August 4

again

I will attempt to take the latest paper on the deterministic particle flow filter discussed in a previous blog post, and strip it down to the essentials.  The authors present a more general, stable and improved version of their previous deterministic particle flow filter, supposedly. This paper is rife with ideas and peculiarly written; for a gentler introduction, please refer to the papers linked to in the previous blog post. Here is the paper:

Exact particle flow for nonlinear filters by Fred Daum, Jim Huang and Arjang Noushin,
Numerical experiments for nonlinear filters with exact particle flow induced by log-homotopy (companion paper)

Friday, July 23, 2010

Deterministic particle filtering

No resampling, rejection, or importance sampling are used. Particles are propagated through time by numerically integrating an ODE. The method is very similar in spirit to Jascha Sohl-Dickstein, Peter Battaglino and Mike DeWeese's Minimum probability flow learning, but applied to nonlinear filtering.


The authors report orders of magnitude speedups for higher dimensional state spaces where sampling rejection would be a problem.

Particle flow for nonlinear filters with log-homotopy by Fred Daum & Jim Huang


There are a couple of papers companion to this one:
Nonlinear filters with particle flow induced by log-homotopy
Seventeen dubious methods to approximate the gradient for nonlinear filters with particle flow


As you may see, the authors have a very peculiar writing style.


However, one very recent paper by Lingji Chen and Raman Mehra points out some flaws in the approach:
A study of nonlinear filters with particle flow induced by log-homotopy
(but see the group meeting announcement above for Fred Daum and Jim Huang's recent answer to this).

Kolia Sadeghi : July 28

I will present work done with Liam, Jeff Gauthier and others in EJ Chichilnisky's lab on locating retinal cones from multiple ganglion cell recordings.  We write down a single hierarchical model where ganglion cell responses are modeled as independent GLMs with space-time-color separable filters and no spike history.  Assuming the stimulus was gaussian ensures that the ganglion cell Spike Triggered Averages are sufficient statistics.  The spatial component is then assumed to be a weighted sum of non-overlapping and appropriately placed archetypical cone receptive fields.  With a benign approximation, we can integrate out the weights and focus on doing MCMC in the space of cone locations and colors only.  As it turns out, this likelihood landscape has many nasty local maxima; we use parallel tempering and a few techniques specific to this problem to ensure ergodicity of the markov chain.

Doing a google scholar search on parallel tempering, also known as replica exchange, or just exchange Monte Carlo, will bring up many papers on this simple technique. Here is a review:
Parallel tempering: Theory, applications, and new perspectives

Thursday, July 22, 2010

Some classic stats papers

These are a bit more old-school, but still of interest:

Some interesting papers from AISTATS 2010

Here are a few potentially interesting papers from AISTATS this year. All pdf's available from

by Botond Cseke, Tom Heskes

by Lauren Hannah, David Blei, Warren Powell

by Jun Li, Dacheng Tao

by Mark Schmidt, Kevin Murphy

by Sajid Siddiqi, Byron Boots, Geoffrey Gordon

by Aarti Singh, Robert Nowak, Robert Calderbank

by Nikolai Slavov

by Bharath Sriperumbudur, Kenji Fukumizu, Gert Lanckriet

by Ryan Turner, Marc Deisenroth, Carl Rasmussen

by James Martens, Ilya Sutskever

by Jimmy Olsson, Jonas Strojby

by Steve Hanneke, Liu Yang

Jittering spike trains carefully


Presented in lab meeting by Alex Ramirez on July 14th, 2010.

by Matthew Harrison and Stuart Geman

In it the authors describe an algorithm that takes a spike train and jitters the spike times to create a new spike train which is maximally random while preserving the firing rate and recent spike-history of the original train.