Abstract:Flexible continuous univariate distributions with bounded support are essential for accurate input modeling in stochastic simulation and decision analysis. Although Bézier distributions provide a powerful family capable of representing complex shapes, their adoption has been hindered by the lack of efficient fitting procedures and modern software implementations. This paper develops a computational framework for fitting Bézier distributions to empirical data via both minimum error and maximum likelihood estimation, leveraging first-order optimization methods and exploiting the geometry of the parameter space. We identify provably (asymptotically) lossless convex restrictions of the feasible set that enable efficient projection operators based on isotonic regression and develop first-order algorithms that reduce computational runtime by three to four orders of magnitude compared to traditional derivative-free methods, while delivering consistent fits across real-world data. When benchmarked against the nonlinear solver IPOPT, our methods prove three orders of magnitude faster on average and more robust, while achieving comparable accuracy. To bridge the gap between theory and practice, we introduce bezierv, an open-source Python package providing a unified interface for fitting, analyzing, and convolving Bézier distributions.
From: Esteban Leiva [view email]
[v1]
Mon, 22 Jun 2026 18:23:31 UTC (143 KB)
此内容由惯性聚合(RSS阅读器)自动聚合整理，仅供阅读参考。 原文来自 — 版权归原作者所有。