astrochron (version 0.4.3)

lowpass: Lowpass Filter Stratigraphic Series

Description

Lowpass filter stratigraphic series using rectangular, Gaussian or tapered cosine window [cosine window is experimental]

Usage

lowpass(dat,padfac=2,fcut=NULL,win=0,demean=T,detrend=F,addmean=T,alpha=3,p=0.25,
        xmin=0,xmax=Nyq,genplot=T,verbose=T)

Arguments

dat

Stratigraphic series for lowpass filtering. First column should be location (e.g., depth), second column should be data value.

padfac

Pad with zeros to (padfac*npts) points, where npts is the original number of data points.

fcut

Cutoff frequency for lowpass filtering.

win

Window type for bandpass filter: 0 = rectangular , 1= Gaussian, 2= Cosine-tapered window.

demean

Remove mean from data series? (T or F)

detrend

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

addmean

Add mean value to bandpass result? (T or F)

alpha

Gaussian window parameter: alpha is 1/stdev, a measure of the width of the Dirichlet kernal. Larger values decrease the width of data window, reduce discontinuities, and increase width of the transform. Choose alpha >= 2.5.

p

Cosine-tapered window parameter: p is the percent of the data series tapered (choose 0-1).

xmin

Smallest frequency for plotting.

xmax

Largest frequency for plotting.

genplot

Generate summary plots? (T or F)

verbose

Verbose output? (T or F)

See Also

bandpass, noKernel, noLow, prewhiteAR, and prewhiteAR1

Examples

# generate example series with periods of 405 ka, 100 ka and 20 ka, plus noise
ex=cycles(freqs=c(1/405,1/100,1/20),noisevar=.1,dt=5)
# lowpass filter using cosine-tapered window
res_ex=lowpass(ex,fcut=.02,win=2,p=.4)