SAPP (version 1.0.4)

pgraph: Graphical Outputs for The Point Process Data Set

Description

Provide the several graphical outputs for the point process data set.

Usage

pgraph(data, mag, threshold=0.0, h, npoint, days, delta=0.0, dmax=0.0, separate.graphics=FALSE)

Arguments

data

point process data.

mag

magnitude.

threshold

threshold magnitude.

h

time length of the moving interval in which points are counted to show the graph.

npoint

number of subintervals in (0,days) to estimate a non parametric intensity under the palm probability measure.

days

length of interval to display the intensity estimate under the palm probability.

delta

length of a subinterval unit in (0,dmax) to compute the variance time curve.

dmax

time length of a interval to display the variance time curve; this is less than (length of whole interval)/4. As the default setting of either delta=0.0 or dmax=0.0, set dmax = (length of whole interval)/4 and delta = dmax/100.

separate.graphics

logical. If TRUE a graphic device is opened for each graphics display.

Value

cnum

cumulative numbers of events time.

lintv

interval length.

tau

=time*(total number of events)/(time end).

nevent

number of events in [tau, tau+h].

survivor

log survivor curve with i*(standard error), i=1,2,3.

deviation

deviation of survivor function from the Poisson.

nomal.cnum

normalized cumulative number.

nomal.lintv

U(i)=-exp(-(normalized interval length)).

success.intv

successive pair of intervals.

occur

occurrence rate.

time

time assuming the stationary Poisson process.

variance

Var(N(0,time)).

error

the 0.95 and 0.99 error lines assuming the stationary Poisson process.

References

Y.Ogata, K.Katsura and J.Zhuang (2006) Computer Science Monographs, No.32, TIMSAC84: STATISTICAL ANALYSIS OF SERIES OF EVENTS (TIMSAC84-SASE) VERSION 2. The Institute of Statistical Mathematics.

Y.Ogata (2006) Computer Science Monographs, No.33, Statistical Analysis of Seismicity - updated version (SASeies2006). The Institute of Statistical Mathematics.

Y.Ogata and K.Shimazaki (1984) Transition from aftershock to normal activity: the 1965 rat islands earthquake aftershock sequence. Bulletin of the seismological society of america, vol. 74, no. 5, pp. 1757-1765.

Examples

  ## The aftershock data of 26th July 2003 earthquake of M6.2
  data(main2003JUL26)
  x <- main2003JUL26
  pgraph(data=x$time, mag=x$magnitude, h=6, npoint=100, days=10)

  ## The residual point process data of 26th July 2003 earthquake of M6.2
  data(res2003JUL26)
  y <- res2003JUL26
  pgraph(data=y$trans.time, mag=y$magnitude, h=6, npoint=100, days=10)