Class of Fitted ETAS Models

Description

A class etas to represent a fitted ETAS model. The output of etas.

Details

An object of class etas represents an ETAS model that has been fitted to a spatio-temporal point pattern (catalog) of earthquakes. It is the output of the model fitter, etas.

The class etas has methods for the following standard generic functions:

Author(s)

Abdollah Jalilian jalilian@razi.ac.ir

See Also

etas,

Examples

  # fitting the ETAS model to an Iranian catalog

  data(iran.quakes)
  summary(iran.quakes)

  # preparing the catalog
  iran.cat <- catalog(iran.quakes, time.begin="1973/01/01",
     study.start="1985/01/01", study.end="2016/01/01",
     lat.range=c(25, 42), long.range=c(42, 63), mag.threshold=4.5)

  ## Not run: 
  plot(ir.cat)
  
## End(Not run)

  # initial parameters values
  param01 <- c(0.4339678,
   0.1988628,
   0.0345206,
   1.6290137,
   1.1286776,
   0.0072539,
   2.1705884,
   0.5706402)

  # fitting the model
  ## Not run: 
  res <- etas(ir.cat, param0=param01, no.itr=5)
  
## End(Not run)

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(ETAS)
Loading required package: maps

 # maps v3.1: updated 'world': all lakes moved to separate new #
 # 'lakes' database. Type '?world' or 'news(package="maps")'.  #


> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ETAS/etas.object.Rd_%03d_medium.png", width=480, height=480)
> ### Name: etas.object
> ### Title: Class of Fitted ETAS Models
> ### Aliases: etas.object methods.etas
> ### Keywords: spatial attribute
> 
> ### ** Examples
> 
>   # fitting the ETAS model to an Iranian catalog
> 
>   data(iran.quakes)
>   summary(iran.quakes)
         date               time           long            lat       
 2014-08-18:  44   00:04:50.70:   2   Min.   :40.01   Min.   :22.13  
 2011-10-23:  42   04:46:33.00:   2   1st Qu.:48.34   1st Qu.:28.32  
 2013-04-09:  31   07:02:27.00:   2   Median :51.88   Median :32.13  
 2013-05-11:  27   00:00:19.42:   1   Mean   :52.20   Mean   :32.80  
 2011-10-24:  19   00:00:41.44:   1   3rd Qu.:56.39   3rd Qu.:37.52  
 1990-06-21:  18   00:00:46.52:   1   Max.   :65.00   Max.   :41.98  
 (Other)   :5789   (Other)    :5961                                  
      mag       
 Min.   :4.000  
 1st Qu.:4.200  
 Median :4.400  
 Mean   :4.468  
 3rd Qu.:4.700  
 Max.   :6.200  
                
> 
>   # preparing the catalog
>   iran.cat <- catalog(iran.quakes, time.begin="1973/01/01",
+      study.start="1985/01/01", study.end="2016/01/01",
+      lat.range=c(25, 42), long.range=c(42, 63), mag.threshold=4.5)
> 
>   ## Not run: 
> ##D   plot(ir.cat)
> ##D   
> ## End(Not run)
> 
>   # initial parameters values
>   param01 <- c(0.4339678,
+    0.1988628,
+    0.0345206,
+    1.6290137,
+    1.1286776,
+    0.0072539,
+    2.1705884,
+    0.5706402)
> 
>   # fitting the model
>   ## Not run: 
> ##D   res <- etas(ir.cat, param0=param01, no.itr=5)
> ##D   
> ## End(Not run)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>