Learn R Programming

STAR (version 0.3-7)

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

Description

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

Usage

data(ShallowShocks)

Arguments

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.

Source

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

Details

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

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

Run this code
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")

Run the code above in your browser using DataLab