epicentres.identify {ssEDA}R Documentation

Identify Epicentre Outliers

Description

Plots earthquake epicentres selected from a given catalogue. The points are plotted on the screen, and then may be identified.

Usage

epicentres.identify(events, mapname = "world2", criteria = TRUE)

Arguments

events defines the events to be plotted. It is an object of class "subset", generally created by subsetcircle, subsetpolygon, subsetrect or subsetsphere.
mapname character string giving the map name. Low resolution maps contained in package maps are "nz" and "world2" (default). High resolution maps contained in the package mapdata are "nzHires", "world2Hires" and "china". The “2” on "world2" denotes the version where the longitudes are all positive.
criteria boolean variable. If TRUE (default) the subsetting parameters are printed below the graph.

Details

The points are plotted on an appropriate map, and points which are considered to be outliers are highlighted with larger circles. After the points have been plotted, the identify command is used for the user to identify points of interest. Use the left mouse button to select points (which are then identified on the map with their sequence number within the list of points plotted). Then use the middle or right button to terminate the identification process.

Value

A vector of indices (into the original catalogue) of the points identified.

See Also

epicentres, identify

Examples

data(NZ55)

a <- subsetrect(NZ55, minmag=6.5)
b <- epicentres.identify(a, mapname = "nz")
if (length(b) > 0) print(NZ55[b,])

[Package ssEDA version 2.2-1 Index]