traceFreq {astrochron}

Frequency-domain minimal tuning: Use interactive graphical interface to trace frequency drift

Description

Frequency-domain minimal tuning: Use interactive graphical interface to trace frequency drift.

Usage

traceFreq(spec,color=2,h=6,w=4,xmin=NULL,xmax=NULL,ymin=NULL,ymax=NULL,ydir=1,
          palette=6,ncolors=100,path=1,pl=0)

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.

color

Line color for tracing. 1 = transparent black; 2 = transparent white; 3 = transparent yellow

h

Height of plot in inches.

w

Width of plot in inches.

xmin

Minimum spatial frequency to plot.

xmax

Maximum spatial frequency to plot.

ymin

Minimum depth/height to plot.

ymax

Maximum depth/height to plot.

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.

path

How do you want to represent the spatial frequency path?: 1=lines and points; 2=lines; 3=points

pl

An option for the color plots: 0=linear scale; 1=plot log of value, 2=normalize to maximum value

See Also

eha and trackFreq

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 amplitude results
out=eha(dat,output=3)

## Interactively track frequency drift
freq=traceFreq(out)

}
 
## End(Not run)