Description Usage Arguments Details Value Warning Author(s) References See Also Examples
The forward complex-valued lifting transform for decomposing a signal of interest is dependent on the trajectory (lifting order) used in the forward lifting transform. This procedure uses trajectory bootstrapping to provide (complex-valued) time-scale information at all times and scales for bivariate series
1 2 |
x1 |
A vector of grid values. Can be of any length, not necessarily equally spaced. |
x2 |
An optional vector of grid values corresponding to f2. Can be of any length, not necessarily equally spaced. If not specified ( |
f1 |
A vector of function values of the first component of a bivariate series, corresponding to x. Must be of the same length as x. |
f2 |
A vector of function values of the second component of a bivariate series, corresponding to x. Must be of the same length as x. |
P |
Number of trajectories to be used in the nondecimated lifting transform. |
nkeep |
Number of scaling points we want at the end of the transform. The usual choice is |
use.same.trajectories |
A boolean variable indicating whether the same set of trajectories should be used for both components of the bivariate signal. |
verbose |
Indicates whether useful messages should be printed to the console during the procedure. |
... |
Any other arguments to be passed to |
Essentially, this function applies the forward complex wavelet lifting transform fwtnppermC
P times on both (x,f1) and (x,f2), each with a different random lifting
trajectory. This results in P sets of complex-valued detail coefficients, along with their associated scales. This information is stored in order to compute the
cross-periodograms for the bivariate series (x,f1,f2). The “degree of asymmetry" in the prediction is also recorded.
This is the ratio between the maximum distance to the removed point to one of its neighbours and the minimum distance from the removed point to one of its neighbours,
see Chapter 5.3 in Sanderson (2010) for more details.
An object of class cnlt
(subclasses biv
and either SG
or DG
).
If both components have the same grid (subclass SG
), a list with components:
x1 |
The sampling grid corresponding to f1 used for the decomposition. |
x2 |
The sampling grid corresponding to f2 used for the decomposition. If the object is of subclass |
det1 |
A list, entry |
det2 |
A list, entry |
lre |
A list, entry |
lreA |
A list, entry |
If both components have different sampling grids, the additional following list components are returned:
lre2 |
A list, entry |
lreA2 |
A list, entry |
Using a large number of trajectories for long datasets could take a long time!
Matt Nunes, Jean Hamilton
Hamilton, J., Nunes, M. A., Knight, M. I. and Fryzlewicz, P. (2018) Complex-valued wavelet lifting and applications.
Technometrics, 60 (1), 48-60, DOI 10.1080/00401706.2017.1281846.
Sanderson, J. (2010) Wavelet methods for time series with bivariate observations and irregular sampling grids. PhD Thesis, University of Bristol, UK.
For the real-valued equivalent procedure, see also
Knight, M. I., Nunes, M. A. and Nason, G. P. (2012) Spectral Estimation for Locally Stationary Time Series with Missing Observations. Stat. Comput.
22 (4), 877–895.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 | # a bivariate series example with same grids
# simulate data, e.g. two sinusoids
dat <- seq(from=1, to=3, by=0.1)
x1 <- cumsum(sample(dat, 200, TRUE))
y1 <-sin(2*pi*(1/25)*x1) + sin(2*pi*(1/50)*x1)+ 1*sin(2*pi*(1/10)*x1)+ rnorm(length(x1), 0,0.2)
y3 <- c(sin(2*pi*(1/25)*x1[1:100]),sin(2*pi*(1/25)*x1[97:196]))+ rnorm(length(x1), 0,0.1)
## Not run:
y1y3.dec<-cnlt.biv(x1, f1=y1, f2=y3, P = 500)
# the complex detail coefficients corresponding to the first timepoint are:
y1y3.dec$det1[[1]]
## End(Not run)
# a bivariate series example with different grids
# load some data in
data(Baidu)
data(Google)
## Not run:
BaiGoo<-cnlt.biv(Baidu$Seconds[1:100], Google$Seconds[1:100], Baidu$Return[1:100],
Google$Return[1:100], P = 500)
# now look at some of the coefficients from the decomposition
# (the complex detail coefficients corresponding to the first timepoint:
BaiGoo$det1[[1]]
BaiGoo$det2[[1]]
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.