Plot Global Storm Activity Index WISA

Description

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

Usage

1
2
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
## 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)