Description Usage Arguments Author(s) References See Also Examples

This function plots wavelet power averages across time of a single time
series, which are provided by an object of class `"analyze.wavelet"`

,
or alternatively of class `"analyze.coherency"`

.
(In the latter case, the series number or name must be specified.)
The vertical axis shows the Fourier periods. The horizontal axis shows the averages.
User-defined minimum and maximum average levels can be applied.
Also, an option is given to individualize the period axis and/or axis of averages
by specifying tick marks and labels.

There is an option to label periods according to significance of averages (if p-values are provided) at given levels of significance. Labels are point symbols along the line of averages which can be assigned individually.

The idea to show significance levels by colors of plotting characters and its
implementation has been adopted from Huidong Tian and Bernard Cazelles
(archived R package `WaveletCo`

).

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ```
wt.avg(WT, my.series = 1, exponent = 1,
show.siglvl = TRUE, siglvl = c(0.05, 0.1),
sigcol = c("red", "blue"), sigpch = 20, sigcex = 1,
minimum.level = NULL, maximum.level = NULL,
label.avg.axis = TRUE,
averagelab = NULL, averagetck = 0.02, averagetcl = 0.5,
spec.avg.axis = list(at = NULL, labels = TRUE,
las = 1, hadj = NA, padj = NA),
label.period.axis = TRUE,
periodlab = NULL, periodtck = 0.02, periodtcl = 0.5,
spec.period.axis = list(at = NULL, labels = TRUE,
las = 1, hadj = NA, padj = NA),
show.legend = TRUE, legend.coords = "topright",
main = NULL,
lwd = 1, col = 1,
lwd.axis = 1,
verbose = FALSE)
``` |

`WT` |
an object of class |

`my.series` |
In case Default: |

`exponent` |
Exponent applied to averages before plotting; the exponent should be positive. Default: |

`show.siglvl` |
Label periods according to significance of averages? (Effective only if p-values are provided.) Default: |

`siglvl` |
a vector of significance levels (of any length and order). Default: |

`sigcol` |
a vector of colors (should be of same length as and correspond to Default: |

`sigpch` |
a vector of plotting "characters" (symbols) to use as labels of significance.
(It should be of same length as and correspond to Default: |

`sigcex` |
a numerical vector working as size of labels of significance.
(It should be of same length as and correspond to Default: |

`minimum.level` |
Minimum plot level of wavelet power averages considered. Default: NULL (referring to minimum level observed). |

`maximum.level` |
Maximum plot level of wavelet power averages considered. Default: NULL (referring to maximum level observed). |

`label.avg.axis` |
Label the axis of averages? Logical. Default: |

`averagelab` |
Label for the axis of averages. Default: |

`averagetck` |
length of tick marks on the axis of averages as a fraction of the smaller of the width or height of
the plotting region; see Default here: |

`averagetcl` |
length of tick marks on the axis of averages as a fraction of the height of a line of text; see Default here: |

`spec.avg.axis` |
a list of tick mark and label specifications for individualized labeling of the axis of averages
(only effective if [ `at` :] locations of tick marks (when`NULL` , default plotting will be applied). Valid tick marks can be provided as numerical and non-negative values only. Default:`NULL` .[ `labels` :] either a logical value specifying whether annotations at the tick marks are the tick marks themselves, or any vector of labels. If`labels` is non-logical,`at` should be of same length. Default:`TRUE` .[ `las` :] the style of axis labels, see`par` . Default:`1` (always horizontal).[ `hadj` :] adjustment of labels horizontal to the reading direction, see`axis` . Default:`NA` (centering is used).[ `padj` :] adjustment of labels perpendicular to the reading direction (this can be a vector of adjustments for each label), see`axis` . Default:`NA` (centering is used).
Mismatches will result in a reset to default plotting. |

`label.period.axis` |
Label the (Fourier) period axis? Logical. Default: |

`periodlab` |
(Fourier) period axis label. Default: |

`periodtck` |
length of tick marks on the period axis as a fraction of the smaller of the width or height
of the plotting region; see Default here: |

`periodtcl` |
length of tick marks on the period axis as a fraction of the height of a line of text; see Default here: |

`spec.period.axis` |
a list of tick mark and label specifications for individualized period axis labeling
(only effective if [ `at` :] locations of tick marks (when`NULL` , default plotting will be applied). Valid tick marks can be provided as numerical and positive values only. Default:`NULL` .[ `labels` :] either a logical value specifying whether annotations at the tick marks are the tick marks themselves, or any vector of labels. If`labels` is non-logical,`at` should be of same length. Default:`TRUE` .[ `las` :] the style of axis labels, see`par` . Default:`1` (always horizontal).[ `hadj` :] adjustment of labels horizontal to the reading direction, see`axis` . Default:`NA` (centering is used).[ `padj` :] adjustment of labels perpendicular to the reading direction (this can be a vector of adjustments for each label), see`axis` . Default:`NA` (centering is used).
Mismatches will result in a reset to default plotting. |

`show.legend` |
Include legend of significance levels into the plot? Logical. Default: |

`legend.coords` |
coordinates to position the legend (as in function Default: |

`main` |
an overall title for the plot. Default: |

`lwd` |
width of line of averages. Default: |

`col` |
color of line of averages. Default: |

`lwd.axis` |
line width of axes. Default: |

`verbose` |
Print verbose output on the screen? Logical. Default: |

Angi Roesch and Harald Schmidbauer; credits are also due to Huidong Tian and Bernard Cazelles

Aguiar-Conraria L., and Soares M.J., 2011. The Continuous Wavelet Transform: A Primer. NIPE Working Paper Series 16/2011.

Carmona R., Hwang W.-L., and Torresani B., 1998. Practical Time Frequency Analysis. Gabor and Wavelet Transforms with an Implementation in S. Academic Press, San Diego.

Cazelles B., Chavez M., Berteaux, D., Menard F., Vik J.O., Jenouvrier S., and Stenseth N.C., 2008. Wavelet analysis of ecological time series. Oecologia 156, 287–304.

Liu Y., Liang X.S., and Weisberg R.H., 2007. Rectification of the Bias in the Wavelet Power Spectrum. Journal of Atmospheric and Oceanic Technology 24, 2093–2102.

Tian, H., and Cazelles, B., 2012. `WaveletCo`

.
Available at https://cran.r-project.org/src/contrib/Archive/WaveletCo/, archived April 2013; accessed July 26, 2013.

Torrence C., and Compo G.P., 1998. A practical guide to wavelet analysis. Bulletin of the American Meteorological Society 79 (1), 61–78.

`analyze.wavelet`

, `wt.image`

, `wt.sel.phases`

, `wt.phase.image`

, `reconstruct`

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 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 | ```
## Not run:
## The following example is adopted from Liu et al., 2007:
series.length <- 6*128*24
x1 <- periodic.series(start.period = 1*24, length = series.length)
x2 <- periodic.series(start.period = 8*24, length = series.length)
x3 <- periodic.series(start.period = 32*24, length = series.length)
x4 <- periodic.series(start.period = 128*24, length = series.length)
x <- x1 + x2 + x3 + x4
plot(x, type = "l", xlab = "index", ylab = "", xaxs = "i",
main = "hourly series with periods of 1, 8, 32, 128 days")
my.data <- data.frame(x = x)
## Computation of wavelet power:
## a natural choice of 'dt' in the case of hourly data is 'dt = 1/24',
## resulting in one time unit equaling one day.
## This is also the time unit in which periods are measured.
my.wt <- analyze.wavelet(my.data, "x", loess.span = 0,
dt = 1/24, dj = 1/20,
lowerPeriod = 1/4,
make.pval = TRUE, n.sim = 10)
## Plot of wavelet power spectrum (with equidistant color breakpoints):
wt.image(my.wt, color.key = "i", main = "wavelet power spectrum",
legend.params = list(lab = "wavelet power levels (equidistant levels)"),
periodlab = "period (days)")
## Note:
## The default time axis shows an index of given points in time,
## which is the count of hours in our example.
## With time elapsed in days
## (starting from 0 and proceeding in steps of 50 days)
## instead of the (default) time index:
index.ticks <- seq(1, series.length, by = 50*24)
index.labels <- (index.ticks-1)/24
## Insert your specification of time axis:
wt.image(my.wt, color.key = "i", main = "wavelet power spectrum",
legend.params = list(lab = "wavelet power levels (equidistant levels)"),
periodlab = "period (days)", timelab = "time elapsed (days)",
spec.time.axis = list(at = index.ticks, labels = index.labels))
## Plot of average wavelet power:
wt.avg(my.wt, siglvl = 0.05, sigcol = "red", periodlab = "period (days)")
## The same plot, but with enhanced symbol size, user-defined period axis,
## and horizontal grid:
wt.avg(my.wt, siglvl = 0.05, sigcol = "red", sigcex = 1.3,
periodlab = "period (days)",
spec.period.axis = list(at = c(1,8,32,128)),
periodtck = 1, periodtcl = NULL,
lwd = 1.5, lwd.axis = 0.25)
## Another style of the plot:
## With user-defined period axis and axis of averages,
## minimum and maximum plot levels of averages:
op <- par(no.readonly = TRUE)
par(cex.lab = 1.3, cex.axis = 1.1)
wt.avg(my.wt, siglvl = 0.05, sigcol = "red", sigcex = 1.3,
minimum.level = 0, maximum.level = 11,
periodlab = "period (days)",
spec.period.axis = list(at = c(1,8,32,128)),
spec.avg.axis = list(at = 0:10),
lwd = 1.5)
par(op)
## Note:
## 'cex.axis' in 'par' controls for the size of axis tick labels,
## while 'cex.lab' controls for the size of axis and legend labels.
## Scaling by 'cex' would also affect 'sigcex'.
## Please see also the examples in our guide booklet,
## URL http://www.hs-stat.com/projects/WaveletComp/WaveletComp_guided_tour.pdf.
## End(Not run)
``` |

WaveletComp documentation built on May 2, 2019, 6:33 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.