timeOptPlot {astrochron}

TimeOptPlot: Generate summary figure for TimeOpt analyses

Description

TimeOptPlot: Generate summary figure for TimeOpt analyses.

Usage

timeOptPlot(dat=NULL,res1=NULL,res2=NULL,simres=NULL,fit=1,fitModPwr,flow=NULL,
            fhigh=NULL,roll=NULL,targetE=NULL,targetP=NULL,xlab="Depth (m)",
            ylab="Proxy Value",fitR=NULL,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.

res1

Data frame containing TimeOpt results: sedimentation rate grid, r^2_envelope, r^2_power, r^2_opt.

res2

Data frame containing the optimal-fitted time series, bandpassed series, envelope, and reconstructed eccentricity model.

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?

fitModPwr

Include the modulation periods in the spectral fit? (T or F)

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_opt value at the optimal sedimentation rate.

verbose

Verbose output? (T or F)

References

S.R. Meyers, 2015, The evaluation of eccentricity-related amplitude modulation and bundling in paleoclimate data: An inverse approach for astrochronologic testing and time scale optimization: Paleoceanography, 30, doi:10.1002/2015PA002850.

See Also

asm, eAsmTrack, testPrecession, timeOpt, and timeOptSim

Examples

 ## Not run: 
# generate a test signal with precession and eccentricity
ex=etp(tmin=1,tmax=1000,dt=1,pWt=1,oWt=0,eWt=1,esinw=TRUE,genplot=FALSE,verbose=FALSE)
# convert to meters with sedimentation rate of 2 cm/kyr
ex[1]<-ex[1]*0.02
# evaluate precession modulations
res1=timeOpt(ex,sedmin=0.5,sedmax=5,numsed=100,fit=1,output=1)
res2=timeOpt(ex,sedmin=0.5,sedmax=5,numsed=100,fit=1,output=2)
simres=timeOptSim(ex,sedrate=2,numsim=2000,fit=1,output=2)
timeOptPlot(ex,res1,res2,simres,flow=0.035,fhigh=0.065,roll=10^3,
 targetE=c(405.6795,130.719,123.839,98.86307,94.87666),
 targetP=c(23.62069,22.31868,19.06768,18.91979),xlab="Depth (m)",
 ylab="Value",fitR=0.832,verbose=T)
 
## End(Not run)