TY - JOUR
T1 - Bayesian inference with reliability methods without knowing the maximum of the likelihood function
AU - Betz, Wolfgang
AU - Beck, James L.
AU - Papaioannou, Iason
AU - Straub, Daniel
N1 - Publisher Copyright:
© 2018 Elsevier Ltd
PY - 2018/6
Y1 - 2018/6
N2 - In the BUS (Bayesian Updating with Structural reliability methods) approach, the uncertain parameter space is augmented by a uniform random variable and the Bayesian inference problem is interpreted as a structural reliability problem. A posterior sample is given by an augmented vector sample within the failure domain of the structural reliability problem where the realization of the uniform random variable is smaller than the likelihood function scaled by a constant c. The constant c must be selected such that 1∕c is larger or equal than the maximum of the likelihood function, which, however, is typically unknown a-priori. For BUS combined with sampling based reliability methods, choosing c too small has a negative impact on the computational efficiency. To overcome the problem of selecting c, we propose a post-processing step for BUS that returns an unbiased estimate for the evidence and samples from the posterior distribution, even if 1∕c is selected smaller than the maximum of the likelihood function. The applicability of the proposed post-processing step is demonstrated by means of rejection sampling. However, it can be combined with any structural reliability method applied within the BUS framework.
AB - In the BUS (Bayesian Updating with Structural reliability methods) approach, the uncertain parameter space is augmented by a uniform random variable and the Bayesian inference problem is interpreted as a structural reliability problem. A posterior sample is given by an augmented vector sample within the failure domain of the structural reliability problem where the realization of the uniform random variable is smaller than the likelihood function scaled by a constant c. The constant c must be selected such that 1∕c is larger or equal than the maximum of the likelihood function, which, however, is typically unknown a-priori. For BUS combined with sampling based reliability methods, choosing c too small has a negative impact on the computational efficiency. To overcome the problem of selecting c, we propose a post-processing step for BUS that returns an unbiased estimate for the evidence and samples from the posterior distribution, even if 1∕c is selected smaller than the maximum of the likelihood function. The applicability of the proposed post-processing step is demonstrated by means of rejection sampling. However, it can be combined with any structural reliability method applied within the BUS framework.
KW - Bayesian model class selection
KW - Bayesian updating
KW - Rejection sampling
KW - Structural reliability
UR - http://www.scopus.com/inward/record.url?scp=85045752635&partnerID=8YFLogxK
U2 - 10.1016/j.probengmech.2018.03.004
DO - 10.1016/j.probengmech.2018.03.004
M3 - Article
AN - SCOPUS:85045752635
SN - 0266-8920
VL - 53
SP - 14
EP - 22
JO - Probabilistic Engineering Mechanics
JF - Probabilistic Engineering Mechanics
ER -