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)

Format

A data.frame with the following variables:

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

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(STAR)
Loading required package: survival
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-12. For overview type 'help("mgcv-package")'.
Loading required package: R2HTML
Loading required package: gss
Loading required package: codetools
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/STAR/ShallowShocks.Rd_%03d_medium.png", width=480, height=480)
> ### Name: ShallowShocks
> ### Title: Shallow Shocks (M >= 6.0) in OFF Tohoku Area for 1885-1980
> ### Aliases: ShallowShocks
> ### Keywords: datasets
> 
> ### ** 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")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>