Learning Lévy density via adaptive RKHS regression with bi-level optimization

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Yang, Luxuan, Lu, Fei, Gao, Ting, Wei, Wei, Duan, Jinqiao
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866908737744142336
author Yang, Luxuan
Lu, Fei
Gao, Ting
Wei, Wei
Duan, Jinqiao
author_facet Yang, Luxuan
Lu, Fei
Gao, Ting
Wei, Wei
Duan, Jinqiao
contents We propose a nonparametric method to learn the Lévy density from probability density data governed by a nonlocal Fokker-Planck equation. We recast the problem as identifying the kernel in a nonlocal integral operator from discrete data, which leads to an ill-posed inverse problem. To regularize it, we construct an adaptive reproducing kernel Hilbert space (RKHS) whose kernel is built directly from the data. Under standard source and spectral decay conditions, we show that the reconstruction error decays in the mesh size at a near optimal rate. Importantly, we develop a generalized singular value decomposition (GSVD)-based bilevel optimization algorithm to choose the regularization parameter, leading to efficient and robust computation of the regularized estimator. Numerical experiments for several Lévy densities, drift fields and data types (PDE-based densities and sample ensemble-based KDE reconstructions) demonstrate that our bilevel RKHS method outperforms classical L-curve and generalized cross-validation strategies and that the adaptive RKHS norm is more accurate and robust than $L^2_ρ$- and $\ell^2$-based regularization.
format Preprint
id arxiv_https___arxiv_org_abs_2512_23621
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Learning Lévy density via adaptive RKHS regression with bi-level optimization
Yang, Luxuan
Lu, Fei
Gao, Ting
Wei, Wei
Duan, Jinqiao
Numerical Analysis
We propose a nonparametric method to learn the Lévy density from probability density data governed by a nonlocal Fokker-Planck equation. We recast the problem as identifying the kernel in a nonlocal integral operator from discrete data, which leads to an ill-posed inverse problem. To regularize it, we construct an adaptive reproducing kernel Hilbert space (RKHS) whose kernel is built directly from the data. Under standard source and spectral decay conditions, we show that the reconstruction error decays in the mesh size at a near optimal rate. Importantly, we develop a generalized singular value decomposition (GSVD)-based bilevel optimization algorithm to choose the regularization parameter, leading to efficient and robust computation of the regularized estimator. Numerical experiments for several Lévy densities, drift fields and data types (PDE-based densities and sample ensemble-based KDE reconstructions) demonstrate that our bilevel RKHS method outperforms classical L-curve and generalized cross-validation strategies and that the adaptive RKHS norm is more accurate and robust than $L^2_ρ$- and $\ell^2$-based regularization.
title Learning Lévy density via adaptive RKHS regression with bi-level optimization
topic Numerical Analysis
url https://arxiv.org/abs/2512.23621