
Last chance! 50% off unlimited learning
Sale ends in
Simulate a combined AR(1) + ETP signal, plot spectrum and confidence levels
ar1etp(etpdat=NULL,nsim=100,rho=0.9,wtAR=1,sig=90,tbw=2,padfac=5,ftest=F,fmax=0.1,
speed=0.5,pl=2,graphfile=0)
Eccentricity, tilt, precession astronmical series. First column = time, second column = ETP. If not entered will use default series from Laskar et al. (2004), spanning 0-1000 kyr.
Number of simulations.
AR(1) coefficient for noise modeling.
Multiplicative factor for AR1 noise (1= eqivalent to ETP variance). If < 0, etp signal will be excluded from the simulations (noise only)
Demarcate what confidence level (percent) on plots?
MTM time-bandwidth product.
Pad with zeros to (padfac*npts) points, where npts is the number of data points.
Include MTM harmonic f-test results? (T or F)
Maximum frequency for plotting.
Set the amount of time to pause before plotting new graph, in seconds.
Plot (1) log frequency-log power or (2) linear frequency-linear power?
Output a pdf or jpg image of each plot? 0 = no, 1 = pdf, 2 = jpeg. If yes, there will be no output to screen. Individual graphic files will be produced for each simluation, for assembling into a movie.
Note: Setting wtAR=1 will provide equal variance contributions from the etp model and the ar1 model. More generally, set wtAR to the square root of the desired variance contribution (wtAR=0.5 will generate an AR1 model with variance that is 25% of the etp model). If you would like to exclusively evaluate the noise (no etp), set wtAR < 0.
Note: You may use the function etp to generate eccentricity-tilt-precession models.
Laskar, J., Robutel, P., Joutel, F., Gastineau, M., Correia, A.C.M., Levrard, B., 2004, A long term numerical solution for the insolation quantities of the Earth: Astron. Astrophys., Volume 428, 261-285.
getLaskar
, and etp
if (FALSE) {
# run simulations using the default settings
ar1etp()
# compare with a second model:
# generate etp model spanning 0-2000 ka, with sampling interval of 5 ka.
ex1=etp(tmin=0,tmax=2000,dt=5)
# run simulations, with rho=-.7, and scaling noise to have 50
ar1etp(etpdat=ex1,rho=0.7,wtAR=sqrt(0.5))
}
Run the code above in your browser using DataLab