We study the multivariate deconvolution problem of recovering the distribution of a signal from independent and identically distributed observations additively contaminated with random errors having known distribution. For errors with ordinary smooth distribution, we recast the multidimensional problem as a one-dimensional problem leveraging the equivalence between the L1-Wasserstein and the max-sliced L1-Wasserstein metrics and derive an inversion inequality relating the L1-Wasserstein distance between two distributions of the signal to the L1-distance between the corresponding mixture densities of the observations. This smoothing inequality outperforms existing inversion inequalities. We apply it to derive L1-Wasserstein rates of convergence for the distribution of the signal. As an application to the Bayesian framework, we consider L1-Wasserstein deconvolution with the Laplace noise in dimension one using a Dirichlet process mixture of normal densities as a prior measure for the mixing distribution. We construct an adaptive approximation of the sampling density by convolving the Laplace density with a well-chosen mixture of Gaussian densities and show that the posterior measure contracts at a nearly minimax-optimal rate, up to a log-factor, in the L1-distance. The rate automatically adapts to the unknown Sobolev regularity of the mixing density, thus leading to a new Bayesian adaptive estimation procedure over the full scale of regularity levels. We illustrate the utility of the inversion inequality also in a frequentist setting by showing that a minimum distance estimator attains the minimax convergence rates for L1-Wasserstein deconvolution in any dimension d≥1, lower bounds being derived here.

Wasserstein convergence in Bayesian and frequentist deconvolution models

Catia Scricciolo
2024-01-01

Abstract

We study the multivariate deconvolution problem of recovering the distribution of a signal from independent and identically distributed observations additively contaminated with random errors having known distribution. For errors with ordinary smooth distribution, we recast the multidimensional problem as a one-dimensional problem leveraging the equivalence between the L1-Wasserstein and the max-sliced L1-Wasserstein metrics and derive an inversion inequality relating the L1-Wasserstein distance between two distributions of the signal to the L1-distance between the corresponding mixture densities of the observations. This smoothing inequality outperforms existing inversion inequalities. We apply it to derive L1-Wasserstein rates of convergence for the distribution of the signal. As an application to the Bayesian framework, we consider L1-Wasserstein deconvolution with the Laplace noise in dimension one using a Dirichlet process mixture of normal densities as a prior measure for the mixing distribution. We construct an adaptive approximation of the sampling density by convolving the Laplace density with a well-chosen mixture of Gaussian densities and show that the posterior measure contracts at a nearly minimax-optimal rate, up to a log-factor, in the L1-distance. The rate automatically adapts to the unknown Sobolev regularity of the mixing density, thus leading to a new Bayesian adaptive estimation procedure over the full scale of regularity levels. We illustrate the utility of the inversion inequality also in a frequentist setting by showing that a minimum distance estimator attains the minimax convergence rates for L1-Wasserstein deconvolution in any dimension d≥1, lower bounds being derived here.
2024
Inversion inequalities
Max-sliced Wasserstein metrics
Minimax rates
Mixtures of Laplace densities
Multivariate deconvolution
Rates of convergence
Sobolev classes
Wasserstein metrics
Adaptation
Density estimation
Dirichlet process mixtures
File in questo prodotto:
File Dimensione Formato  
revision_final_main - ACCEPTED.pdf

embargo fino al 31/12/2025

Descrizione: Main Paper
Tipologia: Documento in Post-print
Licenza: Creative commons
Dimensione 538.09 kB
Formato Adobe PDF
538.09 kB Adobe PDF   Visualizza/Apri   Richiedi una copia
revision_final_supp.pdf

embargo fino al 31/12/2025

Descrizione: Supplementary Material
Tipologia: Documento in Post-print
Licenza: Creative commons
Dimensione 611.81 kB
Formato Adobe PDF
611.81 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11562/1110507
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 0
social impact