ShallowShocks: Shallow Shocks (M >= 6.0) in OFF Tohoku Area for 1885-1980

In STAR: Spike Train Analysis with R

Description

Earthquakes data used by Yosihiko Ogata in his 1988 JASA paper.

Usage

data(ShallowShocks)

Format

A data.frame with the following variables:

year:	year of occurrence.
month:	month of occurrence.
day:	day of occurrence.
hour:	hour of occurrence.
minute:	minute of occurrence.
magnitude:	magnitude on Richter's scale.
type:	type of earthquake: main (shock), foreshock, aftershock; according to Utsu.
Date:	date in days starting from January 1st 1885.
energy.sqrt:	square root of the energy expressed in erg.

Details

Quakes 213 and 214 were given exactly the same dates in Ogata (1988). Quake 214 has here been delayed by 1 minute.

Source

Ogata (1988) Table 1, pp 14-15.

References

Ogata, Yosihiko (1988) Statistical Models for Earthquake Occurrences and Residual Analysis for Point Processes. Journal of the American Statistical Association 83: 9-27.

Examples

data(ShallowShocks)
## Reproduce Fig. 2 of Ogata 1988
layout(matrix(1:3, nrow = 3))
plot(ShallowShocks$Date,
     cumsum(ShallowShocks$energy.sqrt) / 10^13,
     type ="l",
     xlab = "",
     ylab = "",
     main = "Cumulative square root of energy")
plot(ShallowShocks$Date,
     cumsum(1+numeric(dim(ShallowShocks)[1])),
     type ="l",
     xlab = "",
     ylab = "",
     main = "Cumulative number of shocks")
plot(ShallowShocks$Date,
     ShallowShocks$magnitude,
     type = "h",
     ylim = c(5,9),
     xlab = "Time (days)",
     ylab = "",
     main = "Magnitude vs Occurrence time")

Example output

Loading required package: survival
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-20. For overview type 'help("mgcv-package")'.
Loading required package: R2HTML
Loading required package: gss
Loading required package: codetools

STAR documentation built on May 2, 2019, 11:44 a.m.