Functional Estimation of the Marginal Likelihood

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Papaspiliopoulos, Omiros, Stumpf-Fétizon, Timothée, Weare, Jonathan
Format: Preprint
Publié: 2026
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866910014555291648
author Papaspiliopoulos, Omiros
Stumpf-Fétizon, Timothée
Weare, Jonathan
author_facet Papaspiliopoulos, Omiros
Stumpf-Fétizon, Timothée
Weare, Jonathan
contents We propose a framework for computing, optimizing and integrating with respect to a smooth marginal likelihood in statistical models that involve high-dimensional parameters/latent variables and continuous low-dimensional hyperparameters. The method requires samples from the posterior distribution of the parameters for different values of the hyperparameters on a simulation grid and returns inference on the marginal likelihood defined everywhere on its domain, and on its functionals. We show how the method relates to many of the methods that have been used in this context, including sequential Monte Carlo, Gibbs sampling, Monte Carlo maximum likelihood, and umbrella sampling. We establish the consistency of the proposed estimators as the sampling effort increases, both when the simulation grid is kept fixed and when it becomes dense in the domain. We showcase the approach on Gaussian process regression and classification and crossed effect models.
format Preprint
id arxiv_https___arxiv_org_abs_2602_07148
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Functional Estimation of the Marginal Likelihood
Papaspiliopoulos, Omiros
Stumpf-Fétizon, Timothée
Weare, Jonathan
Methodology
Computation
62-08, 65C05, 62F15
G.3
We propose a framework for computing, optimizing and integrating with respect to a smooth marginal likelihood in statistical models that involve high-dimensional parameters/latent variables and continuous low-dimensional hyperparameters. The method requires samples from the posterior distribution of the parameters for different values of the hyperparameters on a simulation grid and returns inference on the marginal likelihood defined everywhere on its domain, and on its functionals. We show how the method relates to many of the methods that have been used in this context, including sequential Monte Carlo, Gibbs sampling, Monte Carlo maximum likelihood, and umbrella sampling. We establish the consistency of the proposed estimators as the sampling effort increases, both when the simulation grid is kept fixed and when it becomes dense in the domain. We showcase the approach on Gaussian process regression and classification and crossed effect models.
title Functional Estimation of the Marginal Likelihood
topic Methodology
Computation
62-08, 65C05, 62F15
G.3
url https://arxiv.org/abs/2602.07148