mvcwt (version 1.3.1)

wmr: Compute the wavelet modulus ratio of multivariate data

Description

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.

Usage

wmr(w, smoothing = 1)

Arguments

w

an object such as returned by mvcwt

smoothing

width of smoothing kernel; larger values give more smoothing

Value

an object of class "mvcwt"

References

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

See Also

mvcwt, image.mvcwt

Examples

# 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)
# }