negloglik_hist.val: Negative Log Likelihood Function (Histogram Data)

In amignan/rseismTLS: Induced Seismicity Traffic Light System

Description

Estimates the negative log likelihood of the Non-Homogeneous Poisson Process (NHPP) derived from the statistical model of Mignan et al. (2017) for histogram data and a specific set of input parameters, for three possible time windows (full sequence, injection phase or post-injection phase).

Usage

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

Arguments

data	a list containing all the necessary histogram data: seism.binned an earthquake-rate histogram data frame of parameters t (time in days) and rate (per time bin) inj.binned matching binned injection profile data frame of parameters t (time in days), dV (flow rate in cubic metre/time bin) and V (cumulative injected volume in cubic metres) b the b-value of the earthquake catalogue (not required for "post-injection") m0 the minimum magnitude cutoff of the earthquake catalogue (not required for "post-injection") ts shut-in time (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
window	the window to be used: "injection", "post-injection", or "full sequence"

Details

Note that the function is here structured for the parameters to be optimized by optim()``, which is done automatically in model_par.mle_hist(). The parvector depends on the selectedwindow, which is dealt with in model_par.mle_hist()' (see Details section there).

Value

the negative log likelihood estimate for the specified par values

References

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_hist, model_par.mle_point, negloglik_point.val, data.bin

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