negloglik_point.val: Negative Log Likelihood Function (Point Data)

In amignan/rseismTLS: Induced Seismicity Traffic Light System

Description

Estimates the negative log likelihood of the statistical model of Mignan et al. (2017) for a specific set of input parameters, for three possible time windows (full sequence, injection phase or post-injection phase).

Usage

negloglik_point.val(data, par, window = "full sequence")

Arguments

data	a list containing all the necessary point data: seism an earthquake catalogue data frame of parameters t (time in days) and m (magnitude) inj matching injection profile data frame of parameters t (time in days), dV (flow rate in cubic metre/day) and V (cumulative injected volume in cubic metres) m0 the minimum magnitude cutoff of the earthquake catalogue ts shut-in time (in days) (not required for "injection") Tmax upper range of the time window (in days) (not required for "injection") lambda0 seismicity rate (per day) at shut-in (only required for "post-injection")
par	a vector of the input parameters a_fb, tau, and b (in this order)
window	the window to be used: "injection", "post-injection", or "full sequence"

Details

See Eq. A2 of Broccardo et al. (2017) for 'full sequence' and Eq. A3 for 'injection'. Note that the function is here structured for the parameters to be optimized by optim()``, which is done automatically in model_par.mle_point(). The parvector depends on the selectedwindow, which is dealt with in model_par.mle_point()' (see Details section there).

Value

the negative log likelihood estimate for the specified par values

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_par.mle_point, model_par.mle_hist, negloglik_hist.val, loglik_point.array

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