SIplot: Plot Global Storm Activity Index WISA

In iWISA: Wavelet-Based Index of Storm Activity

Description

Plot WISA, wavelet index of storm activity, by using ggplot2

Usage

SIplot(x, type = NULL, main = NULL, xlab = NULL, ylab = NULL,
  start = NULL, end = NULL, n.stations = NULL, ...)

Arguments

x	estimation of WISA from SAIndex
type	line width (default=1)
main,	xlab, ylab graphical arguments, see par
xlab	x-axis label
ylab	y-axis label
start	start date of records for magnetic activities.See examples.
end	end date of records for magnetic activities.See examples
n.stations	number of stations included in the study
...	additional arguments

Details

This is an interface for visualizing WISA with ggplot2. SIplot uses as.POSIXct to convert the index object into a data frame with datetime. It essentially uses the mapping geom_line and adds aes(x, y) to visualize the series.

See Also

SAIndex, preindexplot

Examples

## Not run: 
## sample records for one week period
data(record)
coord=matrix(c(124.43, 19.23, 53.77,140.18,68.68,202.00,71.89,293.85),nrow=2,ncol=4)

index.sample<- SAIndex(record, coord, wf="la8")

## example dates
start.date="2001-3-1"
end.date="2001-4-30"
n.station=4
station.names=c("HER","KAK","HON","SJG")

## plot SI
SIplot(index.sample$SI, type=1, start=start.date, end=end.date, 
main="WISA estimation", xlab="Datetime", ylab="iWISA Estimation")

## End(Not run)

iWISA documentation built on May 1, 2019, 7:37 p.m.