On MCMC algorithm for Subset Simulation

Au, Siu-Kui ORCID: 0000-0002-0228-1796 (2016) On MCMC algorithm for Subset Simulation. PROBABILISTIC ENGINEERING MECHANICS, 43. pp. 117-120.

Text pem - short comm.pdf - Author Accepted Manuscript Download (407kB)

Abstract

A new Markov Chain Monte Carlo (MCMC) algorithm for Subset Simulation was recently proposed by imposing a joint Gaussian distribution between the current sample and the candidate. It coincides with the limiting case of the original independent-component algorithm where each random variable is represented by an infinite number of hidden variables. The algorithm is remarkably simple as it no longer involves the explicit choice of proposal distribution. It opens up a new perspective for generating conditional failure samples and potentially allows more direct and flexible control of algorithm through the cross correlation matrix between the current sample and the candidate. While by definition the cross correlation matrix need not be symmetric, this article shows that it must be so in order to satisfy detailed balance and hence to produce an unbiased algorithm. The effect of violating symmetry on the distribution of samples is discussed and insights on acceptance probability are provided.

Item Type:	Article
Uncontrolled Keywords:	Detailed balance, Rare event, Markov Chain Monte Carlo, Monte Carlo, Subset Simulation
Depositing User:	Symplectic Admin
Date Deposited:	12 Sep 2016 08:38
Last Modified:	19 Jan 2023 07:30
DOI:	10.1016/j.probengmech.2015.12.003
Related URLs:	Author Publisher
URI:	https://livrepository.liverpool.ac.uk/id/eprint/3003221