timeOptTemplatePlot: TimeOptTemplatePlot: Generate summary figure for...

In astrochron: A Computational Tool for Astrochronology

TimeOptTemplatePlot: Generate summary figure for TimeOptTemplate analyses

Description

TimeOptTemplatePlot: Generate summary figure for TimeOptTemplate analyses.

Usage

timeOptTemplatePlot(dat=NULL,template=NULL,detrend=T,detrendTemplate=F,flipTemplate=F,
   srMin=NULL,srMax=NULL,res1=NULL,simres=NULL,fit=1,flow=NULL,fhigh=NULL,roll=NULL,
   targetE=NULL,targetP=NULL,xlab="Depth (m)",ylab="Proxy Value",fitR=NULL,output=0,
   verbose=T)

Arguments

dat	Stratigraphic series used for astrochronologic assessment. First column should be depth or height (in meters), second column should be data value.
template	Instantaneous sedimentation rate template to fit. This represents a unitless proportional sedimentation rate history. Default model is a copy of dat, which will be scaled for instantaneous accumulation optimization.
detrend	Remove linear trend from data series? (T or F)
detrendTemplate	Remove linear trend from sedimentation rate template? (T or F)
flipTemplate	Flip direction of sedimentation rate template? (T or F)
srMin	Minimum sedimentation rate for template
srMax	Maximum sedimentation rate for template
res1	Data frame containing TimeOpt results: sedimentation rate grid, and one of the following - r^2_envelope, r^2_power, or r^2_opt.
simres	Data frame containing the r^2_opt value for each Monte Carlo simulation.
fit	Test for (1) precession amplitude modulation or (2) short eccentricity amplitude modulation?
flow	Low frequency cut-off for Taner bandpass (half power point; in cycles/ka).
fhigh	High frequency cut-off for Taner bandpass (half power point; in cycles/ka).
roll	Taner filter roll-off rate, in dB/octave.
targetE	A vector of eccentricity periods to evaluate (in ka). These must be in order of decreasing period, with a first value of 405 ka.
targetP	A vector of precession periods to evaluate (in ka). These must be in order of decreasing period.
xlab	Label for the depth/height axis.
ylab	Label for proxy variable evaluated.
fitR	The r2 value at the optimal sedimentation rate.
output	Which results you like to return to console? (0) no output; (1) return optimal time series, bandpassed series, Hilbert and fitted periods; (2) return space-time map; (3) return the optimal sedimentation rate at each depth or height
verbose	Verbose output? (T or F)

References

S.R. Meyers, 2015, The evaluation of eccentricity-related amplitude modulations and bundling in paleoclimate data: An inverse approach for astrochronologic testing and time scale optimization: Paleoceanography, v.30, 1625-1640.

S.R. Meyers, 2019, Cyclostratigraphy and the problem of astrochronologic testing: Earth-Science Reviews v. 190, 190-223.

See Also

asm, eAsmTrack, testPrecession, timeOpt, timeOptSim, and timeOptTemplate

Examples

## Not run: 
# EXAMPLE: Differential accumulation across bedding couplets
ex=etp(tmin=0,tmax=1000,dt=5,pWt=1,oWt=0,eWt=1,esinw=TRUE)
ex2=diffAccum(ex,0.01,.05)
ex2=linterp(ex2)

# evaluate with the timeOptTemplate approach
res1=timeOptTemplate(ex2,sedmin=1,sedmax=4,difmin=.5,difmax=6,numsed=100,output=1,ncores=2)
simres=timeOptTemplateSim(ex2,sedmin=1,sedmax=4,difmin=.5,difmax=6,numsed=100,numsim=1000,
 output=2,ncores=2)

timeOptTemplatePlot(ex2,template=ex2,srMin=1,srMax=4,res1=res1,simres=simres,
  targetP=c(23.62069,22.31868,19.06768,18.91979),
  targetE=c(405.6795,130.719,123.839,98.86307,94.87666))

## End(Not run)

astrochron documentation built on June 8, 2025, 10:46 a.m.