trackFreq {astrochron}

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,palette=6,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 ouput by function eha.

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.

palette

What color palette would you like to use? (1) rainbow, (2) grayscale, (3) blue, (4) red, (5) blue-white-red, (6) viridis

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

 ## Not run: 
# Check to see if this is an interactive R session, for compliance with CRAN standards.
# YOU CAN SKIP THE FOLLOWING LINE IF YOU ARE USING AN INTERACTIVE SESSION.
if(interactive()) {

# 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.8
freq=trackFreq(out,0.8)

}
 
## End(Not run)