model_posterior.distr: Posterior distribution estimation

In amignan/rseismTLS: Induced Seismicity Traffic Light System

Description

Estimates the posterior distributions of the 3 parameters of the induced seismicity model of Mignan et al. (2017), model_rate.val(), including joint posterior distributions.

Usage

model_posterior.distr(prior, joint_prior, LL)

Arguments

prior	the list of model parameter increments as defined in model_prior.distr: bi the vector of b increments ai the vector of a_fb increments taui the vector of tau increments
joint_prior	the joint prior computed from model_joint_prior.distr: joint.prior_norm the 3-dimensional array of the normalized joint prior distribution
LL	the log likelihood distribution computed from loglik_point.array

Details

Read Broccardo et al. (2017) for details.

Value

a list of the posterior distributions:

• a.post the vector of the a_fb posterior distribution

• b.post the vector of the b posterior distribution

• tau.post the vector of the tau posterior distribution

• b_a.post the array of the (b, a_fb) joint distribution

• b_tau.post the array of the (b, tau) joint distribution

• a_tau.post the array of the (a_fb, tau) joint distribution

• joint.post_norm the 3-dimensional array of the normalized joint posterior distribution

References

Broccardo M., Mignan A., Wiemer S., Stojadinovic B., Giardini D. (2017), Hierarchical Bayesian Modeling of Fluid‐Induced Seismicity. Geophysical Research Letters, 44 (22), 11,357-11,367, doi: 10.1002/2017GL075251

Mignan A., Broccardo M., Wiemer S., Giardini D. (2017), Induced seismicity closed-form traffic light system for actuarial decision-making during deep fluid injections. Sci. Rep., 7, 13607, doi: 10.1038/s41598-017-13585-9

See Also

model_prior.distr, model_joint_prior.distr, loglik_point.array

amignan/rseismTLS documentation built on July 4, 2020, 1:04 p.m.