Guillaume Chennetier
banner
gchennetier.bsky.social
Guillaume Chennetier
@gchennetier.bsky.social
PhD Student in probability and statistics (Polytechnique x EDF)

Working on rare event simulation for piecewise deterministic Markov processes

#probability #statistics
The following paper by Thomas Galtier et al. gives the probability density function (pdf) of a PDMP trajectory of fixed duration and its associated reference measure : www.esaim-ps.org/articles/ps/... 12/n
On the optimal importance process for piecewise deterministic Markov process | ESAIM: Probability an...
PS : ESAIM: Probability and Statistics, publishes original research and survey papers in the area of Probability and Statistics
www.esaim-ps.org
January 8, 2024 at 2:24 PM
The idea is to generate the PDMP trajectories according to an alternative (importance) distribution that realizes the event of interest more frequently, and to rectify the bias of the estimator by introducing a suited likelihood ratio. 11/n
January 8, 2024 at 2:23 PM
A standard Monte Carlo method (generating a large number of PDMP trajectories modeling the system, and returning the proportion of simulations realizing the event of interest) is inefficient when the event is rare. We therefore turn to an importance sampling method. 10/n
January 8, 2024 at 2:23 PM
Critical failure occurs when physical variables exceed a given threshold, and is only possible when some key combinations of components fail at the same time and are not repaired in time. It is therefore a rare event (of the order of 10^-5 at most in our test cases). 9/n
January 8, 2024 at 2:23 PM
The distribution of waiting times before breakdown or repair depends on the state of the components and continuously on the physical variables. We are interested in the probability that the system will go into critical failure before a fixed date. 8/n
January 8, 2024 at 2:23 PM
The system's continuous variables (temperature, pressure, liquid levels, etc.) evolve according to deterministic physical laws that depend on the state of the various system components (pumps, valves, generators, etc.), which are subject to random breakdowns and repairs. 7/n
January 8, 2024 at 2:22 PM
PDMPs model a wide range of phenomena, and are particularly popular in biology. They also form the basis of PD-MCMC methods in Bayesian inference for simulating high-dimensional distributions. We use them to model complex hybrid dynamic industrial systems. 6/n
January 8, 2024 at 2:22 PM
A PDMP is a "hybrid" continuous-time stochastic process. It jumps at random times to new states like a classical non-diffusive Markovian process, but follows deterministic (non-trivial) trajectories between jumps. 5/n
January 8, 2024 at 2:22 PM
I give some details below on PDMPs, rare event simulation, importance sampling, the committor function, the cross entropy method, recycling schemes and our family of approximations of the committor function. 4/n
January 8, 2024 at 2:22 PM
Our contribution: a new family of importance distributions based on committor function approximation for PDMPs modeling complex high-dimensional industrial systems, and a central limit theorem for the adaptive importance sampling procedure of PDMPs with recycling scheme. 3/n
January 8, 2024 at 9:29 AM
Our goal: to estimate the (very small) probability that a PDMP (piecewise deterministic Markov process) trajectory will reach a critical region of its state space before a fixed date, using an importance sampling method. 2/n
January 8, 2024 at 9:28 AM
Some keywords not in the title of the paper: rare event simulation, reliability assessment, cross entropy minimization with recycling scheme, Monte Carlo, stochastic processes, committor function approximation. 1/n
January 8, 2024 at 9:22 AM