| ShallowShocks {STAR} | R Documentation |
Earthquakes data used by Yosihiko Ogata in his 1988 JASA paper.
data(ShallowShocks)
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. |
Quakes 213 and 214 were given exactly the same dates in Ogata (1988). Quake 214 has here been delayed by 1 minute.
Ogata (1988) Table 1, pp 14-15.
Ogata, Yosihiko (1988) Statistical Models for Earthquake Occurrences and Residual Analysis for Point Processes. Journal of the American Statistical Association 83: 9-27.
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")