Plot M8 Series

Description

Plot the seven M8 series plus the eighth critical series.

Usage

## S3 method for class 'M8':
plot(x, title="M8 Series and TIPs", ...)

Arguments

Details

This produces a 4x2 grid of graphs. The first seven are graphs of the seven series, and the eighth is the critical series. When the training mode is not "moving", the 90th and 75th percentiles for the first 6 and 7th series are plotted, respectively. These percentiles represent the thresholds for each of the seven series. When the training mode is "moving" the threshold values are not constant over time, and are not plotted.

The plot of the “critical series” has a reference line marked at zero. A TIP is called if this series exceeds zero in two consecutive time periods. Points are marked on all plots where a TIP is in effect, the symbols having the following interpretations:

• Successful TIP (STIP): pch = 5
• Near Successful TIP (STIP-): pch = 2
• False TIP (FTIP): pch = 4
• Continuing TIP (CTIP): pch = 3
• TIP caused by quake of target magnitude (c.e.): pch = 0

Note that the plotting symbols may be different on different plotting devices, see the argument pch in the graphics function par.

Value

NULL

Author(s)

Li Dongfeng, 1997; modified by Maaike Vreede, 1998.

See Also

M8, M8.series, M8.TIP

Examples

#   Requires the package ssNZ
require(ssNZ)

as.catalogue(subsetrect(NZ, minday=julian(1,1,1965),
                        maxday=julian(1,1,2000), minmag=4.5),
             catname="NZ45")

decluster.M8(NZ45, cutoff.mag=4.5, decluster.name="NZ.mainshocks")

temp <- M8(NZ.mainshocks, M0=7.0, minday=julian(y=1965, x=1, d=1),
           start.series=julian(y=1975, x=1, d=1), centrelong=176, 
           centrelat=-39, end.training=julian(y=2000, x=1, d=1),
           plotit=FALSE)

savpar <- par(no.readonly=TRUE)
plot(temp, title="M8 Analysis of NZ events")
par(savpar)