spect {Rsac} | R Documentation |
This is essentially a wrapper function for spectrum
from
the stats package. Thus, it is important the user be
familiar with the documentation of the functions listed in the
‘See Also’ section. This is simply a quick way to compute the
spectra of "rsac"
class objects.
spect(x, plot = TRUE, main = NULL, ...)
x |
an object of class "rsac" . |
plot |
logical for if the spectrum should be plotted. |
main |
optional character vector for the plots' titles. If not
provided, the code creates the character vector from the station,
component, and network in the header files, but these are not always
reliable. The length of the vector must equal the length of
x . |
... |
further parameters. See spectrum and
plot.spec . |
a vector of objects, each object is of "rsac"
.
Eric M. Thompson <eric.thompson@tufts.edu>
spectrum
, plot.spec
,
spec.pgram
, spec.ar
data(srosa) data(sumatra) eqs <- c(srosa[1], sumatra[11]) spect(eqs) spect(eqs, spans = 101) spect(eqs, spans = 201) spect(eqs, spans = c(101, 151)) spect(eqs, spans = c(101, 151), log = "no") getdelta <- function(X){ delta <- gcp(s = list(lat = X$stla, lon = X$stlo), e = list(lat = X$evla, lon = X$evlo))$a} del <- sapply(X = sumatra, getdelta) sum <- sumatra[c(1, 3, 4, 9)] spect(sum, plot = TRUE, spans = c(101, 151), log = "no", xlim = c(0, 0.15))