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")