R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Clustering Part of Conditional Intensity Function of the ETAS...

lambda

R Documentation

Clustering Part of Conditional Intensity Function of the ETAS Model

Description

A function to compute the clustering part of the conditional intensity function of the ETAS model at specified time and location.

Usage

  lambda(t, x, y, param, object)

Arguments

`t`	A numeric value. The time that the conditional intensity is to be computed at.
`x`	A numeric value. The x-coordinate of the location that the conditional intensity is to be computed at.
`y`	A numeric value. The y-coordinate of the location that the conditional intensity is to be computed at.
`param`	Vector of model paramters.
`object`	An object of class `"catalog"` containing an earthquake catalog dataset.

Details

For a given t, x and y, this function computes

sum[t[i] < t] k(m[i]) g(t - t[i]) f(x - x[i], y - y[i]|m[i]).

Value

A numeric value.

Author(s)

Abdollah Jalilian jalilian@razi.ac.ir

References

Zhuang, J., Ogata, Y. and Vere-Jones, D. (2005). Diagnostic analysis of space-time branching processes for earthquakes. Lecture Note in Statistics: Case Studies in Spatial Point Process Models (Baddeley, A., Gregori, P., Mateu, J., Stoica, R. and Stoyan, D.), Springer-Verlag, New York, 185, 276–292.

Zhuang, J., Ogata, Y. and Vere-Jones, D. (2002). Stochastic declustering of space-time earthquake occurrences. Journal of the American Statistical Association, 97, 369–380.

Examples


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

  param <- c(0.46, 0.23, 0.022, 2.8, 1.12, 0.012, 2.4, 0.35)

  ## Not run: 
  lambda(15706, 40.12, 34.5, param, iran.cat)
## 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/lambda.Rd_%03d_medium.png", width=480, height=480)
> ### Name: lambda
> ### Title: Clustering Part of Conditional Intensity Function of the ETAS
> ###   Model
> ### Aliases: lambda
> ### Keywords: spatial math earthquake modeling
> 
> ### ** Examples
> 
> 
>   iran.cat <- catalog(iran.quakes, time.begin="1973/01/01",
+      study.start="1996/01/01", study.end="2016/01/01",
+      lat.range=c(25, 42), long.range=c(42, 63), mag.threshold=4.5)
> 
>   param <- c(0.46, 0.23, 0.022, 2.8, 1.12, 0.012, 2.4, 0.35)
> 
>   ## Not run: 
> ##D   lambda(15706, 40.12, 34.5, param, iran.cat)
> ## End(Not run)
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>