RSEIS (version 4.2-4)

autoreg: Auto-Regressive Spectrum Estimate

Description

Auto-Regressive Spectrum Estimate

Usage

autoreg(a, numf = 1024, pord = 100, PLOT = FALSE, f1 = 0.01, f2 = 50)

Value

LIST:

amp

amplitudes

freq

frequencies, Hz

Arguments

a

signal

numf

number of frequency points to calculate

pord

order

PLOT

logical, TRUE=plot

f1

low frequency

f2

high frequency

Author

Jonathan M. Lees<jonathan.lees.edu>

See Also

fft, mtapspec, plt.MTM0

Examples

data(CE1)
Xamp <- CE1$y[CE1$x>5.443754 & CE1$x<5.615951]

ZIM <- autoreg(Xamp , numf=length(Xamp) , pord = 100, PLOT=FALSE,  f1=.01, f2=50)