Log In Sign Up

Predictive risk estimation for the Expectation Maximization algorithm with Poisson data

by   Paolo Massa, et al.

In this work, we introduce a novel estimator of the predictive risk with Poisson data, when the loss function is the Kullback-Leibler divergence, in order to define a regularization parameter's choice rule for the Expectation Maximization (EM) algorithm. To this aim, we prove a Poisson counterpart of the Stein's Lemma for Gaussian variables, and from this result we derive the proposed estimator showing its analogies with the well-known Stein's Unbiased Risk Estimator valid for a quadratic loss. We prove that the proposed estimator is asymptotically unbiased with increasing number of measured counts, under certain mild conditions on the regularization method. We show that these conditions are satisfied by the EM algorithm and then we apply this estimator to select its optimal reconstruction. We present some numerical tests in the case of image deconvolution, comparing the performances of the proposed estimator with other methods available in the literature, both in the inverse crime and non-inverse crime setting.


page 11

page 12

page 15


An Expectation Maximization Framework for Preferential Attachment Models

In this paper we develop an Expectation Maximization(EM) algorithm to es...

Coupled Bootstrap Test Error Estimation for Poisson Variables

Test error estimation is a fundamental problem in statistics and machine...

A Stochastic Path-Integrated Differential EstimatoR Expectation Maximization Algorithm

The Expectation Maximization (EM) algorithm is of key importance for inf...

Reconstruct Kaplan–Meier Estimator as M-estimator and Its Confidence Band

The Kaplan–Meier (KM) estimator, which provides a nonparametric estimate...

Low-rank matrix denoising for count data using unbiased Kullback-Leibler risk estimation

This paper is concerned by the analysis of observations organized in a m...

Monotonic Alpha-divergence Minimisation

In this paper, we introduce a novel iterative algorithm which carries ou...