plot.imftest: Graphical IMF significance test

In PierreMasselot/Library--emdr: Tools for EMD-regression

Description

Graphically displays the result of an IMF significance test by plotting the log energy vs. log period of IMFs.

Usage

## S3 method for class 'imftest'
plot(
  x,
  level = 0.05,
  select.var = NULL,
  select.imf = NULL,
  alternative = x$alternative,
  nline = 500,
  col = "black",
  col.signif = NULL,
  pch = 4,
  pch.signif = 8,
  imf.names = NULL,
  col.names = NULL,
  col.lines = "black",
  lty.lines = 2,
  lwd.lines = 1,
  add = FALSE,
  ...
)

Arguments

x	An imftest x resulting from a call to link{imf.test}.
level	The significance level at which to draw the rejectance region lines.
select.var	Character of numeric vector giving a subset of variables to plot.
select.imf	Character of numeric vector giving a subset of IMFs to plot.
alternative	The alternative hypothesis, either "one.sided" or "two.sided".
nline	Numeric value giving the resolution of the rejectance region lines.
col	The colors of the points.
col.signif	The color to use for significant IMFs exclusively. If NULL, col is used.
pch	The plotting character of IMFs.
pch.signif	Specific pch for significant IMFs only.
imf.names	Character vector of IMF labels to plot next to the points.
col.names	Color of the IMF labels.
col.lines	Colors of the significance lines.
lty.lines	Significance lines type.
lwd.lines	Significance lines width.
add	Logical value. If TRUE, the plot is added to an existing one.
...	Additional graphical parameters. See par.

References

Wu, Z., Huang, N.E., 2004. A study of the characteristics of white noise using the empirical mode decomposition method. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 460, 1597–1611.

Examples

   library(dlnm)
   
   # Decompose both temperature and relative humidity with NA-MEMD
   # Adding two noise variables 
   X <- chicagoNMMAPS[,c("temp", "rhum")]
   set.seed(123)
   mimfs <- memd(X, l = 2) # Takes a couple of minutes
   cmimfs <- combine.mimf(mimfs, list(10:11, 12:13), 
     new.names = c("C10", "C11"))
   
   # Apply the test of Wu and Huang (2004)
   testres <- imf.test(cmimfs)
   plot(testres)

PierreMasselot/Library--emdr documentation built on June 19, 2021, 8:58 a.m.