Learn R Programming

remote (version 1.2.1)

plot: Plot an Eot* object

Description

This is the standard plotting routine for the results of eot. Three panels will be drawn i) the predictor domain, ii) the response domain, iii) the time series at the identified base point

Usage

"plot"(x, y, pred.prm = "rsq", resp.prm = "r", show.bp = FALSE, anomalies = TRUE, add.map = TRUE, ts.vec = NULL, arrange = c("wide", "long"), clr = NULL, locations = FALSE, ...)
"plot"(x, y, pred.prm = "rsq", resp.prm = "r", show.bp = FALSE, anomalies = TRUE, add.map = TRUE, ts.vec = NULL, arrange = c("wide", "long"), clr = NULL, locations = FALSE, ...)

Arguments

x
either an object of EotMode or EotStack as returned by eot
y
integer or character of the mode to be plotted (e.g. 2 or "mode_2")
pred.prm
the parameter of the predictor to be plotted. Can be any of "r", "rsq", "rsq.sums", "p", "int" or "slp"
resp.prm
the parameter of the response to be plotted. Can be any of "r", "rsq", "rsq.sums", "p", "int" or "slp"
show.bp
logical. If TRUE a grey circle will be drawn in the predictor image to indicate the location of the base point
anomalies
logical. If TRUE a reference line will be drawn a 0 in the EOT time series
add.map
logical. If TRUE country outlines will be added to the predictor and response images
ts.vec
an (optional) time series vector of the considered EOT calculation to be shown as the x-axis in the time series plot
arrange
whether the final plot should be arranged in "wide" or "long" format
clr
an (optional) color palette for displaying of the predictor and response fields
locations
logical. If x is an EotStack, set this to TRUE to produce a map showing the locations of all modes. Ignored if x is an EotMode
...
further arguments to be passed to spplot

Methods (by class)

  • x = EotStack,y = ANY: EotStack

Examples

Run this code
data(vdendool)

## claculate 2 leading modes
nh_modes <- eot(x = vdendool, y = NULL, n = 2, 
                standardised = FALSE, 
                verbose = TRUE)

## default settings 
plot(nh_modes, y = 1) # is equivalent to

## Not run: 
# plot(nh_modes[[1]]) 
# 
# plot(nh_modes, y = 2) # shows variance explained by mode 2 only
# plot(nh_modes[[2]]) # shows cumulative variance explained by modes 1 & 2
# 
# ## showing the loction of the mode
# plot(nh_modes, y = 1, show.bp = TRUE)
# 
# ## changing parameters
# plot(nh_modes, y = 1, show.bp = TRUE,
#      pred.prm = "r", resp.prm = "p")
#         
# ## change plot arrangement
# plot(nh_modes, y = 1, show.bp = TRUE, arrange = "long") 
# 
# ## plot locations of all base points
# plot(nh_modes, locations = TRUE)
# ## End(Not run)

Run the code above in your browser using DataLab