Saved in:
Bibliographic Details
Main Authors: Liao, Aaron R., Oguri, Kenshiro, Carpenter, Michele D.
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2604.09375
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866915930243596288
author Liao, Aaron R.
Oguri, Kenshiro
Carpenter, Michele D.
author_facet Liao, Aaron R.
Oguri, Kenshiro
Carpenter, Michele D.
contents Accurate representation of non-Gaussian distributions of quantities of interest in nonlinear dynamical systems is critical for estimation, control, and decision-making, but can be challenging when forward propagations are expensive to carry out. This paper presents an approach for estimating probability density functions of states evolving under nonlinear dynamics using Seminonparametric (SNP), or Gallant-Nychka, densities. SNP densities employ a probabilists' Hermite polynomial basis to model non-Gaussian behavior and are positive everywhere on the support by construction. We use Monte Carlo to approximate the expectation integrals that arise in the maximum likelihood estimation of SNP coefficients, and introduce a convex relaxation to generate effective initial estimates. The method is demonstrated on density and quantile estimation for the chaotic Lorenz system. The results demonstrate that the proposed method can accurately capture non-Gaussian density structure and compute quantiles using significantly fewer samples than raw Monte Carlo sampling.
format Preprint
id arxiv_https___arxiv_org_abs_2604_09375
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Data-Efficient Non-Gaussian Semi-Nonparametric Density Estimation for Nonlinear Dynamical Systems
Liao, Aaron R.
Oguri, Kenshiro
Carpenter, Michele D.
Optimization and Control
Statistics Theory
Machine Learning
Accurate representation of non-Gaussian distributions of quantities of interest in nonlinear dynamical systems is critical for estimation, control, and decision-making, but can be challenging when forward propagations are expensive to carry out. This paper presents an approach for estimating probability density functions of states evolving under nonlinear dynamics using Seminonparametric (SNP), or Gallant-Nychka, densities. SNP densities employ a probabilists' Hermite polynomial basis to model non-Gaussian behavior and are positive everywhere on the support by construction. We use Monte Carlo to approximate the expectation integrals that arise in the maximum likelihood estimation of SNP coefficients, and introduce a convex relaxation to generate effective initial estimates. The method is demonstrated on density and quantile estimation for the chaotic Lorenz system. The results demonstrate that the proposed method can accurately capture non-Gaussian density structure and compute quantiles using significantly fewer samples than raw Monte Carlo sampling.
title Data-Efficient Non-Gaussian Semi-Nonparametric Density Estimation for Nonlinear Dynamical Systems
topic Optimization and Control
Statistics Theory
Machine Learning
url https://arxiv.org/abs/2604.09375