proxy_ensemble {ananke} | R Documentation |
Layer-Counted Proxy Records Uncertainties
Description
Represents layer-counted proxy records as sequences of probability distributions on absolute, error-free time axes.
Usage
proxy_ensemble(positions, ...)
## S4 method for signature 'numeric'
proxy_ensemble(
positions,
proxy_values,
proxy_errors,
proxy_step,
time_values,
time_errors,
calendar,
from = NULL,
to = NULL,
by = NULL,
n = 30,
progress = getOption("ananke.progress"),
verbose = getOption("ananke.verbose")
)
Arguments
positions |
A positive |
... |
Currently not used. |
proxy_values |
A |
proxy_errors |
A |
proxy_step |
A length-one |
time_values |
A |
time_errors |
A |
calendar |
An |
from |
A length-one |
to |
A length-one |
by |
A length-one |
n |
An |
progress |
A |
verbose |
A |
Value
A ProxyRecord
object.
Note
This function is currently experimental.
Author(s)
N. Frerebeau
References
Boers, N., Goswami, B. & Ghil, M. (2017). A Complete Representation of Uncertainties in Layer-Counted Paleoclimatic Archives. Climate of the Past, 13(9): 1169-1180. doi:10.5194/cp-13-1169-2017.
See Also
Other proxy tools:
proxy_plot
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)