Fixed interpretation of arrow directions in the documentation
Merged PR from Dirk Eddelbuettel for defining STRICT_R_HEADERS in Rcpp
Updated rcpp_wt_bases_paul.cpp to fix array
out-of-bounds reading issue
Updated plot.biwavelet help file with correct phase arrow interpretation
Updated package by removing benchmarks in vignettes for CRAN submission
yaxis tickmarks should now be accurate (no more rounding issues)
Fixed documentation for wtc function
Fixed return NULL issues with Windows platforms
Fixed plot.biwavelet so that the COI extends all
the way to the bottom of the plot (max of periods)
Fixed plot.biwavelet so that the arrow.cutoff
argument applies to rsq values for wtc and pwtc objects
Fixed plot.biwavelet so that the arrow.cutoff
argument applies to rsq values
Build vignettes
Function phase.biwavelet now plots the regions whose significance
exceeds arrow.cutoff. If the object being plotted does not have
a significance field, regions whose zvalues exceed the arrow.cutoff
quantile will be plotted.
Fixed C++ warning about unsupported dynamically sized arrays
Fixed handling of lag1 coefficients in wtc.
Fixed handling of axis preferences in plot.biwavelet.
Fixed handling of time in check.data.
Fixed x-axis in plot.biwavelet.
Fixed displacement of COI, contours and phase arrows in plot.biwavelet when adding a color bar.
Fixed check.datum; hopefully for the last time.
Faster wt bases and row quantile (Rcpp implementations):
The param parameter for all rcpp_wt_bases_* must be
within interval (0..10).
The rcpp_row_quantile function requires
a matrix as a parameter (use as.matrix() for vectors).
Fixed Rcpp implementation of the wt.bases functions,
i.e., rcpp_wt_bases_*.cpp.
Replacing int type with double type for scale parameter
which caused small scales to be rendered incorrectly.
Fixed interpretation of phase differences in plot.biwavelet help file
Added unit tests for approx 78% of the code.
Implemented a parallelized Monte Carlo simulation function
wtc_sig_parallel which is 2 to 4 times faster on a 4-core CPU
than the original wtc.sig. The speedup is noticeable on:
large simulations nrads >= 800,
multiple simulations,
multi-core systems with 4+ cores.
However, parallelization involves a significant heat-up phase because all the workers need to be started and they need to load all the required packages. This will be addresses in future versions of biwavelet.
Added a speed-optimized version of convolve2D.
Replaced standard arima.sim function with a pair of functions
get_minroots and ar1_ma0_sim. These functions are still
implemented in R. We can reimplement them later in C.
Reimplemented wt.bases morlet, paul and dog in C.
Removed unused function meshgrid.
close all progress bars after use
Function wtc can now handle non-finite values when computing
the quantiles of the rsq values from the Monte Carlo simulations
Added ability to handle custom color palettes in plot.biwavelet.
Users can now specify any color scheme using the fill.cols argument.
Fixed limited padding issue, which could lead to weird edge effects. Current padding level is identical to that of Torrence & Compo (1998).
Changed the default tol value from 0.95 to 1 in the plot function.
Added semi-transparent COI.
Fixed check.datum function so that it does not assume a sampling frequency of 1.
Added ability to set zlim in plot.biwavelet.
Improved the implementation of phase.plot to allow for much better looking phase arrows (thanks Huidong Tang).
Made function wt faster by avoiding excessive padding (thanks Huidong Tang).
Made check.datum tolerate slight inconsistencies in the size of timesteps.
Added arguments in plot.biwavelet and phase.plot to control the length
of the phase arrows and the size of the arrow heads independently.
Fixed code in check.data to test for constant step size in the data.
Function pwtc can be used to perform partial wavelet coherence between two time series
y and x1 by controlling for (or partialling-out) a third time series x2.
Users can now specify the density of the phase arrows using the plot function.
Fixed bug in wt affecting the significance region (thanks Patrick Kilduff and Flora Cordoleani).
Users can now specify the color, line width and line type for
the COI, significance contours and phase arrows using the plot function.
Removed misleading examples showing how to compute the 'bias-corrected' wavelet coherence.
There is no bias for the wavelet coherence function, so using the default type argument in the
plot function is recommended.
Fixed typos in the documentation of plot.biwavelet and xwt (thanks Lei Cheng).
As of biwavelet version 0.14, the bias-corrected wavelet and cross-wavelet spectra are automatically computed and plotted by default using the methods described by Liu et al. (2007) and Veleda et al. (2012). This correction is needed because the traditional approach for computing the power spectrum (e.g., Torrence and Compo 1998) leads to an artificial and systematic reduction in power at lower periods.
Plotting function now accepts traditional plotting flags such as xaxt and yaxt to control x and y tickmarks.