astrochron (version 0.3.1)

eha: Evolutive Harmonic Analysis

Description

Evolutive Harmonic Analysis using the Thomson Multitaper Method

Usage

eha(dat,tbw=2,pad,fmin,fmax,step,win,demean=T,detrend=T,siglevel=0.90,
    sigID=F,ydir=1,output=0,pl=1,xlab,ylab,genplot=2,verbose=T)

Arguments

dat

Stratigraphic series to analyze. First column should be location (e.g., depth), second column should be data value.

tbw

MTM time-bandwidth product (

pad

Pad with zeros to how many points? Must not factor into a prime number >23. Maximum number of points is 200,000.

fmin

Smallest frequency for analysis and plotting.

fmax

Largest frequency for analysis and plotting.

step

Step size for EHA window, in units of space or time.

win

Window size for EHA, in units of space or time.

demean

Remove mean from data series? (T or F)

detrend

Remove linear trend from data series? (T or F)

siglevel

Significance level for peak identification/filtering (0-1)

sigID

Identify signficant frequencies on power, amplitude, and probabilty plots. Only applies when one spectrum is calculated. (T or F)

ydir

Direction for y-axis in EHA plots (depth,height,time). -1 = values increase downwards (slower plotting), 1 = values increase upwards

output

Return output as new data frame? 0=no; 1=all results; 2=power; 3=amplitude; 4=probability; 5=significant frequencies (only for one spectrum); 6=significant frequencies and their probabilities (only for one spectrum)

pl

Plot logarithm of spectral power (1) or linear spectral power (2)?

xlab

Label for x-axis. Default = "Frequency"

ylab

Label for y-axis. Default = "Location"

genplot

Plotting options. 0= no plots; 1= power, amplitude, f-test, probability; 2=data series, power, amplitude, probability; 3= data series, power, normalized amplitude, filtered (at siglevel) and normalized amplitude

verbose

Verbose output? (T or F)

See Also

extract, trackFreq and traceFreq

Examples

## as an example, evaluate the modelA
data(modelA)

## interpolate to even sampling interval of 0.075 m
ex1=linterp(modelA, dt=0.075)
  
## perform EHA with a time-bandwidth parameter of 2, using an 7.95 meter window, 0.15 m step, 
## and pad to 1000 points
## set labels for plots (optional)
eha(ex1,tbw=2,win=7.95,step=0.15,pad=1000,xlab="Frequency (cycles/m)",ylab="Height (m)")

## for comparison generate spectrum for entire record, using time-bandwidth parameter of 3, and 
## pad to 5000 points
## start by making a new plot
pl(1)
eha(ex1,tbw=3,win=38,pad=5000,xlab="Frequency (cycles/m)")