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

View source: R/tvvar.R

plot.tvspcR Documentation

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 class '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)

TSSS documentation built on Sept. 29, 2023, 9:07 a.m.