asm {astrochron}

Average Spectral Misfit

Description

Calculate Average Spectral Misfit with Monte Carlo spectra simulations, as updated in Meyers et al. (2012).

Usage

asm(freq,target,fper=NULL,rayleigh,nyquist,sedmin=1,sedmax=5,numsed=50,
    linLog=1,iter=100000,output=F,genplot=T)

Arguments

freq

A vector of candidate astronomical cycles observed in your data spectrum (cycles/m). Maximum allowed is 500.

target

A vector of astronomical frequencies to evaluate (1/ka). These must be in order of increasing frequency (e.g., e1,e2,e3,o1,o2,p1,p2). Maximum allowed is 50 frequencies.

fper

A vector of uncertainties on each target frequency (1/ka). Values should be from 0-1, representing uncertainty as a percent of each target frequency. The order of the uncertainties must follow that of the target vector. By default, no uncertainty is assigned.

rayleigh

Rayleigh frequency (cycles/m).

nyquist

Nyquist frequency (cycles/m).

sedmin

Minimum sedimentation rate for investigation (cm/ka).

sedmax

Maximum sedimentation rate for investigation (cm/ka).

numsed

Number of sedimentation rates to investigate in ASM optimization grid. Maximum allowed is 500.

linLog

Use linear or logarithmic scaling for sedimentation rate grid spacing? (0=linear, 1=log)

iter

Number of Monte Carlo simulations for significance testing. Maximum allowed is 100,000.

output

Return output as a new data frame? (T or F)

genplot

Generate summary plots? (T or F)

Details

This function will caculate the Average Spectral Misfit between a data spectrum and astronomical target spectrum, following the approach outlined in Meyers and Sageman (2007), and the improvements of Meyers et al. (2012).

Value

A data frame containing: Sedimentation rate (cm/ka), ASM (cycles/ka), Null hypothesis significance level (0-100 percent), Number of astronomical terms fit.

References

S.R. Meyers and B.B. Sageman, 2007, Quantification of Deep-Time Orbital Forcing by Average Spectral Misfit: American Journal of Science, v. 307, p. 773-792.

S.R. Meyers, B.B. Sageman and M.A. Arthur, 2012, Obliquity forcing of organic matter accumulation during Oceanic Anoxic Event 2: Paleoceanography, 27, PA3212, doi:10.1029/2012PA002286.

See Also

eAsm, eAsmTrack, testPrecession, timeOpt, and timeOptSim

Examples

## These frequencies are from modelA (type '?astrochron' for more information). 
## They are for an 8 meter window, centered at 22 meters height. Units are cycles/m . 
freq <- c(0.1599833,0.5332776,1.5998329,2.6797201,3.2796575,3.8795948,5.5194235,6.5459830)
freq <- data.frame(freq)

## Rayleigh frequency in cycles/m
rayleigh <- 0.1245274

## Nyquist frequency in cycles/m
nyquist <- 6.66597

## orbital target in 1/ky. Predicted periods for 94 Ma (see Meyers et al., 2012)
target <- c(1/405.47,1/126.98,1/96.91,1/37.66,1/22.42,1/18.33)

## percent uncertainty in orbital target
fper=c(0.023,0.046,0.042,0.008,0.035,0.004)

asm(freq=freq,target=target,fper=fper,rayleigh=rayleigh,nyquist=nyquist,sedmin=0.5,sedmax=3,
    numsed=100,linLog=1,iter=100000,output=FALSE)