plot.sealevel {oce} | R Documentation |
Plot a summary diagram for sealevel data.
## S4 method for signature 'sealevel' plot(x, which=1:3, adorn=NULL, drawTimeRange=getOption("oceDrawTimeRange"), mgp=getOption("oceMgp"), mar=c(mgp[1]+0.5,mgp[1]+1.5,mgp[2]+1,mgp[2]+3/4), marginsAsImage=FALSE, debug=getOption("oceDebug"), ...)
x |
an object of class |
which |
a numerical or string vector indicating desired plot types, with
possibilities 1 or |
adorn |
vector of expressions to be executed for the panels in turn, e.g. to adorn the plots. If the number matches the number of panels, then the strings are applied to the appropriate panels, as they are drawn from top-left to bottom-right. If only a single expression is provided, it is used for all panels. (See “Examples”.) |
drawTimeRange |
boolean that applies to panels with time as the horizontal axis, indicating whether to draw the time range in the top-left margin of the plot. |
mgp |
3-element numerical vector to use for |
mar |
value to be used with |
marginsAsImage |
boolean, |
debug |
a flag that turns on debugging, if it exceeds 0. |
... |
optional arguments passed to plotting functions. |
Creates a plot for a sea-level dataset, in one of two varieties.
Depending on the length of which
, either a single-panel or
multi-panel plot is drawn. If there is just one panel, then the value of
par
used in plot.sealevel
is retained upon exit, making it
convenient to add to the plot. For multi-panel plots, par
is
returned to the value it had before the call, and so adorn
must be
used to add to individual panels.
None.
Dan Kelley
The example refers to Hurricane Juan, which caused a great deal of damage to Halifax in 2003. Since this was in the era of the digital photo, a casual web search will uncover some spectacular images of damage, from both wind and storm surge. A map of the path of Hurricane Juan across Nova Scotia is at http://ec.gc.ca/ouragans-hurricanes/default.asp?lang=En&n=222F51F7-1 (link checked April 15, 2015). Landfall, very near the site of this sealevel gauge, was between 00:10 and 00:20 Halifax local time on Monday, Sept 29, 2003.
The documentation for sealevel-class
explains the
structure of sealevel objects, and also outlines the other functions
dealing with them.
library(oce) data(sealevel) ## local Halifax time is UTC + 4h; see [2] on timing juan <- as.POSIXct("2003-09-29 00:15:00", tz="UTC")+4*3600 plot(sealevel, which=1, xlim=juan+86400*c(-7,7)) abline(v=juan, col='red')