astrochron (version 1.5)

timeOptTemplatePlot: 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

if (FALSE) {
# 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))
}