R: Plot MTM structure

plt.MTM0 {RSEIS} | R Documentation

## Plot MTM structure

### Description

Plot MTM structure

### Usage

    plt.MTM0(frange, prange, plxy, M, freqs, amp, a, dof = dof, Fv = Fv, COL = 2)

### Arguments

`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

### Value

Graphical Side Effect

### Author(s)

Jonathan M. Lees<jonathan.lees.edu>

### References

Lees, J. M. and Park, J., 1995: Multiple-taper spectral analysis: A stand-alone C-subroutine, Computers \& Geology, 21(2), 199-236.

### See Also

MTM.drive

### Examples

    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 <- list(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[[1]], amp[[1]])

     plt.MTM0(frange, prange, plxy, M, freqs, amp, a,
            dof=mydof, Fv=myFv, COL=4)

[Package RSEIS version 3.7-4 Index]