Penalized likelihood estimation of probability density functions using compositional splines

2026-08-24

AI Summary

本文提出了一种基于中心对数比变换和ZB样条基函数的惩罚最大似然框架，直接从原始观测数据估计概率密度函数，减少两阶段方法带来的偏差。

Abstract

Probability density functions are commonly estimated through preliminary smoothing or aggregation procedures, e.g., histograms or kernel density estimation, before subsequent functional representation and functional data analyses. Such a two-stage approach can lead to additional approximation bias and weaken the direct connection between the observed data and the underlying distributional structure. In this paper, we propose a penalized maximum likelihood framework for direct estimation of probability density functions from raw observations within the framework of Bayes Hilbert spaces while preserving the compositional geometry of densities. The methodology is based on the centered log-ratio (clr) transformation, an isometric isomorphism between Bayes Hilbert spaces and the standard Lebesgue space of square integrable functions with zero integral, enabling efficient spline representations. The clr transformed densities are represented using ZB-spline basis functions, while their smoothness is controlled through quadratic penalties imposed on the spline coefficients. The proposed framework is developed for univariate and bivariate densities. In the latter case, it naturally incorporates the orthogonal decomposition into independent and interactive parts together with the corresponding geometric marginals. The performance is evaluated in a simulation study involving multiple complex scenarios and compared with kernel smoothing. Finally, the applicability of the framework is illustrated using empirical geochemical data.

Problem

• probability density functions
• approximation bias
• Bayes Hilbert spaces
• compositional geometry
• centered log-ratio

Innovation

• penalized maximum likelihood
• Bayes Hilbert spaces
• centered log-ratio (clr) transformation
• ZB-spline basis functions
• quadratic penalties