timeOptSim {astrochron} | R Documentation |
Monte Carlo simulation for TimeOpt
Description
Perform Monte Carlo AR1 simulations to evaluate significance of TimeOpt results, as in Meyers (2015).
Usage
timeOptSim(dat,numsim=2000,rho=NULL,sedrate=NULL,sedmin=0.5,sedmax=5,numsed=100,
linLog=1,limit=T,fit=1,r2max=1,fitModPwr=T,flow=NULL,fhigh=NULL,roll=NULL,
targetE=NULL,targetP=NULL,detrend=T,ncores=2,output=0,genplot=T,
check=T,verbose=T)
Arguments
dat |
Stratigraphic series for astrochronologic assessment. First column should be depth or height (in meters), second column should be data value. |
numsim |
Number of Monte Carlo AR1 simulations. |
rho |
AR1 coefficient to use in simulations. By default this will be estimated from the stratigraphic series. |
sedrate |
Sedimentation rate for investigation (cm/ka). This option is for compatibility with prior versions of timeOptSim. Please use sedmin, sedmax, numsed. |
sedmin |
Minimum sedimentation rate for investigation (cm/ka). |
sedmax |
Maximum sedimentation rate for investigation (cm/ka). |
numsed |
Number of sedimentation rates to investigate in optimization grid. |
linLog |
Use linear or logarithmic scaling for sedimentation rate grid spacing? (0=linear, 1=log) |
limit |
Limit evaluated sedimentation rates to region in which full target signal can be recovered? (T or F) |
fit |
Test for (1) precession amplitude modulation or (2) short eccentricity amplitude modulation? |
r2max |
Which maximum in r2 do you want to use for the time model calibration? (1) r^2_opt, (2) r^2_spectral, (3) r^2_envelope |
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. |
detrend |
Remove linear trend from data series? (T or F) |
ncores |
Number of cores to use for parallel processing. Must be >=2 |
output |
Which results would you like to return to console? (0) no output; (1) p-value; (2) simulation r2 results |
genplot |
Generate summary plots? (T or F) |
check |
Conduct compliance checks before processing? (T or F). In general this should be activated. |
verbose |
Verbose output? (T or F) |
Details
TimeOpt is an astronomical testing algorithm for untuned (spatial) stratigraphic data. The algorithm identifies the sedimentation rate(s) that simultaneously optimizes: (1) eccentricity amplitude modulations within the precession band, and (2) the concentration of spectral power at specified target astronomical periods. The statistical significance of the r^2_opt is determined via Monte Carlo simulation using timeOptSim.
The present version of timeOptSim improves upon the original significance testing method of Meyers (2015), by conducting simulations across the entire sedimentation grid. This approach more rigorously protects against inflation of the p-value due to multiple testing. Parallel processing has been implemented to address the greater computational demand that is required.
See timeOpt for more information on the basic methodology.
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.
See Also
asm
, eAsm
, eAsmTrack
, testPrecession
, timeOpt
, and timeOptPlot
Examples
## Not run:
# generate a test signal with precession and eccentricity
ex=etp(tmin=1,tmax=1000,dt=5,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 with timeOptSim. be patient, this may take a while to run.
timeOptSim(ex,sedmin=0.5,sedmax=5,numsed=100)
## End(Not run)