Description Usage Arguments Value References Examples

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

1 2 |

`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 nonparametric 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 an interval to display the variance time curve; |

`separate.graphics` |
logical. If |

`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. |

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

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

Ogata, Y. and Shimazaki, K. (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.

1 2 3 4 5 6 7 8 9 | ```
## The aftershock data of 26th July 2003 earthquake of M6.2
data(main2003JUL26)
x <- main2003JUL26
pgraph(x$time, 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(y$trans.time, y$magnitude, h = 6, npoint = 100, days = 10)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.