italy.quakes | R Documentation |
A data frame with 2158 rows and 6 columns giving occurrence date, time, longitude, latitude, magnitude and depth of earthquakes occurred since 2005-04-16 till 2013-11-01 in Italy and its vicinity (6.15-19E and 35-48N) with magnitude greater than or equal to 3.
data(italy.quakes)
An object of class "data.frame"
containing the following columns:
date
Occurrence date in the format "yyyy/mm/dd"
time
Occurrence time in decimal days after the first earthquake
long
Latitude of epicenter in decimal degrees
lat
Latitude of epicenter in decimal degrees
mag
Magnitude of each earthquake by JMA (Japan Meteorological Agency)
depth
Depth of each earthquake
Data are retrieved from the Italian Seismological Instrumental and Parametric Data Base (ISIDE) available at http://iside.rm.ingv.it by Marcello Chiodi and Giada Adelfio in the etasFLP package.
Marcello Chiodi and Giada Adelfio (2015). etasFLP: Mixed FLP and ML Estimation of ETAS Space-Time Point Processes. R package version 1.3.0. https://CRAN.R-project.org/package=etasFLP.
# creat an earthquake catalog
italy.cat <- catalog(italy.quakes, dist.unit="km")
## Not run:
plot(italy.cat)
## End(Not run)
# set initial parameter values
mu <- 1
k0 <- 0.005
c <- 0.005
alpha <- 1.05
p <- 1.01
D <- 1.1
q <- 1.52
gamma <- 0.6
# reparametrization: transform k0 to A
A <- pi * k0 / ((p - 1) * c^(p - 1) * (q - 1) * D^(q - 1))
param0 <- c(mu, A, c, alpha, p, D, q, gamma)
## Not run:
italy.fit <- etas(italy.cat, param0)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.