analysis.APositiveAValueEstimator#

class seismostats.analysis.APositiveAValueEstimator#

Returns the a-value of the Gutenberg-Richter (GR) law using only the earthquakes with magnitudes \(m_i \ge m_{i-1} + dmc\).

Source:: Following the idea of positivity of van der Elst 2021 (J Geophysical Research: Solid Earth, Vol 126, Issue 2). Note: This is not a-positive as defined by van der Elst and Page 2023 (JGR: Solid Earth, Vol 128, Issue 10).

Examples

>>> import numpy as np
>>> from seismostats.analysis import APositiveAValueEstimator

>>> magnitudes = np.array([2.1, 2.3, 2.0, 2.0, 2.1, 2.2, 2.1, 2.3,
...                        2.0, 2.0])
>>> times = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9, 10])

>>> my_estimator = APositiveAValueEstimator()
>>> my_estimator.calculate(
...     magnitudes=magnitudes, mc=2.0, delta_m=0.1, times=times)

>>> my_estimator.a_value

0.9542425094393249

Attributes

`a_value`	The a value of the Gutenberg-Richter law.
`b_value`	The b-value of the Gutenberg-Richter law.
`delta_m`	The bin size of the discretized magnitudes.
`m_ref`	The reference magnitude.
`magnitudes`	The magnitudes used to estimate the a-value.
`mc`	The completeness magnitude.
`n`	Number of magnitudes used to estimate the a-value.
`scaling_factor`	The scaling factor.
`std`	Uncertainty of the a-value estimate, implemented as a placeholder.
`value`	The a value of the Gutenberg-Richter law.

Methods

calculate

param magnitudes:: Array of magnitudes.