Tossing half-coins and other partial coins: signed probabilities and Sibuya distribution
Nikolai Leonenko, Igor Podlubny
TL;DR
This work addresses the challenge of simulating signed probability distributions that arise in partial coins and fractional-order processes by leveraging the Sibuya distribution and a generating-function decomposition. The authors provide explicit constructions for nonnegative distributions $g$ and $h$ such that $f g = h$ for the $1/n$-coin, establishing positivity of the corresponding coefficients and enabling Monte Carlo simulation via cumulated-mass methods. They develop a concrete numerical framework, derive recurrence-based coefficient evaluations, and implement a MATLAB toolbox to perform partial-coin and multi-coin simulations, including biased cases. The results extend signed-distribution simulation to infinite-support cases and have potential implications for fractional differentiation, uncertainty quantification, and related stochastic modeling tasks.
Abstract
A method for the numerical simulation of signed probability distributions for the case of tossing $1/n$-th of a coin is presented and illustrated by examples.
