timeOptPlot: TimeOptPlot: Generate summary figure for TimeOpt analyses

In astrochron: A Computational Tool for Astrochronology

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,sedmin=0.5,sedmax=5,numsed=100,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)

astrochron documentation built on Sept. 30, 2024, 9:14 a.m.