plt.MTM0 {RSEIS} | R Documentation |
Plot MTM structure
plt.MTM0(frange, prange, plxy, M, freqs, amp, a, dof = dof, Fv = Fv, COL = 2)
frange |
frequency range |
prange |
point range |
plxy |
log x,y axes |
M |
structure from MTM |
freqs |
frequencies |
amp |
amplitude |
a |
list(y=original data, dt=deltat) |
dof |
degrees of freedom |
Fv |
F-values |
COL |
color |
Graphical Side Effect
Jonathan M. Lees<jonathan.lees.edu>
Lees, J. M. and Park, J., 1995: Multiple-taper spectral analysis: A stand-alone C-subroutine, Computers & Geology, 21(2), 199-236.
MTM.drive
data(CE1) plot(CE1$x, CE1$y, type='l') len = length(CE1$y) len2 = 2*next2(len) Mspec = mtapspec(CE1$y, CE1$dt, klen=len2, MTP=list(kind=1,nwin=5, npi=3,inorm=0) ) f=Mspec$freq M = 1 f1 = 0.01 f2 = 100 plxy = '' flag = f>=f1 & f <= f2; freqs = f[flag] mydof = NULL myFv = NULL amp = Mspec$spec[1:length(f)] amp = list(amp[flag]) a = list(y=CE1$y, dt=CE1$dt) frange = range(freqs, na.rm = TRUE) prange = range(amp , na.rm = TRUE) ### plot(freqs, amp[[1]]) plt.MTM0(frange, prange, plxy, M, freqs, amp, a, dof=mydof, Fv=myFv, COL=4)