wmr: Compute the wavelet modulus ratio of multivariate data
In mvcwt: Wavelet Analysis of Multiple Time Series

Description Usage Arguments Value Author(s) References See Also Examples

Computes the wavelet modulus ratio described in Keitt (2008). A value of one indicated perfect synchrony among time series and a value of zero, perfect compensation.

1	wmr(w, smoothing = 1)

`w`	an object such as returned by `mvcwt`
`smoothing`	width of smoothing kernel; larger values give more smoothing

an object of class "mvcwt"

Timothy H. Keitt

Keitt, T. H. 2008. Coherent ecological dynamics induced by large-scale disturbance. Nature 454:331-4. doi:10.1038/nature06935.

mvcwt, image.mvcwt

## Not run: 
data(lrlake)
x = subset(lrlake, Basin == "Treatment", LRL.Day) / 365.25
y = subset(lrlake, Basin == "Treatment", -(1:8))
w = mvcwt(x, y, min.scale = 0.25, max.scale = 4)
mr = wmr(w)
image(mr, reset.par = FALSE)
contour(mr, bound = NA, add = TRUE)
## End(Not run)