The data are classified with the "epoch"
variable,
which separates the series into pre-seismic and seismic
data; this is defined relative to the predicted P-wave
arrival time from a traveltime model.
The original dataset contained NA
values, which were
imputed using zoo::na.locf
, which fills NA
with the last previous observation.
pspectrum
, hfsnm
,
magnet
TauP.R
for an R-implementation of the traveltime
calculations: