gutenberg {geostats}

create a Gutenberg-Richter plot

Description

Calculate a semi-log plot with earthquake magnitude on the horizontal axis,and the cumulative number of earthquakes exceeding any given magnitude on the vertical axis.

Usage

gutenberg(m, n = 10, ...)

Arguments

m

a vector of earthquake magnitudes

n

the number of magnitudes to evaluate

...

optional arguments to the generic points function.

Value

the output of lm with earthquake magnitude as the independent variable (mag) and the logarithm (base 10) of the frequency as the dependent variable (lfreq).

Examples

data(declustered,package='geostats')
gutenberg(declustered$mag)