TSSS (version 1.3.4-5)

plot.tvspc: Plot Evolutionary Power Spectra Obtained by Time Varying AR Model

Description

Plot evolutionary power spectra obtained by time varying AR model returned by tvspc.

Usage

# S3 method for tvspc
plot(x, tvv = NULL, dx = 2, dy = 0.25, ...)

Arguments

x

an object of class "tvspc".

tvv

time varying variance as returned by tvvar.

dx

step width for the X axis.

dy

step width for the Y axis.

...

further graphical parameters may also be supplied as arguments.

Examples

# seismic data
data(MYE1F)
v <- tvvar(MYE1F, trend.order = 2, tau2.ini = 6.6e-06, delta = 1.0e-06,
           plot = FALSE )

z <- tvar(v$nordata, trend.order = 2, ar.order = 8, span = 20,
          outlier = c(630, 1026), tau2.ini = 6.6e-06, delta = 1.0e-06,
          plot = FALSE)

spec <- tvspc(z$arcoef, z$sigma2, span = 20, nf = 400)
plot(spec, tvv = v$tvv, dx = 2, dy = 0.10)