astrochron (version 0.3.1)

trackFreq: Frequency-domain minimal tuning: Use interactive graphical interface and sorting to track frequency drift

Description

Frequency-domain minimal tuning: Use interactive graphical interface and sorting algorithm to track frequency drift.

Usage

trackFreq(spec,threshold=NULL,pick=T,fmin=NULL,fmax=NULL,dmin=NULL,dmax=NULL,xmin=NULL,
          xmax=NULL,ymin=NULL,ymax=NULL,h=6,w=4,ydir=1,ncolors=100,genplot=T,verbose=T)

Arguments

spec

Time-frequency spectral results to evaluate. Must have the following format: column 1=frequency; remaining columns (2 to n)=power, amplitude or probability; titles for columns 2 to n must be the location (depth or height). Note that this format is o

threshold

Threshold level for filtering peaks. By default all peak maxima reported.

pick

Pick the peaks of interest using a graphical interface? (T or F). Only activated if genplot=T.

fmin

Minimum frequency for analysis.

fmax

Maximum frequency for analysis.

dmin

Minimum depth/height for analysis. NOT ACTIVATED YET!

dmax

Maximum depth/height for analysis. NOT ACTIVATED YET!

xmin

Minimum frequency for PLOTTING.

xmax

Maximum frequency for PLOTTING.

ymin

Minimum depth/height for PLOTTING.

ymax

Maximum depth/height for PLOTTING.

h

Height of plot in inches.

w

Width of plot in inches.

ydir

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

ncolors

Number of colors to use in plot.

genplot

Generate summary plots? (T or F)

verbose

Verbose output? (T or F)

See Also

eha and traceFreq

Examples

# Generate example series with 3 terms using function 'cycles'.
# Then convert from time to space with sedimentation rate that increases from 1 to 5 cm/ka, using
# function 'sedramp'.
# Finally interpolate to median sampling interval using function 'linterp'.
dat=linterp(sedRamp(cycles(freqs=c(1/100,1/40,1/20),start=1,end=2500,dt=5)))

# EHA anlaysis, output probability results
out=eha(dat,output=4)

## Isolate peaks with probability >= 0.9
#freq=trackFreq(out,0.9)