proxy_plot {ananke} | R Documentation |
Plot Layer-Counted Proxy Records Uncertainties
Description
Plot Layer-Counted Proxy Records Uncertainties
Usage
## S4 method for signature 'ProxyRecord,missing'
plot(
x,
calendar = get_calendar(),
iqr = TRUE,
xlab = NULL,
ylab = NULL,
col = grDevices::hcl.colors(12, "YlOrRd", rev = TRUE),
col.mean = "black",
col.iqr = col.mean,
lty.mean = 1,
lty.iqr = 3,
lwd.mean = 2,
lwd.iqr = lwd.mean,
...
)
Arguments
x |
A |
calendar |
An |
iqr |
A |
xlab , ylab |
A |
col |
A list of colors such as that generated by |
col.mean , col.iqr |
A specification for the line colors. Only used if
|
lty.mean , lty.iqr |
A specification for the line types. Only used if
|
lwd.mean , lwd.iqr |
A specification for the line widths. Only used if
|
... |
Further parameters to be passed to |
Value
plot()
is called it for its side-effects: it results in a graphic
being displayed. Invisibly returns x
.
Author(s)
N. Frerebeau
See Also
Other proxy tools:
proxy_ensemble()
Examples
## Get NGRIP records
data("ngrip2010", package = "folio")
ngrip2010 <- subset(ngrip2010, !is.na(MCE))
ngrip2010 <- ngrip2010[nrow(ngrip2010):1, ] # Sort in chronological order
## Replicate fig. 3d from Boers et al. (2017)
## /!\ This may take a while... /!\
ngrip_record <- proxy_ensemble(
positions = ngrip2010$depth,
proxy_values = ngrip2010$delta,
proxy_errors = 0.01,
proxy_step = 0.01,
time_values = ngrip2010$age,
time_errors = ngrip2010$MCE,
calendar = b2k(), # /!\
by = 20,
n = 30
)
plot(ngrip_record)
[Package ananke version 0.1.0 Index]