Localisation

Adresses

Aix-Marseille Université
Institut de Mathématiques de Marseille (I2M) - UMR 7373
Site Saint-Charles : 3 place Victor Hugo, Case 19, 13331 Marseille Cedex 3
Site Luminy : Campus de Luminy - Case 907 - 13288 Marseille Cedex 9

Séminaire

A Monte Carlo EM for the Poisson Log-Normal Model




Date(s) : 20/11/2023   iCal
14h00 - 15h00

The Poisson log-normal (PLN) model is a generic model for the joint distribution of count data, accounting for covariates. It is also an incomplete data model. A classical way to achieve maximum likelihood inference for model parameters 𝜃 is to resort to the EM algorithm, which aims at maximising, with respect to 𝜃, the conditional expectation, given the observed data 𝑌, of the so-called complete log-likelihood. Unfortunately, the evaluation of the latter is intractable in the case of the PLN model because the conditional distribution of the latent vector conditionally on the corresponding observed count vector has no closed form and none of its moments can be evaluated in an efficient manner. Variational approaches have been studied to tackle this problem but lack from statistical guarantees. Indeed the resulting estimate does not enjoy the general properties of MLEs. In particular, its (asymptotic) variance is unknown, so no test nor confidence interval can be derived easily from the variational inference. Starting from already available variational approximations, we define a first Monte Carlo EM algorithm to obtain maximum likelihood estimators of this model. We then extend this first algorithm to the case of a composite likelihood in order to be able to handle higher dimensional count data. Both methods are statically grounded and provide confidence region for model parameters.

Emplacement
Saint-Charles - FRUMAM (2ème étage)

Catégories

Tags :

Leave a comment

Secured By miniOrange